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Das Anwachsen der Grenzschicht in und hinter einer Expansionswelle 


Von E. Becker 


1. Einleitung. Langs einer ebenen Wand, die an anfangs ruhendes Gas grenzt, laufe eine ebene 
Verdiinnungswelle mit wandparalleler Wellennormale (Abb. 1). Die Dichte 9), der Druck p, und 
die Schallgeschwindigkeit a) des Gases vor der Welle fallen beim Durchgang der Welle an jeder 
Stelle auf die kleineren Werte 0,, p,, a, hinter der Welle ab. Gleichzeitig wird das Gas parallel 
zur Wand in der zur Wellenausbreitung entgegengesetzten Richtung mit der Geschwindigkeit U 
in Bewegung gesetzt; U nimmt hinter der Welle einen konstanten Wert — U, an. Beim Fort- 
schreiten flacht sich die Welle laufend ab, Wellenkopf und Wellenende riicken immer weiter aus- 
einander. 

An der ebenen Wand entsteht nach Durchgang des Wellenkopfes eine Grenzschicht, in der die 
Reibungseinfliisse maBgebend sind. In der vorliegenden Arbeit werden fiir das Anwachsen der 
Grenzschichtdicke in laminarer und turbulenter Strémung einfache Naherungsformeln gegeben 
werden, die sich auf den Impulssatz griinden. 


Eine Kenntnis dieses Wachstums ist bei verschie- Oy» Py Oy Qo: Po, 20 
denen Anordnungen von Interesse, z. B. beim U=-U, U=0 
StoBwellenrohr+? und bei dem von H. Ludwieg 
vorgeschlagenen Rohrwindkanal®. AuBerdem a Bey a, 
: : ° 7 dimes, 0 
lassen sich durch einfache Veranderungen in den : — 
nachfolgenden Rechnungen auch Grenzschichten ry 
bei Strémungsvorgangen in flachem Wasser be- [ | a T= ~ 2 
handeln, was im Hinblick auf die Anwendbarkeit EEE 4 ‘ 
che a ee | x =(a,-l4)t 
der Flachwasseranalogie in der Gasdynamik einige tg 
i L=Ayt =i 
Bedeutung hat}+, : 79-0 
2. Grundgleichungen und Vernachlassigungen. x=0 
; Ae 2 7=7 
Zur Zeit t = 0 sei die ebene Verdiinnungswelle vee ; F LS OE er 
od. xpansionswelle an einer ebenen and; oor atensystem 
an der Stelle x =0 der ebenen Wand konzen- 3 und Bezeichnungen. : 


triert, fiir t > 0 laufe sie nach rechts in das Ge- 

biet x >0 hinein, der Wellenkopf mit der Geschwindigkeit a), das Wellenende mit a,—U, (Abb. 1). 
Dieser Fall ist gerade im Rohrwindkanal verwirklicht (wenn man hier einmal davon absieht, 
daB die Rohrwand dort im Gegensatz zu unserer ebenen Wand gekriimmt ist). Bekanntlich 
sind Geschwindigkeit U und Dichte 9 in der Welle an der Stelle x und zur Zeit t Funktionen einer 
einzigen Variablen 7, wobei 


ay Cee cap | x (1) 
ay t dy t 


ist; 7 bedeutet also den mit apt dimensionslos gemachten Abstand vom Wellenkopf. Fiir ein Gas 
mit ~ = 1,4 gilt: 
U § 
i 2 
a 2 6 1) ( ) 


(das negative Vorzeichen steht, weil die Geschwindigkeit nach links gerichtet ist) und 


s=(-a ® 


1 C. Donaldson u. R. D. Sullivan, The effect of wall friction on the strength of shock waves in tubes and 
hydraulic jumps in channels, NACA TN 1942 (1949). cian 

2 P. W. Huber u. D. R. McFarland, Boundary layer growth and shock attenuation in a shock tube with 
roughness, NACA TN 3627 (1956). 

3 H. Ludwieg, Z. Flugwiss. 3 (1955), S. 206. ae aD 

4 H. Bémelburg, Die praktische Anwendung der Wasseranalogie in quantitativer Form auf spezielle 
Probleme der Gasdynamik, Mitt. Max-Planck-Inst. Strémungsf. Nr. 10 (1954). 
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Der Berechnung der Grenzschicht wird der Impulssatz fitr instationare, kompressible Grenz- 


schichten zugrunde gelegt, wie er sich nach Routinemethoden aus den Grenzschichtgleichungen — 


durch Integration herleiten laBt: 


. 0 ee Ol a -—e 
Us | (o— 0.) de FOU 0) Ot ie Oe ae (4) 


0 
Dabei ist z die Koordinate senkrecht zur Wand, 0. und U sind die Werte von Dichte und Ge- 


; : x ‘ 
schwindigkeit in weiter Entfernung von der Wand, 1,, ist die Wandschubspannung, 0* und U | 


sind folgendermafen definiert: 


— | t Q it dz (Verdrangungsdicke) , (5) 
Qs 
0 
ies [ Bs ( dz (Impulsverlustdicke) , (6) | 


0 
wobei 9 die Dichte, u die Geschwindigkeit im Abstand z von der Wand bedeuten. 


Fiir die weitere Rechnung werden folgende Annahmen gemacht: 
I. Die Machzahl M der stationaren Strémung hinter der Welle sei nicht zu groB und bleibe 


jedenfalls unter dem Wert 1 (beim Rohrwindkanal ist dies der Fall, da hier M im Rohr den | 


Wert 1 nicht iiberschreiten kann. Man darf dann wohl annehmen, daB iiber die ganze Grenz- 


schichtdicke mit geniigender Naherung @ = Q, ist. Bei stationarer Strémung ist dies be- | 


kanntlich der Fall, da dort der durch 0 4 0. gekennzeichnete Einflu8 der Kompressibilitat fiir 
kleine Machzahlen vernachlassigbar ist. Unsere in- 


Sue) 80 Po, %o stationare Grenzschicht kann nun aber hinsichtlich 
= Druckanderungen tiber die Grenzschichtdicke als 
a Laos quasistationar behandelt werden, wenn nur die Grenz- 


schichtdicke klein ist gegeniitber der Ausdehnung 


in der sich Druckanderungen ther die Grenzschicht- 
dicke ausgleichen kénnen, ist nadmlich dann viel 
kleiner als die charakteristische Zeit fiir Druck- 
anderungen auBerhalb der Grenzschicht; jene ist 


Abb. 2. Konzentrierte Welle an einer ebenen Wand; Schallgeschwindigkeit, diese durch das Verhaltnis 
bed pane ica cet cigar as der Ausdehnung der Welle in x-Richtung zur Schall- 
geschwindigkeit gegeben!. — Damit entfallt das erste Glied in (4). Fiir das entsprechende Problem 
bei Flachwasserstrémungen ist die Vernachlassigung tibrigens exakt richtig, da sich die Dichte 
des Wassers ja nicht dndert. 


II. Das Verhaltnis von Impulsverlustdicke und Verdrangungsdicke sei konstant und die | 


Wandschubspannung hange nur von der Grenzschichtdicke ab (also von 6* bzw. #). Diese an 
sich recht einschneidende Annahme dirfte zu nicht allzu groBen Fehlern fiihren, da in der betrach- 
teten Strémung in Strémungsrichtung durchgehend Druckabfall herrscht?. Sicherlich geben trotz 
dieser Vernachlassigung die weiter unten hergeleiteten Naherungsformeln die hier wesentlichen 
Effekte der Instationaritat und des starken Abfalls der Dichte mit wachsender Entfernung vom 


Wellenkopf richtig wieder, da dieser EinfluB bedeutend gréfer ist als ein durch die Annahme IT 


zu erwartender Fehler. 


3. Laminare und turbulente Grenzschicht hinter einer konzentrierten Welle. Zur spateren An- 
wendung wird zunachst eine vereinfachte Strémung betrachtet: Die Verdiinnungswelle sei wahrend | 
hrer ganzen Bewegung auf eine scharfe Unstetigkeit konzentriert, deren Fortpflanzungsgeschwin- | 


' Die Voraussetzung @ = Qo bedeutet auch T = T,, iiber die ganze Grenzschichtdicke. In diesem Fall 
findet kein Warmeiibergang an die Wand statt. Bei gut warmeleitender Rohrwand ist diese Voraussetzung 
verletzt und fiir gré8ere Machzahlen sind dann Abweichungen von den folgenden Ergebnissen zu erwarten. 


Fur ideal warmeleitende Wand, d. h. konstante Wandtemperatur T)kann man bei Annahme eines Temperatur- | 
profils und entsprechender Beriicksichtigung des ersten Gliedes in (4) auch den Warmeiibergang beriicksichtigen. | 
Diese Rechnung wurde vom Verf. ebenfalls durchgefiihrt und wird demnachst i. a. Zusammenhang veroffentlicht. | 


2 Vel. hierzu Ziff. 7. 


der Welle in x-Richtung. Die charakteristische Zeit, — 


durch das Verhaltnis der Grenzschichtdicke zur | 
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digkeit mit V bezeichnet wird (Abb. 2). Man kann den ganzen Vorgang auf eine stationare Grenz- 
schichtstrémung zuriickfiihren mit nicht verschwindender Wandgeschwindigkeit, indem man die 
Geschwindigkeit —V iiberlagert. Dieser Fall wurde vom Verfasser 1951 nach dem Impulssatz 
behandelt!. Eine exakte Behandlung, wie sie von H. Mirels? gegeben wurde, fihrt fiir laminare 
Grenzschichten auf eine modifizierte Blasiusgleichung. Unsere Annahme II ist hier gut erfillt, 
da die Geschwindigkeitsprofile affin sind. 

Als Vorbereitung auf die Rechnungen in Ziff. 4 wird hier folgender, an sich recht umstandlicher 
Weg eingeschlagen: In konsequenter Auslegung der Annahme II (Ziff. 2) wahlt man ein Geschwin- 
digkeitsprofil, namlich 


u z z 
ap rueeeg, Hites Oe. 2 = 0.5 (7) 
dann wird 
1 
= 
O* = 5 Ors (8) 
6 
Tip 2 
o U2 QUo (10) 
(u = dynamische Zahigkeit). Der Impulssatz (4) ergibt dann mit U = — U, und den Abkiirzungen 
= 1 
P - (11) 
folgende Gleichung: 
1 dp Leeyens 12 
5 ee oe 
wobei »; = //;/0, ist. Die Anfangsbedingung ist 
lh 
= 0 fii = —t, 13 
yp tie U, (13) 


deb eo =O) fans Ve: 
Die Gleichung (12) kann nach der iiblichen Methode durch Zuriickfiihrung auf folgendes System 
gewohnlicher Differentialgleichungen gelést werden?: 


‘A charakteristische Grundkurven 
1 1 


Tae acl eee p=2v,. 14 Neate 
t=>, § i?) P= Ay (14) \ coh we \ 
Dabei bedeuten Punkte Differentiation nach einem Para-  \ e \aie \ \ 5-[fe)t 
meter, der sich leicht durch Division je zweier der Glei- h \ \ 7 
chungen (14) mit folgendem Ergebnis eliminieren 1aBt: ee 
cs 2 
oe ee 15 
7 ae (15) 
Al Abb. 3. s/t-Ebene fiir eine konzentrierte Welle 
OP = 12) V1. (16) mit charakteristischen Grundkurven. 
dt 


Gleichung (15) definiert eine Schar von Geraden in der (s, t)-Ebene (,,charakteristische Grund- 
kurven“), Gleichung (16) gibt das Wachstum von ¢ langs dieser Geraden (Abb. 3). Die Lésung 
y(s, t) ist in diesem ganz einfachen Falle offenbar eine durch s = (V/U,) t gehende Ebene. Man 
kann die Gleichung dieser Ebene sofort der geometrischen Anschauung entnehmen. Sie ist 


2 So 12, ( Z] ; ifm (17) 
i 4 
Durch Einfiihrung der Zeit t, wobei 
t=t— (18) 


1 E. Becker, Merkwiirdige Erscheinungen beim Wasserschwall, Dipl.-Arbeit, Gottingen 1951; siehe auch 


. Becker, Ing.-Arch. 21 (1953), S. 42. geo ; 
; 2H. Mirels, fener ey ee layer behind shock advancing into stationary fluid, NACA TN 3401 (1955) 


3 Courant-Hilbert, Methoden der mathematischen Physik. Bd. II, Kap. 1 u. 2, Berlin 1937, 
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ist, also derjenigen Zeit, die an der Stelle x vergangen ist, seit die Welle diesen Ort uberstrichen 
hat, laBt sich aus (17) folgender Ausdruck fiir 6 entnehmen: 


6=/Rmt Ay] (19) 
an f 3) =? 1 (20) 


peat ipa 
l4+a7 


Fir U,/V = 0 erhalt man 6 = V 129, ca-bzw) 0" — /4/3 J», . Dies entspricht der Lésung fiir 
die ruckartig mit der Geschwindigkeit U, in Bewegung gesetzte ebene Wand im inkompressiblen 


Medium, die sich auch exakt berechnen laBt (nach Rayleigh”). Der Abminderungsfaktor f(U,/V), _ 


der in Abb. 4 dargestellt ist, tragt der endlichen Wellengeschwindigkeit V Rechnung. 
In dem anderen Grenzfall U,/V—> co erhalt man mit Vt = x:6 =/302/U, . Dies ent- 
spricht der Grenzschicht an einer mit U, angestromten ebenen Platte in der Entfernung x von der 


Vorderkante. 
Fir die turbulente Grenzschicht werden folgende Naherungsformeln benutzt!: 


ff we ee 
if zs) fir 0<2<6, (21) 
1 
joe (22) 
8 
p= 16 23 
i (23) 
ie. U #\-14 
aia = 00125. eA . (24) 
a Durch eine ganz analoge Rechnung 
wie oben erhalt man hier 
Qe y, \15 U 
fee 0,308 (74) Nch( 7): (25) 
4 
Y Der Wert 
0 HF 08 72 76 20 eee 
Abb. 4. Funktionen f; und f, Cl. (20) und Gl. (26). 6 = 0,303 (ay ) U,t 


entspricht der Grenzschicht an einer ruckartig in Bewegung gesetzten ebenen Wand. Der Ab- 
minderungsfaktor f, ist ebenfalls in Abb. 4 dargestellt, er ist durch 


ee : 


ie +55) 


(26) 
DVS 
gegeben. Fiir U,/V— co erhalt man hier genauso wie bei der laminaren Strémung die Grenz- 


schicht an einer stationdar mit U, angestrémten ebenen Platte, wenn x = Vr die Entfernung von 
der Vorderkante bedeutet. 


4, Laminare Grenzschicht in der Expansionswelle. Der Impulssatz (4) wird nun (unter Vernach- 
lassigung des ersten Gliedes auf der linken Seite) auf die laminare Grenzschicht in der Expansions- 
welle angewandt. Unter Benutzung von (8), (9) und (10) und nach Einfiihrung der abhangigen 


Variablen 
U 2 
g= ( “ (27) 


A Qo 


ergibt eine langere, aber elementare Rechnung aus dem Impulssatz (4) 


min tobe aes eee (28) | 


6a) Ot 18 agt 15 dy ox Cy Oy By \\ Gp 


1 H, Schlichting, Grenzschichttheorie, 1. Aufl., S. 63. Karlsruhe 1951. 
2 H. Schlichting, a. a. O., S. 393. 


XXV. Band 1957 Becker: Das Anwachsen der Grenzschicht in und hinter einer Expansionswelle 159 


Nun ist es vorteilhaft, neue unabhangige Variablen in (28) einzufiihren, namlich die durch (1) 
definierte GréBe 


Gate 


Die Koordinate x werde beibehalten, der konsequenten Schreibweise wegen aber in € umbenannt: 
§ =x. Fir die Ableitungen gilt dann 


a ee en 10 (1—vn/ a 
Ox —sE ey on dy) at a On * a2) 
Aus (28) ergibt sich dann mit (2) und (3) 
an ik We Op tea le 7\6 
see ee omar n& A Qo Mo ( _ :) ches e 


Diese Gleichung ist mit der Anfangsbedingung pm = 0 fiir 7 = 0 zu lésen. Das dem System (14) 
in Ziff. 3 entsprechende Gleichungssystem ist fiir die partielle Differentialgleichung (30) 


6 1 — 
Veaee O2) (31) 
Ses (32) 
5 l1—y7 fy JB 7 \5 
2 = -(1——) 2577. 
a aera are agra ( 4 n (33) 
Aus (31) und (32) folgt als Gleichung fiir die Schar der charakteristischen Grundkurven 
dy J—y7 
a= t=" 2), (34) 
mit der Lésung 
c(l1—n) 


= Baayen” ae 
C ist der Scharparameter; eine einzelne Kurve der Schar ist durch einen bestimmten Wert von C 
festgelegt. Langs einer solchen charakteristischen Grundkurve gilt fiir g die aus (32) und (33) 
folgende, gewéhnliche Differentialgleichung: 

dy 7 Mo m\8 25 7? g 

= States : 36 

dq * 339? Ho 41 Qo ( q (3 —2) 1—y (36) 
Auf einer charakteristischen Grundkurve ist € aber durch (35) gegeben und nach Einsetzen in 
(36) kann man durch einfache Quadratur sofort y bestimmen. Das Ergebnis dieser Rechnung, 
das die Randbedingung » = 0 fiir 7 = 0 erfiillt, ist 


' Cane ae 
Ba LU] Méo eel (n| lee 
— "0 95 € (3 —2 n)/ dn. (37) 


Ay Oo . (3 — 2 )8 


Unter dem Integral steht der Ausdruck j4/j/j. Hierfiir darf man in unserer Naherung sicher die 
Sutherlandsche Formel benutzen!, die die Zahigkeit eines Gases in Abhangigkeit von der Temperatur 


gibt: 
1/2 
ys ia gE yee ; (38) 
Mo Ty 0,505 
T/T, 


In der Verdiinnungswelle gilt nun 


el aaa (39) 


Setzt man (39) in (38) ein und beschrankt sich auf eine lineare Naherung, so wird 


el 0,29 174 (40) 


Mo 


1 4, H. Shapiro, The dynamics and thermodynamics of compressible fluid flow, Bd. II, S. 1059, New York 
1954. 
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In (37) wird weiterhin C nach (35) durch € und 7 ersetzt. Wenn man dann beachtet, daB nach 
der Definition 


a 41) 
] —7 ae Ay t ( ) 
ist und da @ durch (27) definiert ist, so ergibt sich nach einiger Umrechnung 
& = v9 t 8\(1) (42) 
mit 
" 
(1 — 0,28 7) (1 a 
eal: —F 
(3 — 2 ) 
(4) = 36 7 10 (8—27)* n° dy, « (43) 
ie 
6 | 
0 


Die Funktion g(7) ist in Abb. 5 im Intervall 0 < 7 <1 aufgetragen’ (j = 1 entspricht gerade 
der Machzahl M = 1; denn zwischen 7 und M besteht die Beziehung M = 5 n|(6—n)). 
In der Umgebung von 7 = 0 laBt 
VA sich g, linear annahern: 


60 gin) ~ 47. (44) 
Daher wird fiir sehr kleine 7 
6 ~Ayytn. (45) 

Fiihrt man wieder nach (18) die Zeit t 
ein (mit V =a,), die seit Voriiber- 
streichen des Wellenkopfes an einer 
bestimmten Stelle x verstrichen ist, so 
erhalt man aus (45) 
“a O Art. (46) 
Zz Dies entspricht der Rayleighschen Lisung 
Be an einer mit zeitlich linear anwachsender 
Geschwindigkeit, also konstanter Be- 

Py, schleunigung in Bewegung gesetzten 
41 ebenen Wand. Die Lésung laBt sich 
ZA. auf einfachem Wege direkt aus der Be- 
"9 r wegungsgleichung unter Vernachlassi- 
gung der konvektiven Glieder herleiten 
bzw. direkt aus dem in diesem Falle 
LN Qo verkiirzten Impulssatz 


0 
+ (¢9 U 5%) =, (47) 


Mit (8) und (10) und mit U/U =1/t 
2 (bei linearer Zeitabhangigkeit von U) 
erhalt man 


00? 07 
= +2—=R», (48) 


IH (7) 


0 G2 04 LAT) 10 


Abb. 5. Funktion g)(7j) mit der Naherung (44). 


und hieraus ergibt sich sofort (46), wenn man ¢ durch T ersetzt. 

Die Naherung (44) fiir g,(7) ist zusammen mit g(7) in Abb.5 aufgetragen. Die Funktion 
gi(7) weicht von der Naherung nach oben, d.h. im Sinne erhéhter Grenzschichtdicke, ab. Der 
EinfluB der endlichen Wellengeschwindigkeit allein lieBe nach Ziff. 3 bzw. Abb. 4 eine Abweichung 
nach unten, d. h. in Richtung kleinerer Grenzschichtdicken, vermuten, da f; < list. Diesem Effekt 
wirkt aber der Abfall von Dichte und kinematischer Zahigkeit mit wachsender Entfernung vom 
Wellenkopf (wachsendem 7) entgegen und treibt die Grenzschichtdicke iiber die Rayleighsche 


Naherung hinaus. 


5. ZusammenschluB der Grenzschicht in der Expansionswelle mit der Grenzschicht hinter der 
Welle. Das Wellenende sei durch den Wert 7, festgelegt. Die Grenzschichtdicke 5, am Wellenende 
ist dann durch 6 = yt g(7,;) gegeben. Hinter der Expansionswelle ist der Verlauf der Grenz- 


} Fur die Berechnung der Funktion gi(7) sowie auch von gt(7) danke ich Herrn Dr. E. Martensen. 
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schichtdicke durch (17) bestimmt, wenn man dort V geeignet wahlt, namlich derart, daB am Ort 
des Wellenendes die aus (17) folgende Grenzschichtdicke mit 5, iibereinstimmt. Dies wird auf 
foigende Weise erreicht: Gleichung (17) ergibt 


= Doye(I at alle: a (49) 
Am Wellenende ist x/(a)t) = 1— », und a,/U, = 6/(5 1,), also 


(=m) 5 aca (50) 


Oo 12 OE 


Wenn dies mit 0 iibereinstimmen soll, mu8 folgende Gleichung gelten: 


@ GF 


12» |1— (1m) 557] FE (= rosie) (51) 


Diese Gleichung kann fiir vorgegebenes 77, nach U,/V aufgelést werden. Unter Benutzung des 
Ausdrucks (20) fiir f, erhalt man 


Vo 
U,; Ne 12». a "eumh) (52) 
= = es ‘ 5) 
V 6 (Ce — mn) 
as Vy gum) + PBS o 


Damit liegt auch die Grenzschicht hinter der Expansionswelle fest. 


6. Turbulente Grenzschicht in der Expansionswelle. Analog zu der nunmehr etwas ausfiihrlicher 
dargelegten Abschatzung fiir die laminare Grenzschicht lat sich auch die turbulente Grenzschicht 
naherungsweise erfassen. Die Rechnungen verlaufen ganz den obigen Rechnungen im laminaren 
Fall entsprechend, so daB sie hier nur in grofen Schritten wiedergegeben werden sollen: 

Man geht mit den Beziehungen (20) bis (23) in den Impulssatz (4) ein (in dem wieder das erste 
Glied auf der linken Seite vernachlassigt wird) und fiihrt als abhangige Variable 


pai sy" (53) 


29 A 


ein. Die unabhangigen Variablen werden wie bei der laminaren Grenzschicht als € = x und 
y = 1—vx/(a)t) gewahlt. — An sich gelten die Beziehungen (20) bis (23) nur fiir turbulente 
Grenzschichten ohne Druckanderung in guter Naherung. Aber entsprechend dem in Ziff. 2 zur 
Begriindung der Annahme II Gesagten darf man annehmen, daf diese Naherung auch in unserem 
Falle noch ausreicht, um die hier wesentlichen Effekte der Instationaritat und der Dichteanderung 
mit wachsender Entfernung vom Wellenkopf zu erfassen. Eine langere, elementare Umformung 
von (4) ergibt dann fir y die der Gleichung (30) entsprechende Beziehung 


30 Op L007 155 (1 We He Eee fos ao 54, 
54 OF * 54 EH Ua 6 (i Cy 4% cy) 


pee Nee (55) 


gegeben, und langs dieser Kurven gilt fiir p 


; 5 2 i) UES fap, SE 
eee ee as Gt es EE SRS ae (56) 
dy o4: uy) “ \ Ho a : 


5417 


Die Lésung von (56) wird unter Beriicksichtigung von (40), (53) und (55) und der Anfangsbedin- 
gung gm = 0 fir 7 = 0 


y NT} 
O°" = agt (2) “Bi(") oy) 
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ry 154 N- 


19 \19 7 \° 
as he 1 == 0,087) (tae 


g(n) = 0,197 25 173 nf dy « (58) 


n\4 si Ng 


Die Funktion g,(7) ist in Abb. 6 dargestellt. Fiir 7 ~ 0 wird g(7) angenahert durch 


Bi(7) ~ 0,066 7" . (59) 


Man erhalt dann durch Einfihrung von 7 (siehe Ziff. 3 bzw. 4) fiir 6°/4 den Ausdruck 


8°!4 ~ 0,075 [U(t) P4 v2 t . (60) 


Dies entspricht der Lésung fiir die turbulente Grenzschicht an einer mit konstanter Beschleunigung 


in Bewegung gesetzten ebenen Wand. Die Lésung (60) 1aBt sich analog wie bei der laminaren | 
Grenzschicht sofort aus dem verkiirz- 


I ten Impulssatz herleiten. 
Die Naherung (59) ist zusammen 
mit g,(7) in Abb. 6 aufgetragen. Im 


G06 


allein eine Verminderung der Grenz- 
schichtdicke gegeniiber den durch die 
Naherung gegebenen Werten zu er- 


G04 


ker sein mu als bei laminarer Stré- 
mung. Die Verdickung durch den Ab- 
fall von 9 und y mit wachsendem 7 
reicht hier nicht aus, um g,(7) iiber 
die Naherungskurvezuheben, sondern 
beide Effekte kompensieren sich der- 
art, daB im ganzen Bereich die Nahe- 
rung (59) recht gut gilt. Diese Uber- 
einstimmung ist aber, auBer in der 
Umgebung von 7 =0, rein zufalliger 
Art. 

Der AnschluB der durch obige 


G03 


G02 + 


007 


07 
Formeln gegebenen _ turbulenten 
0 Ge Gt G6 98 40 Grenzschicht an der Stelle 7=7, an 
Abb. 6. Funktion g,(7) mit der Niherung (59). die Grenzschicht hinter der Welle 


kann auf entsprechende Weise wie 


bei der laminaren Grenzschicht hergestellt werden; die Grenzschicht hinter der Welle ist dabei | 


durch (25) gegeben. Die Durchfiihrung der Rechnung wird hier iibergangen und nur das Ergebnis 
fiir U,/V mitgeteilt: 


_ (U,\3/4 VY) \r/4 

lag) a) te 3 

V ie % \1/4 Were U; 3/4 6 d= 7%) ( ) 
9 \>, gi(m) + 9, a 5 


Hiermit ist die turbulente Grenzschicht auch hinter der Expansionswelle bekannt. 


7. Ausblick. Man kann nun noch iiberschlaglich nachpriifen, wenigstens fiir laminare Strémung, 
ob die in Ziff. 2 gemachte Voraussetzung, das Verhaltnis 6*/) sei konstant und die Wandschub- 
spannung hinge nur von der Grenzschichtdicke ab, einigermafen giiltig ist. Bei einer genaueren 
Berechnung muf} man diese Annahme sicher fallen lassen und mindestens einen Formparameter 


Gegensatz zur laminaren Strémung | 
(Abb.5) liegt die exakte Kurve hier | 
unter der Naherung. Dies ist einzu- || 
sehen, wenn man bedenkt, daB durch | 
die endliche Wellengeschwindigkeit — 


warten ist, die hier wegen f,<f, star- _ 


XXV. Band 1957 Becker: Das Anwachsen der Grenzschicht in und hinter einer Expansionswelle 163 


einfiihren. Dieser Formparameter ist beim Pohlhausenverfahren! bekanntlich definiert als 


_ 1 ee &/1aU , aU 
iS 0 U Ox » ae eal (62) 


Setzt man hier fiir 6 die Lésung (42) ein, so erhalt man fiir 2 den Ausdruck 


_ (L—1/6)* ge(n) 
A 1—0,28y 7 * : i) 


Aus (63) folgt, daB 2 im Bereich 0 << 1 in den Grenzen 4,0 > A > 2,7 bleibt. Fiir solche Werte 
von A weichen aber 6*/% sowie die Wandschubspannung noch nicht wesentlich von den Werten fiir 
A = 0 (keine Druckanderung in Strémungsrichtung) ab, so daB es berechtigt erscheint, den EinfluB 
von / in erster Naherung iiberhaupt zu vernachlassigen; dies entspricht dann der Annahme in 
Ziff. 2. 

Fiir héhere Genauigkeitsanspriiche kann man die Naherungsformel (42) durch Iteration ver- 
bessern, indem man von vornherein mit dem Pohlhausenpolynom in den Impulssatz eingeht, in 
erster Naherung aber / = 0 setzt. Dann erhalt man das Ergebnis (42) (mit etwas anderen Werten 
fiir g,(7), da unser Ansatz (7) nicht mit dem Pohlhausenpolynom identisch ist). In nachster Naherung 
setzt man den nach (63) berechneten Parameter A als bekannte GréBe in den Impulssatz ein und 
bestimmt die Grenzschichtdicke von neuem. Eine Abschatzung hat ergeben, da sich bei dieser 
Iteration die Grenzschichtdicke um nicht mehr als héchstens 10% gegeniiber dem aus (42) folgen- 
den Wert andert. Auf die Durchfiihrung der Iteration wurde daher der Einfachheit halber und in 
Anbetracht der ohnehin begrenzten Genauigkeit des Verfahrens verzichtet (z. B. bringt die An- 
nahme 0 = ¢@,, fiir gréBere Machzahlen Fehler von derselben GréBenordnung). 

Zum SchluB sei auf folgendes hingewiesen: Auf grundsatzlich ahnliche Art wie oben die von 
Anfang an laminare bzw. turbulente Grenzschicht berechnet wurde, kann man auch den Fall be- 
handeln, da die Grenzschicht anfangs laminar ist und erst in einiger Entfernung vom Wellenkopf 
in die turbulente Grenzschicht umschlagt. Die Anfangsbedingungen fiir die turbulente Grenzschicht 
sind hier etwas komplizierter und der ZusammenschluB von laminarer und turbulenter Grenzschicht 
ist nicht so einfach zu vollziehen wie etwa oben der ZusammenschluB gleichartiger Grenzschichten 
an der Grenze der beiden verschiedenen Wellengebiete. Durch Zurickgreifen auf die hier benutzte 
Charakteristikentheorie lieBe sich aber auch dieses Problem in jedem einzelnen Falle lésen, worauf 
jedoch in dieser, zur ersten Orientierung iiber die zu erwartenden Effekte gedachten Arbeit ver- 
zichtet wird. 


(Eingegangen am 24. Juli 1956) 
Anschrift des Verfassers: Dr. E. Becker, Gottingen, Bunsenstr. 10. 


1 Vel. H. Schlichting, a. a. O., Kap. XII, b. 
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Biegetheorie der Rotationsschale mit flacher, kreisférmiger Erzeugender 
Von R. Windels 


1. Einleitung. In den letzten Jahren ist eine groe Anzahl von Veréffentlichungen tiber Schalen 
erschienen, und es ist heute méglich, fiir jede Schalenform die zur Ermittlung des Spannungs- und 


Verschiebungszustandes erforderlichen Ausgangsgleichungen ohne groBen Rechenaufwand zu ge- | 
winnen. Es sei hier auf eine Ubersicht von Zerna! verwiesen, der Hinweise auf das vorhandene | 


Schrifttum gibt. 


Die hauptsichlichen Schwierigkeiten der Schalentheorie treten bei der Lésung der Gleichungs- | 


systeme auf. Schon bei der geometrisch einfachsten Schalenform, der Kreiszylinderschale, sind sie 
nicht unerheblich, und die Ermittlung des Spannungsverlaufes erfordert langwierige Berechnungen.. 
Bei doppelt gekriimmten Schalen ist es in den meisten Fallen schon fiir den Membranspannungs- 


zustand notwendig, ein Naherungsverfahren, wie zum Beispiel das Differenzenverfahren, anzu- | 


wenden. 


Es scheint wiinschenswert, fiir weitere Schalenformen Lésungsmethoden anzugeben. In dieser | 
Arbeit wird der Versuch unternommen, das Spannungsbild fiir eine Rotationsschale mit flacher, | 


kreisférmiger Erzeugender nach der Biegetheorie zu ermitteln. 


ro 


Qh 


pee (et, frotations- 
ie z | 


tore Fs 
17 | Rotations- 
. & achse 
Why 


Z 
Abb. 1. Schnittebene mit Rotationsachse Abb. 2. Schnitt senkrecht zur Rotations- Abb. 3. Schale mit negativem 
(Meridianschnitt). achse (Breitenkreisschnitt). Kriimmungsmab. 


2. Schalenform und Bezeichnungen. Wenn die Einheitsvektoren i,, ij, i; ein kartesisches Rechts- | 


system aufspannen, kann der Ortsvektor der Mittelflache einer Rotationsschale 
t =r [t, sin 0 + @ (iz sin p — i, cos g)] (1) 


geschrieben werden. Darin sind (siehe Abb. 1 und 2) # und @ zwei Parameter, die ein Koordinaten- 
system auf der Schalenflache bilden. Fir = konst. erhalt man die Schar der Breitenkreise, fiir 


gy = konst. die Schar der Meridiankurven. Der Abstand eines Punktes der Mittelflache von der | 


Rotationsachse betragt 
To =1 Qo- (2) 
Fiir den hier betrachteten Sonderfall, bei dem die rotierende Kurve ein Kreis mit dem Halbmesser r 
ist, wird 
0 =o +— cos? . (3) 
Der Mittelpunkt des rotierenden Kreises hat den Abstand @ r von der Rotationsachse. Bei Schalen 
mit negativem Kriimmungsmaf bleibt r eine positive GréBe, die positive Richtung des Winkels 
ist dann aus Abb. 3 zu ersehen, und fiir gy gilt, wie auch in allen folgenden Gleichungen, das untere 
Vorzeichen. 
Die Schale wird durch zwei Binderscheiben ) = konst. und zwei Randglieder g = konst. be- 
grenzt. Die Untersuchung soll sich auf den Membranspannungszustand und seine Stérungen durch 
die von den Randgliedern ausgehenden Biegespannungswellen beschranken. Wie die Berechnung 


der Kreiszylinderschale zeigt *, sind die Stérungen des Membranspannungszustandes an den Bin- 
dern von untergeordneter Bedeutung. 


! W. Zerna, Beton- und Stahlbetonbau 48 (1953) S. 88. * U. Finsterwalder, Ing.-Arch. 4 (1933) S. 43. 
° I’. Dischinger, Beton und Eisen 34 (1935) S. 257. 


} 
} 
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Der Binderabstand wird mit | bezeichnet, und die Lage der Binder wird so gewahlt, daB der 
Breitenkreis ? = 0 von jedem Binder gleich weit entfernt ist. Als wesentliche Einschrankung wird 
vorausgesetzt, daf die Schale sich zwischen den Bindern nur wenig wolbt. Sie soll in der )-Richtung 
als flach angesehen werden kénnen, eine Forderung, die erfiillt wird, wenn die Binder nicht mehr 
als } ~ + 0,3 (~ 20°) vom Koordinatenursprung entfernt sind. 

Die Verschiebungen der Schalenmittelflache in den Tangentenrichtungen des Meridians und des 
Breitenkreises werden mit u und v bezeichnet, die Verschiebung in Richtung der Schalennormalen 
mit w (positiv nach innen, vgl. Abb. 4a). Der Abstand eines Schalenpunktes von der Mittelflache 
sei 2 (positiv nach auBen). 

Auf die Schale wirkt die Flachen- 
last p, deren Komponenten in Rich- ? Ba) 
tung der Meridiantangente p,, in Rich- 
tung der Breitenkreistangente p,, und 
in Richtung der Schalennormalen p, a N 
(positiv nach innen) sind (vgl.Abb.4b). aa M 

An einem Schalenelement greifen aa, ba oe 
die in Abb. 4c,d dargestellten Langs- eas a ie 
krafte N,, N,, die Schubkrafte NV,,, — ¢ Mp d My 

N,,, die Biegemomente M,, M,, die 
Drillungsmomente M,,,, M, , und die 
Querkrafte Q,, Q, an. 
Ableitungen ack @ werden durch einen Strich, Ableitungen nach g durch einen Punkt an- 
-gedeutet. 

Die folgenden Betrachtungen sollen fiir eine anisotrope Schale gelten, fiir die Dischinger! die 
Dehnungssteifigkeiten D,, D,,, D, sowie die Biegesteifigkeiten B,, B,,, B, einfiihrt. Es werden 

die Abkiirzungen 


Abb. 4. Schalenelement. 


xp? xp? 


Ds Dx gy 
ise 7D,” ( 

eee ieee ee 

x By ? xp By 


benutzt. Fiir den Sonderfall der isotropen Schale mit der Schalenstirke t folgt 


Et } 
D D, Die D; = — ? 
Ee 
, ges eT al (5) 
6, = 6,9 = | 


eet J 


3. Gleichgewichtsbedingungen. Die Aufstellung der Gleichgewichtsbedingungen fiihrt man 
aweckmaBigerweise auf eine in der Literatur aN: Art durch, z. B. nach Fliigge®. Dabei 
empfiehlt es sich, die folgenden reduzierten Schnittkrafte einzufiihren: 


== = 1 

N, = N,» N= Ng, 

; ; 2 

= — 1 

M,, = 0M,» M, sre ea (6) 
0 


Q, ae C0 Oz, Oy a, ms 2 


Dann erhalten die Gleichgewichtsbedingungen folgende exakte Form: 
Nz + Nz, + o sin dN, £O,/r +r ooP,=9, 
, nT: 1 . i} — 
Neg + No F7-sin 8 (Ney + Nos) = 


=00 rPp =O, 


r 09 


SS — 
= 
Jj 
ob) 
— 


PCOR Nee OF +Q,/r +17 0 p,=9, 


1 Siehe FuBnote 3 von S. 164. 
2 W. Fligge, Statik und Dynamik der Schalen, Berlin 1934. 
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| 
M, + M;,+ o, sin? M,—O, =0, 1 
/ er ES 1 . oa 
Meo, onc ae | (7b) 
ree ~~ 008 O My, —1 (Nzp— Ny) =). | | 
J 


4, Elastizitatsgesetz. Das Elastizitatsgesetz kann ebenfalls auf einem in der Literatur gezeigten 
Weg gefunden werden. Zerna! gibt alle erforderlichen Gleichungen in der nach der klassischen 
Elastizitatstheorie exakten Art an. Es kénnen jedoch, wie Green und Zerna? zeigen, weitere Ver- 
einfachungen eingefiihrt werden, ohne daB die Genauigkeit 
der Berechnung dadurch geringer wird. Diese Vereinfachun- 
gen bestehen im wesentlichen darin, daf} in konsequenter 
Weise diejenigen Glieder, die die Schalenstarke t als Fak- 
tor enthalten, gegeniiber allen itibrigen Gliedern in den. 
Gleichungen vernachlassigt werden. In der erwahnten Ar- | 
beit wird die Vektor- und Tensorrechnung angewendet, wo- | | 
durch eine groBe Allgemeingiiltigkeit erreicht wird ; anderer- | 
seits mégen dem mit dieser Schreibweise nicht vertrauten 
Leser Schwierigkeiten entstehen. Daher sollen im folgen-| 
den die erforderlichen Gleichungen in anschaulicher Weise) 
hergeleitet werden, was bei dem hier betrachtetem spe-. 
ziellen Fall keine Umstiande bereitet. Zunachst werden) 
noch einmal die in der Schalentheorie tiblichen Annahmen 
genannt: 


1) Die Schalenstarke ist klein gegeniiber allen anderen 
\ Abmessungen. | 
Abb. s< Moridiausebuitt: 2) Die Verschiebungen sind klein gegeniil er den Schalen- | 
abmessungen (auBer der Schalenstarke), Gré®en héherer 

Ordnung kénnen neben solchen erster Ordnung vernachlassigt werden. | 
3) Die Spannungen und Dehnungen in Richtung der Schalennormalen sind klein und werden 
vernachlassigt. 


4) Die Schalennormalen bleiben auch Normalen der verformten Mittelflache. 
Es werden zwei Schalenschnitte betrachtet, ein Meridianschnitt mit dem Halbmesser r (Abb. 5) | 


und ein Schnitt in Richtung der zweiten Hauptkriimmung mit dem Halbmesser rz = r,/cos #. Lm 
Abstand z von einem Punkt A auf der Schalenmittelflache liegt der Punkt P. Nach der Schalen- 


verformung wird die Lage der beiden Punkte mit A und P bezeichnet. Die Verschiebungen des 
Punktes P seien u,, v,, w,. | 


Aus dem Meridianschnitt (Abb. 5) folgt 


Teese Rat) a) Par Bs aula 1 (8) 


r Pip = B= Ob r l—w/(r + 2) 


Werden entsprechend der ersten und zweiten Annahme die Groen z/r und w/(r + 2) gegen 1 
vernachlassigt, so erhalt man daraus 


w’ 

u,=u+z—, (9) 

und auf demselben Wege aus dem Schnitt in Richtung der zweiten Hauptkrimmung | 
w 

auBerdem mit der dritten Annahme | 
w, = w. (11) 


} W. Zerna, Ing.-Arch. 17 (1949) S, 149. 


nes ee eenu. W. Zerna, Quart. Journ. Mech. and Applied Math. 3 (1950) S. 9 und Theoretical Elasticity, | 
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Die Dehnungen fiir einen Punkt P in Abstand z von der Mittelflache werden 


eon =~ (iw cos #) =e, +2 ne : 
; oo Lee (12) 
eee eal 
_worin die Dehnungen der Mittelflache 
| e, =—(w' Fw), ep = (0 —weos 8), Yep = w poo (13) 


sind. 

Das Hookesche Gesetz laBt sich fiir den zweidimensionalen Spannungszustand ¢,, 0, T.,, =T 

_=T wie folgt schreiben: cea 
E E 


Oi, as BE gy th oes)’ 0p = i eos eI ReN ie v 
Ili afi y 


px 


HE yp 
72 9 VxQ,2° 


(14) 


_ Durch Integration iiber die Schalenstarke erhalt man die Schnittkrafte. Wird dabei fiir die Breite 
eines Schalenelementes entsprechend der ersten Annahme 


| Pa 2 


=14361 as 


r 


pBesetzt (ebenso fiir 79), so ergibt sich mit Hilfe der Cleichungen (4), (5), (12), (14) das Elastizitats- 
| gesetz in der Form 


) N,= J o,dz= D, (e, + @ &) ; 
Ne jio,d:—D., era p ,)’, 


Nop = Npz = [t de =Dip— "yuo, 
M,=—fo,rdz=—B,(% +H); (16) 
M, =— Jo, 2ds=—B,( +m), 
; My = My, =— feds =— Boyz" [E+ Bele 
| Die beiden ersten Gleichungen lauten nach ¢,, ¢, aufgelést: 
= = alp— ae 2 pel pe Hp): (17) 


- Damit sind alle erforderlichen Ausgangsgleichungen aufgestellt. Man erkennt, daf die eingefiihrten 
Naherungen eine paarweise Gleichheit der Schubkrafte N,, = N,, und der Drillungsmomente 
M,,, = M,, bewirken. 


5. Membranspannungszustand. Beim Membranspannungszustand werden in den Gleichgewichts- 
_bedingungen (7) die Momente und Querkrafte vernachlassigt, und man erhalt 


Nz + Nj o,sin dN, +rop,=9, 
== 1 
Nee NG, tne (Ng Nos) LP (18) 


0, cos ON, + N, +7 Oop, =9. 
Aus der letzten Gleichgewichtsbedingung (7) folgt die Gleichheit der Schubkrafte N,, = Ny,. 
Die Lisung dieses Gleichungsystem wird von Riidiger’ und Zerna® gezeigt, darum soll hier nicht 
naher darauf eingegangen werden. 


1 D. Riidiger, Ing.-Arch. 22 (1954), S. 336. 
2 W. Zerna, Ing.-Arch. 17 (1949), S. 223. 
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Fiir eine Belastung infolge Eigengewicht 
=+peosgsin?, Po= P sing, P. = p cos 7 cos # (19) 


ergeben sich die Schnittkrafte zu 


— } 
set Pea ee 
=a 1 
Np =—2rplyt aay Ol cos, (20) 
1 tg O ; | 
a oe, ‘at 
Ne ee A rece ae ip ® sing. | 


In diesen Gleichungen bedeuten 
& = 0 (Hsin F + cos 8) -+ sin? 9 + C, sind +C,, 


D' =cos?(+sn0+o00+C). | 


Von der Richtigkeit der Gleichungen (20) kann man sich durch Einsetzen in (18) tiberzeugen. Die 
beiden Integrationskonstanten C,, C, lassen sich zur Erfiillung von zwei Randbedingungen be- | 
nutzen. Sollen zum Beispiel an zwei Binderscheiben r sin J) =-+ 1/2 die Langskrafte N, ver-| 


schwinden, erhalt man 


CG =o; C, = — 0 (YJ sin Vy + cos 9) sin? Dae (22), 


(21), 


Der Biegespannungszustand erfordert die Entwicklung der Krafte als Fourier-Reihen, daher | 
sollen auch die Membrankrafte infolge einer als Fourier-Reihe dargestellten Eigengewichtsbelastung 
ermittelt werden. Die spater erlauterten, fiir flache Schalen giltigen Naherungen 


i 


sind? ~ 0, cos? ~1, 0 ~~ =—et!l (23) | 


} 


ermoglichen die Berechnung der Schnittkrafte aus den Gleichungen (18) fiir die Belastung (EKigen- | 


gewicht) : 


‘pa—05 Po = Pn Sing cos 0, , Pz = Pn ©O8 P cos H, (24) | 
mit 
n I nm It 
B, = 59,9 =18, art ee he aie ae (25) 
worin #, den Wert von ? an den Binderscheiben darstellt (r sin ®g = -| 1/2). Die Richtigkeit der 


nachfolgend angegebenen Schnittkrafte kann wieder durch Einsetzen in die Gleichungen (18) | 
kontrolliert werden: 


eo Data yy 

Vi EL 1G, 008 P 098 On» 

Va 2 ¥ 1/0 | 

N,= —? Paap Pie: or ean (26), 
ae HPA, ; 

NS Baie, ein ee | 


Die Schalenverformungen infolge der Belastung (24) erhalt man nach einigen Zwischenrechnungen 
aus den Gleichungen (6), (13), (17), (26), wobei die Naherungen (23) beibehalten ia | 


py h — = a 0, 
= Daw PE Tepe 2 Oo: mye + 08 ( 1 + HQ) (0? F 1a) | 
+ 4(1 + u)/6,,,] cos y sin @,, , 
if Pn = oa a 0. 
Y= FD, =) FE Toop 2 eo #)/5. — @o(E 1 + u/e0) (A? + 1/0) 
+412 0 (1 + 1)/dz9] sing cos 8, . 
= Tan =» =. 2 i 
Y= DO BELA 2 les — HA )I6, + 0 (A? — los) (A? = 1/00) 
+- 41? (1 a {)/dx 01 CORY E08 a a 


7) 


— 
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Der Normalenverschiebung w sind entsprechend den Gleichungen (16) Momente zugeordnet, die 
jedoch im allgemeinen so klein bleiben, daf sie in den Gleichgewichtsbedingungen vernachlassigt 
werden kénnen. Dann bildet der Membranspannungszustand eine Partikularlésung des durch die 
Gleichungen (7), (16) gebildeten Systems von Differentialgleichungen. 


6. Differentialgleichungen fiir den Biegespannungszustand. Der Membranspannungszustand er- 
laubt es in den meisten Fallen nicht, an den Schalenrandern die Verformungsbedingungen zwischen 
Schale und Randglied bzw. Binder zu erfiillen. Diese Aufgabe muB von dem Biegespannungs- 


- zustand tibernommen werden. Es sollen hier nur diejenigen Randstérungen untersucht werden, 


die von den Randgliedern y = konst. ausgehen. Die an den Bindern = konst. entstehenden 
Stérungen sind von untergeordneter Bedeutung. 

Die Ausgangsgleichungen fiir den Biegespannungszustand sind die Gleichgewichtsbedingungen 
(7), in denen die Belastungskomponenten p,, p,, p, entfallen, und das Elastizitatsgesetz (16). 
Wie schon erwahnt, sollen hier nur solche Schalen betrachtet werden, die in der 0-Richtung flach 
sind, bei denen der Halbmesser Q9 7 sich also nur wenig verandert. Dann kénnen die Naherungen 
(23) eingefiihrt werden. Die Naherungen sin } ~ 0 und 0, ~ 0) mégen sehr grob erscheinen. Es 
ist aber zu beachten, daB fiir ® = 0 alle Gleichungen exakt bleiben, und das zu berechnende Span- 
nungsbild interessiert im wesentlichen gerade fiir den Bereich in der Nahe von ? = 0. Die mit 
Hilfe der Naherungen (23) ermittelten Schnittkrafte werden also in Schalenmitte (# = 0) nur 
geringfiigig von denjenigen abweichen, die die Gleichungen (7), (16) fiir alle Werte # exakt erfiillen. 

Green und Zerna? zeigten, da unter Beachtung der ersten Annahme (vgl. Abschnitt 4) in den 
beiden ersten Gleichgewichtsbedingungen (7) die Querkrafte vernachlassigt werden kénnen. Diese 
Vereinfachung bedeutet keine grébere Naherung als die zur Aufstellung des Elastizitatsgesetzes 
benutzten. 

Die letzte Gleichgewichtsbedingung (7) steht im Widerspruch zu der sich aus dem Elastizitats- 
gesetz ergebenden Gleichheit der Schubkrafte und Drillungsmomente. Sie wird daher nicht be- 
achtet. 

Die Gleichgewichtsbedingungen (7) und das Elastizitatsgesetz (16) erhalten bei Beriicksichtigung 
der genannten Naherungen und der Gleichungen (4), (6) folgende Form: 


Ne+Nye=0, NiegtNe=0, ON,+N.+-Gt+O=0, (28) 
= M+ M;,, 0. MM, 4M. (29) 
N—D,6,0,(6.+ 7%), N,=D,°28=, N.,=N,.=D, yO) 
i, eta) 
M, =— B,B,(w" + uw/e)>> | 
Te de ee. pee ! (31) 
Mey = My, = — By Boo aie © | 


Die Langs- und Schubkrafte kénnen durch eine Spannungsfunktion @ so ersetzt werden, da die 
ersten beiden Gleichungen (28) erfiillt sind: 


N= pO; Neo, Os Nig =— 0". (32) 
Werden nun diese Werte in die dritte Gleichung (28) unter Beachtung von (29), (31) eingefiihrt, 
ergibt sich die Differentialgleichung 


RT ee le ye Sie nate e ge 
Beene a 2 piogan pw (0, Ol aS a =) (33) 


worin 


B=B8,o +50 +6.—2 8,0) - (34) 
Aus (13) folgt die Identitat 
& + 03 £5 — Oo Veo + (Qo w+ w) cx 0, (35) 


1 Siche FuBnote 2 von S. 166. 
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und daraus mit (6), (16), (17), (32) eine zweite Differentialgleichung : 


6,08 @'" 12 002 @"" + B'" + (@ w" + 0) Ww’) (1 — p?) ar =, (36) 
worin 
Ox ia Ox 
es (2 = 6, (37) 
ist. 


Die beiden Gleichungen (33), (36) bilden ein System homogener partieller Differentialgleichungen 
vierter Ordnung. Die gesuchten Funktionen sind w und @, aus denen die Schnittkrafte nach (6), 
(29), (31), (32) ermittelt werden kénnen. Die Verschiebungen ergeben sich aus (13), (16), (17). 


7. Lésung der Differentialgleichungen. Durch den Ansatz 
w = W,e™? cos 8, , 
@ = F_e™? cos, 
folgt aus den Differentialgleichungen (33), (36) das homogene Gleichungssystem 


(,— 2 Bm? +m) W, + (LF oomt) e F, = 0, 


nk = Echt D 39) | 
(3,20 m? + mt) F, — (1 FQ m*) (1-44) 2 W, = 0 ae | 
mit der Abkirzung 
m = ae (40) 


fiir das nur dann eine nicht triviale Lésung existiert, wenn die Nennerdeterminante verschwindet, 
wenn also 


(8, —2B m® + mi) (0, —2.5m® + mt) + e8 (1 Gym)? = 0 (41) 
wird, worin 
= j/ r? D1 — p*) (42) 
054" Bo 
bedeutet. Gleichung (41) ist von achter Ordnung in m und besitzt als Lésungen die acht Wurzeln 
ih = — Mh, = — Hy ts 
Wa See es aa (43) 
Mm; = —M, = — Ha +12, 
mM, = — mg = — %z — 1 Ve 
mit den reellen Konstanten 
= x. — v 
moe te oge ae 


Die Wurzeln sind paarweise konjugiert komplex und paarweise entgegengesetzt gleich. Ihre Be- 
stimmung ist fiir den allgemeinen Fall der anisotropen Schale recht umstandlich, kann aber fir 
die spater behandelten Sonderfalle wesentlich vereinfacht werden. 


Aus den Gleichungen (39) leitet sich fiir die komplexen Konstanten W,, F. die Beziehung 


— nD. re Se ae — — D ™. mi — | 
We Qor 1 + @)m ee aes r A? 08 Owen 20 EHS (45) | 


r —P. be him Sa D, (1 — p?) 1 + 0 m? ? 


ab, die nach Multiplikation beider Gleichungen miteinander folgende Form erhalt: 


W =iKk)/=»—2+m F_rK a B 
) eaapee = Ky +iy)F, (46) 
mit 

ro 


~ YBp Ds — 1) 


(47) 


bent 8 (38) | 
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Die reellen Konstanten y, y lassen sich aus den Gleichungen (45) baw. (46) fiir jeden Wert m be- 
rechnen: ; 
1G we 
Aes 4 0 
x+y é Bx — 2 B m? 4 
1 by —2m® + ma! 


Laas an Le o;m?* 7 (48) 


an at) 
| rt 


dx — 2 dm? + m! 
Be =2 Pm 2 in* 


Darin bedeutet in Ubereinstimmung mit Gleichung (42) 


Ytip=t 


Te ea [Dz 2). 4. 
eed ae 


Mit den Wurzeln (43) folgt fiir y, y aus (45) 
M1, Ms > — YX, — Yy » 
M,™Mg ?— 4%, + Y» 
M3, Mz? + Xo, + o> | 
M4, Mg : + Xo, — Yo- J 


(59) 


Fir die praktische Berechnung wird von den komplexen Konstanten W,, F,, auf reelle Kon- 
stante W,, F iibergegangen. Bei allen folgenden Schritten soll dabei eine Beschrankung auf die 
ersten vier Wurzeln (43) vorgenommen werden, die fiir die Berechnung einer Randstérung aus- 
reichen. Zur Beseitigung der Stérung an einem zweiten Rand stehen die iibrigen Wurzeln zur Ver- 
fiigung, fiir die alle folgenden Rechenoperationen ebenfalls durchgefiihrt werden miissen. 

Nach Einfiihrung der neuen Konstanten 


W,=i(W,—W,), F,=i(F,—F,), (61) 

W,= (W; + W,), F, = (F; + Fy), . 

W,=i(W,— W,), F, =i(F; — F,) 
lassen sich die Gleichungen (38) in der Form 

w=([e“? (W, cosy, p + W, sin», ¢) 

J 

+e? (W, cos, gp + W,sin y, ¢)| cos 8, , 

92 

| @ = [e~? (F, cos ¥, @ + PF, sin V; 2) ( 
| +e? (Fs cos v9 p + Fy sin y, p)] cos 0, 


schreiben, worin von den acht Konstanten W,, F nur vier frei wahlbar sind. Die tbrigen lassen 


sich mit Hilfe der aus den Gleichungen (46), (50), (51) abgeleiteten Beziehungen berechnen: 


| W, =(W,+W,) =—K[(Fi+ Fn +i(i— Fy] 
=—K(yF, +4, Ff); 
| W, =—K iy, Fo—¥1 Fy) > (53) 


| Wy = K (x2 Fs + Yo Fi)» 

| W, = K (x2 Fy — Yo Fs) - 

Jetzt kénnen alle SchnittgroBen angegeben werden. Die Schnittkrafte erhalt man mit den 
; Gleichungen (6), (29), (31), (32) aus den Funktionen w, @ (38) und die Verschiebungen ergeben 
- sich aus den Gleichungen (13), (17) zu 


Wht i r ri es 0”) Ee | dé. | 


Dy (Lp) , % | 
— 


QoT 2 i ee NE PET ao dd . | 
Do U — #) Le ae 3)? 00 [| (5; (1 — p?) ey 


\ 


| 13 


(54) 
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8. Die isotrope Schale. Fiir den Sonderfall der isotropen Schale werden entsprechend den Glei- 
chungen (5), (34), (37) 
6, =0=1, ps=p=1, D7 =D, BeBe 


” 
Dann kann aber Gleichung (41) in der Form 


(lint)! + &8 (1 F Gym)? = (55) | 


geschrieben werden. Daraus erhalt man mit den Hilfsgréfen 


Bera see 1 tl eee 472 56) 
a= V Vise a +( Fy G0) a ees ( 


die Werte 


a (57) | 
= 1 /V 1 
= — 1 ( — 409 €&) —l1—e 
rs iz | ( ney ae: 5 © 00 2 ea) 
AuBerdem folgt aus den Gleichungen (48) 
Als ties ally (58) | 


(2) (2) 
so daB die Gleichungen (53) weitgehend vereinfacht werden. 
Der Wert ¢ der Gleichung (42) nimmt die Form 


e = |/RO= (59) 


Q) A? t 


an und liegt bei Stahlbetonschalen im allgemeinen zwischen ¢= 1,5 fiir kurze Schalen und ¢=6,0 
fiir lange Schalen. 


9. Die in y-Richtung durch Rippen verstarkte Schale. Sehr haufig miissen bei weitgespannten 
Schalen in der g-Richtung Verstarkungsrippen angeordnet werden, um die auftretenden Biege- 
momente aufnehmen zu kénnen oder um die Knicksicherheit der Schale zu vergréBern. Schon 


Dischinger? bemerkt, daB bei einem nicht zu groBem Rippenabstand das Tragwerk als eine aniso- | 


trope Schale aufgefaBt werden kann, die sich durch eine grofe Biegesteifigkeit B, auszeichnet. 
Die Biegesteifigkeiten B,, B,,, bleiben im Verhaltnis so klein, daB sie gegeniiber B, vernachlassigt 
werden kénnen. Die Dehnungssteifigkeit D, vergréBert sich durch die Anordnung von Rippen 
nur geringfiigig und naherungsweise wird daher D, =D, =D,,, gesetzt. 


Mit diesen Naherungen kénnen nach (34), (37) B=0, 6=1 gesetzt werden, so daB aus Glei- 
chung (41) die Gleichung 


mt (1 — mt)? + 28 (1 gym)? = 0 (60) 
folgt, deren Lésung mit den HilfsgréBen 


1 1 = i= 1 = 
a= Vn —earte (L Fz) tz0—#e (61) 


>) 


die Werte 
ie aes 1 1 2 ii Pi 3 l 
> ZV ta) +(keeasa) + baa 
! 1 2 1 = (62) 
aa a - : ; 
seas ebay -lad 
ergibt. 


10. Die Kreiszylinderschale. Alle bisher gebrachten Gleichungen bleiben nach dem Grenziiber- 


gang 0)—> 0 baw. r0)—> ry fiir die Kreiszylinderschale giltig. Bei Isotropie ergibt sich dabei die von | 


Zerna* gefundene Form. 


1 Siehe FuBnote 3 von S. 164. 
2 W. Zerna, Ing.-Arch. 20 (9a2) Siro oie 
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11. Zusammenfassung. Bei der Auflésung der Gleichungssysteme von doppelt gekriimmten 
Schalen treten im allgemeinen so grofe mathematische Schwierigkeiten auf, daB eine Berechnung 
der Schnittkrafte nicht méglich ist. In der vorliegenden Arbeit wird gezeigt, daB bei Einfiihrung 
einiger Naherungen der Spannungszustand einer Rotationsschale mit flacher, kreisférmiger Er- 
zeugender ermittelt werden kann. Sowohl bei der Aufstellung der Gleichgewichtsbedingungen als 
auch beim Elastizitatsgesetz werden in allen Gleichungen diejenigen Glieder, die die Schalenstarke 
enthalten, vernachlassigt, soweit sie neben anderen Gliedern vorkommen. Diese Naherung ist nicht 
gréber als die schon seit jeher bei Schalenberechnungen eingefiihrten Vereinfachungen des Elastizi- 
tatsgesetzes. Bei flachen Schalen kénnen in den Gleichungen weitere Glieder vernachlassigt werden, 
und es gelingt, das Problem auf zwei homogene partielle Differentialgleichungen vierter Ordnung 
mit konstanten Koeffizienten fiir eine Spannungsfunktion ® und die Normalenverschiebung w 
zurickzufiihren. Die charakteristische Gleichung wird fiir eine beliebige anisotrope Schale an- 
gegeben. Fiir die Sonderfalle der isotropen Schale und der in der Rotationsrichtung durch Rippen 
verstarkten Schale kénnen die Wurzeln der charakteristischen Gleichung als einfache Ausdriicke 
geschrieben werden. Die Schnittkrafte und Verschiebungen der Schale werden aus einer Spannungs- 
funktion und der Normalenverschiebung berechnet. 


(Eingegangen am 28. Juli 1956) 
Anschrift des Verfassers: Dr.-Ing. R. Windels, Hamburg 39, Krohnskamp 74. 
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Die Matrizentheorie der Statik 
Von J. H. Argyris 


1. Einleitung. Die folgende Arbeit! bringt einen Uberblick und eine betrachtliche Erweiterung 
der Matrizentheorie der Statik der elastischen Tragwerke, die vom Verfasser in einer friheren 
Abhandlung? entwickelt worden ist. ; 

Wir sind uns schon seit einigen Jahren bewuBt, daB keine der gewohnlichen statischen Methoden 
wirklich geeignet ist, die Spannungsverteilung und die Nachgiebigkeitsmatrizen der hochgradig 
statisch unbestimmten Systeme der modernen Luftfahrtkonstruktionen zu bestimmen. Ahnliche 
Schwierigkeiten treten auch in anderen Anwendungsgebieten der Statik auf. Die Iterationsver- 
fahren kénnen in gewissen Fallen niitzlich sein, sind aber im allgemeinen zu langwierig und haben 
sich nicht bei den membran- und schalenférmigen Tragwerken der Luftfahrt bewahrt. 

Diese Schwierigkeiten kénnen wir mit der Matrizenformulierung der Statik in Verbindung mit 
dem elektronischen Digitalautomaten iiberwinden. Die Matrizenformulierung erlaubt es nicht 
nur, die Rechnungen viel iibersichtlicher zu gestalten, sondern ist auch die ideale Schreibweise 
fiir den Digitalautomaten’. AufSerdem sind die theoretischen Ableitungen der Matrizentheorie so 
durchsichtig und elegant, daB neue und praktisch wertvolle Beziehungen, die in der gewéhnlichen 
Schreibweise unméglich oder nur schwierig erkennbar waren, sich jetzt sehr einfach ergeben. 

Die Theorie wird hier zugleich fiir das Kraft- und Deformationsverfahren entwickelt. Es zeigt 
sich, daB die beiden Methoden dual sind: jede Beziehung in dem einen Verfahren hat eine ent- 
sprechende Beziehung in dem anderen Verfahren, die sich durch eine einfache ,,Ubersetzung™ 
ergibt. Wir verweisen auch auf die frithere Arbeit”, wo die Vor- und Nachteile der beiden Methoden 
eingehend behandelt sind. Fiir kontinuierliche Tragwerke wie Flugzeugschalen ist das Kraftver- 
fahren im allgemeinen vorzuziehen. 

Im Abschnitt 2 werden die Einheitslast- und Einheitsverschiebungsgesetze, die eine besonders 
elegante Ableitung des Kraft- und Deformationsverfahrens erméglichen, wiedergegeben. Wir be- 
schranken uns hier auf kleine Verschiebungen und Verzerrungen, fiir die die gewéhnlichen Gleich- 
gewichtsbedingungen und Verzerrungsausdriicke gelten. Andererseits kann das Elastizitatsgesetz 
nichtlinear sein. Fiir Tragwerke nehmen die Einheitsgesetze eine besonders einfache Form ein. 

Der Abschnitt 3 bringt eine erweiterte Formulierung des Einheitslast- und Einheitsverschie- 
bungsgesetzes fiir den speziellen Fall von linearelastischen Systemen. AnschlieBend wird die Nach- 
giebigkeit und Steifigkeit eines Tragwerkes eingehend behandelt. Die eigentliche statische Be- 
rechnung von Systemen nach dem Kraft- und Deformationsverfahren wird danach im Abschnitt 4 
entwickelt. Alle Belastungsarten sowie Temperatureinfliisse sind beriicksichtigt. Es zeigt sich, 
daB wir die Nachgiebigkeit (oder Steifigkeit) des Tragwerks auBerst einfach als Endergebnis der 
Spannungsermittlung erhalten kénnen. Damit haben wir alle elastischen Unterlagen um auch 
die dynamischen Berechnungen und insbesondere die Eigenschwingungsuntersuchungen durch- 
zufiihren. Auf alle Falle ist es unnétig, die Nachgiebigkeit oder Steifigkeit ge- 
trennt von der statischen Berechnung zu bestimmen. Mit den hier gegebenen Methoden 
ist es méglich, in der Praxis die Untersuchungen von komplizierten Systemen wie Flugzeugfliigeln®* 
und Riimpfen systematisch fiir alle statischen und dynamischen Falle durchzufihren. Selbst- 
verstandlich ist auch die Behandlung von Balken- und Rahmenproblemen auch bei Beriicksich- 
tigung der Schubverformung ein einfaches Problem. Es ist beachtenswert, da® wir in allen stati- 
schen Berechnungen nur drei prinzipiell einfache Matrizen und eine Spaltenmatrix der Lasten 
bendtigen. 

Die allgemeine Theorie von Abschnitt 4 kann auf jedes Tragwerk, also auch auf eines mit 
Ausschnitten angewandt werden. Die praktische Berechnung derartiger Systeme wird aber doch 
im allgemeinen komplizierter und weniger iibersichtlich als fiir das entsprechende Tragwerk ohne 
Ausschnitte. Dieses vermeidet eine besondere Methode (Abschnitt 5), die schon in der vorerwahnten 
Arbeit? erschien, und die es erméglicht, die Spannungsverteilung des Tragwerks mit Ausschnitten 


* Erweiterte Fassung des Vortrages des Verfassers an der 1956 Pfingsttagung der Gesellschaft fiir an- 
gewandte Mathematik und Mechanik in Stuttgart. 

* J. H. Argyris, Aircraft Engineering, 26 (1954) S. 347, 383 und 27 (1955) S. 42, 80, 125, 145. Siehe 
auch J. H. Argyris u. S. Kelsey, Wissenschaftliche Gesellschaft fiir Luftfahrt, Jahrbuch 1956. 

3 P. M. Hunt, Aircraft Engineering, 28 (1956) S. 70, 111, 155. 
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von der statischen Berechnung des kontinuierlichen Systems abzuleiten. Das Verfahren ist ver- 
bliffend einfach und eignet sich auch fiir Fachwerke und Rahmen, deren regularer Aufbau nur 
durch einige fehlende Glieder gestért ist. 

In der Praxis kommt es immer wieder vor, daB einzelne Teile eines Tragwerks nach der Voll- 
endung der statischen Berechnung modifiziert werden. Um eine Wiederholung der ganzen elastischen 
Untersuchung zu vermeiden, wird im Abschnitt 6 eine Methode hergeleitet, um die Spannungs- 
verteilung des modifizierten Systems aus der des urspriinglichen Systems zu ermitteln. Dieses 
neue Verfahren ist eine Verallgemeinerung der Methode von Abschnitt 5 und in den meisten Fallen 
viel kiirzer als die direkte Berechnung des modifizierten 
Tragwerks. jeeenng 

Arbeiten itiber die Matrizentheorie der Statik insbeson- 
dere der Kraftmethode sind schon friiher! veréffentlicht 
worden. Es scheint aber, da keine dieser Abhandlungen 
so allgemein und doch so einfach wie die vorliegende 
Theorie ist. 


2. Allgemeine Grundlagen. a) Einfiihrung. Wir be- 
handeln in diesem Abschnitt die Grund gesetze der Matrizen- 
theorie der Statik. Wie erwahnt, beschranken wir uns auf 
den Fall kleiner Verschiebungen und Verzerrungen. Das 
Elastizitatsgesetz kann aber nichtlinear sein, so lange es 
eindeutig ist (Abb. 1). Um die Dualitat des Kraft- und 
Deformationsverfahrens besonders klar hervorzuheben, 
fihren wir die Gesetze und deren Ableitungen in zwei 
entsprechenden Spalten, und zwar links immer das Kraftverfahren, auf jeder Seite vor. Das 
zu jeder Beziehung duale Gesetz lat sich unter Benutzung folgenden ,,Worterbuches“ leicht 
ableiten. 


Dehnung 


Abb, 1. Nichtlineares Elastizitatsgesetz. 


Kraftverfahren Deformationsverfahren 


Kraft (Spannung) Verschiebung (Verzerrung) 


| 
| 
Verschiebung (Verzerrung) | Kraft (Spannung) 
| 
| 


apr F Verschiebung Ae ‘ Kraft 
Nachgiebigkeit = a errr Steifigkeit = Wate chichane 


Die Dualitat von Nachgiebigkeit und Steifigkeit wird schon von C. B. Biezeno und R. Grammel? 
eingehend behandelt. 

Der Begriff der Kraft ist hier allgemein gefaBt. Damit kann er sich auf eine Krafigruppe oder 
ein Moment beziehen. Ahnlich sind eine Verschiebungsgruppe oder Verdrehung auch in der Bezeich- 
nung Verschiebung enthalten. 


In einem Kontinuum sind die Spannungen und Verzerrungen eines Elementes dx dy dz als 
Spaltenmatrizen 


C5 — (Oz, Ovy O,, Oxy Ove ip und € = Lone Evy Exz Exy Eye ey (1) 
geschrieben. Man beachte, daB {+--+ } immer eine Spaltenmatrix bedeutet. Die allgemeine Be- 
zeichnung von Matrizen erfolgt mit eckigen Klammern [:**::- | 


In den folgenden Untersuchungen, die sich auf Tragwerke beziehen, sind die Begriffe Spannung 
bzw. Verzerrung allgemein zu verstehen und schlieSen auch innere Krafte, Momente usw. bzw. 
Langenanderungen, Drehungen usw. ein. Die Bezeichnung der Spannungen und Verzerrungen in 
Tragwerken erfolgt mit S, v. 


b) Die Einheitsgesetze 


Das Einheitslastgesetz | Das Einheitsverschiebungsgesetz 
Wir bezeichnen mit | Wir bezeichnen mit 
€ den Verzerrungszustand und mit r eine asso- | 6 den Spannungszustand und mit R eine asso- 
ziierte Verschiebung eines Kérpers infolge ge-  ziierte Kraft in einem Korper infolge gegebener 
gebener Lasten, Temperaturdehnungen oder ir- | Verschiebungen, Temperaturspannungen oder 
gendwelcher eingepragter Verzerrungen, wie sie | irgendwelcher eingepragter Zwangsspannungen, 


1 B., Langefors, Aeronautical Sciences, 19 (1952), S. 451; H. Falkenheiner, La Recherche Aeronautique 17 


1950). : 
e fs B. Biezeno u. R. Grammel, Technische Dynamik, Kap. II, Ziff. 9 und 10, 1, Aufl, Berlin 1939, 
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z. B. aus Fertigungsfehlern oder Lagerverschie- 
bungen entstehen! (Abb. 2a). 

Die Beziehung zwischen € und r ist eine kine- 
matische und kann von den Vertraglichkeits- 
bedingungen abgeleitet werden. Im folgenden 
benennen wir die Vertraglichkeitsbedingungen 
als kinematische Vertraglichkeitsbedingungen. 


Zustand 6, 
Abb. 2a. Das Einheitslastgesetz. 


Wir fiihren auch ein: 


6, ein (virtueller) Spannungszustand infolge 
einer Kraft R = 1 in Richtung von r, der nur 
den statischen, aber nicht den _ kinemati- 
schen Vertraglichkeitsbedingungen zu geniigen 
braucht. Der Spannungszustand ist so gewahlt, 
daB keine Verschiebungen des Zustandes € in 
Richtung der mit R = 1 im Gleichgewicht be- 
findenden Kraften stattfinden (Abb. 2a). 

Wir bezeichnen o, als einen statisch vertrig- 
lichen Spannungszustand. In allen folgenden 
Entwicklungen bezichen sich iiberstrichene Gré- 
Ben auf statisch vertragliche Spannungsvertei- 
lungen. 

Das Einheitslastgesetz ersetzt nun die kine- 
matische Beziehung zwischen € und r unter Be- 
nutzung der statisch vertraglichen 6, ; es schreibt 
sich? 


le7 fo, eaV. (2a) 


wo sich die Integration tiber das ganze Volumen 
erstreckt®, Gleichung (2a) 1aBt sich am einfach- 
sten vom Prinzip der virtuellen Krafte ab- 
leiten. 


wie sie z. B. aus Fertigungsfehlern oder Lager- 
verschiebungen entstehen! (Abb. 2b). 

Die Beziehung zwischen o und R ist eine 
statische und kann von den Gleichgewichts- 
bedingungen abgeleitet werden. Im folgenden 
benennen wir die Gleichgewichtsbedingungen 
als statische Vertraglichkeitsbedingungen. 


{’ x 


Zustand O 


Zustand €; 
Abb. 2b. Das Einheitsverschiebungsgesetz. 


Wir fiihren auch ein: 


€, ein (virtueller) Verzerrungszustand infolge 
einer Verschiebung r = 1 in Richtung von R, 
der nur den kinematischen, aber nicht den stati- 
schen Vertraglichkeitsbedingungen zu geniigen 
braucht. Der Verzerrungszustand ist so ge- 
wahlt, daB keine Krafte des Zustandes o in 
Richtung der mit r = 1 assoziierten Verschie- 
bungen wirken (Abb. 2b). 

Wir bezeichnen €, als einen kinematisch ver- 
traglichen Verzerrungszustand. In allen folgen- 
den Entwicklungen beziehen sich unterstrichene 
GréBen auf kinematisch vertrigliche Verzer- 
rungsverteilungen. 

Das Einheitsverschiebungsgesetz ersetzt nun 
die statische Beziehung zwischen o und R unter 
Benutzung der kinematisch vertriaglichen €,; 
es schreibt sich? 

1h el aide (2b) 
wo sich die Integration iiber das ganze Volumen 
erstreckt®, Gleichung (2b) 1a8t sich am einfach- 


sten vom Prinzip der virtuellen Verschiebungen 
ableiten. 


* Wir gebrauchen im folgenden die kollektive Bezeichnung: gegebene Beanspruchung fiir diese verschiedene 


Belastungsarten. 
2 Vel. FuBnote 2 S. 174. 


° Ein Apostroph bezeichnet die transponierte oder gespiegelte Matrix. 


 zerrungszustand ist; es geniigt daB 


\ 
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Da das Gesetz (2a) nur kinematische, aber | 
nicht statische Beziehungen ersetzt, ist eigent-_ 
lich nicht notwendig, daB € ein wahrer Ver- 


€=€ 


ist, wo € ein kinematisch vertraglicher Verzer- 


-rungszustand infolge der gegebenen Beanspru- 


chung ist. 


Da das Gesetz (2b) nur statische, aber nicht 
kinematische Beziehungen ersetzt, ist eigentlich 
nicht notwendig, das o ein wahrer Spannungs- 
zustand ist; es geniigt das 


o=«a 


ist, Wo G ein statisch vertraglicher Spannungs- 
zustand infolge der gegebenen Beanspruchung 
ist. 


c) Anwendung auf Tragwerke. Wir interessieren uns hier besonders fiir die Anwendung 
der Matrizentheorie der Statik auf Tragwerke. Jedes Tragwerk wird fiir die praktische Berechnung 
in eine endliche Anzahl von Elementen zerlegt. Wir setzen voraus, daB die Spannungen und Ver- 
zerrungen in jedem Elemente vollstandig bestimmt sind, wenn gewisse Krafte oder Spannungen S 
und Verschiebungen oder Verzerrungen v an den Randern des Elements bekannt sind. Das ein- 
fachste Beispiel eines Elementes ist ein Stab in einem idealen Fachwerk; hier geniigt es, eine 
Spannung oder Verzerrung in jedem Elemente vorzuschreiben, um das ganze Spannungs- oder 
Verzerrungsfeld zu kennen. In einem kontinuierlichen Tragwerk wie einem Flugzeugfliigel unter- 
teilen wir das System mit einem Netz von Linien und benennen wieder die Teile des Tragwerkes, 
die zwischen zwei sich schneidenden Paaren von benachbarten Netzlinien liegen, als Elemente. 
Die Schnittpunkte der Netzlinien! werden als Knotenpunkte bezeichnet. Es wird im allgemeinen 
angenommen, da die 4uBeren Krafte an den Knotenpunkten wirken?, und daB die Langsspannun- 
gen linear? zwischen benachbarten Punkten variieren. Wir kénnen hier nicht auf die naheren 
Einzelheiten des Verfahrens eingehen, verweisen aber auf die friihere Arbeit’. 


Wir betrachten jetzt ein Tragwerk, das aus / Elementen besteht. Aufere Krafte R,::-R,, oder 
Verschiebungen r,:*-r,, kbnnen an m Knotenpunkten des Systems vorgeschrieben sein. 


Die folgenden Bezeichnungen werden eingefiihrt: 


R={R,'R,,}  AuBere Krafte, 

Sa Ores} Spannungen oder Krafte an den / Elementen infolge R, 3 
iets at came Verschiebungen in Richtung von R, | (3) 
© = {v,°"-0,} Verzerrungen der / Elemente infolge r. 


Abb. 3b. Beispiel fiir die Berechnung 
einer b-Matrix. 


Abb. 3a. Beispicl fiir die Berechnung einer @-Matrix, 


Wenn mehr als eine Spannung oder Verzerrung in den Elementen vorgeschrieben ist, bezeichnen 

S,, v, usw. Untermatrizen. Der Verschiebungszustand (r, v) und der Spannungszustand 
. . i 

(R, S) kénnen unabhangig von einander sein. Man beachte, daB R, S, ” und v Spalten- 


matrizen sind. , 


Das Einheitslastgesetz 
WennS ein statisch vertraglicher Spannungs- 
zustand, der den Kraften R entspricht, ist, dann 
kénnen wir immer 


S=bR (4a) 


1 Die Netzlinien brauchen nicht orthogonal zu sein. 


Das Einheitsverschiebungsgesetz 

Wenn v ein kinematisch vertriaglicher Ver- 
zerrungszustand, der Verschiecbungen Yr ent- 
spricht, ist, dann kénnen wir immer 


v=ar (4b) 


2 Diese Annahmen sind nicht notwendig, siehe Abschnitt 4 a). 


3 Vel. FuBnote 2 von S, 174. 
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setzen, wo D eine rechteckige Matrize, die sich 
lediglich aus statischen Beziehungen berechnet, 
ist (siehe auch die Definition von 6,). Man be- 
achte, daB (4a) auch fiir nichtlineare Systeme 
giiltig ist. 

Als Beispiel einer b Matrix betrachten wir 
den statischen vertraglichen Spannungszustand 


in Abb. 3a. Hier ist 
R, R, 

=I Om Vial 

==ll 0 2, 

| 0 0 3 

2 1 4. 

a 0 0 5 
b = 0 1 6 (5a) 

2 0 7 

0 0 8 

0 0 9 

i —y2| 1 


Uber jeder Spalte ist die zugehérige Kraft an- 
gegeben. Die Zahl in kursivem Druck neben 
jeder Spalte bezieht sich auf die Numerierung 
der Stabe in Abb. 3. 

Wir kénnen nun leicht zeigen, daB das Ein- 
heitslastgesetz (2a) fiir Tragwerke in folgender 
einfacher Form geschrieben werden kann: 


r=b'v. (6a) 


An Stelle von © kann man einen kinematisch 
vertraglichen Zustand WU einfiithren und erhalt 
unter Benutzung von (6b) 


setzen, wo @ eine rechteckige Matrize, die sich 
lediglich aus kinematischen Beziehungen be- 
rechnet, ist (siehe auch die Definition von €,). 
Man beachte, daB (4b) auch fiir nichtlineare 
Systeme giiltig ist. 

Als Beispiel einer @ Matrix betrachten wir | 
den kinematisch vertraglichen Verzerrungs- 
zustand in Abb. 3b. Hier ist 


i Up 
[ee O Ossetian! 
0 0 2 
0 0 3 
0 0 4 
1 0 5 

a= 0 1 6 (5b) 

1/y2 0 i 
ae ete 
0 1/y2 9 
0 0 10 


Uber jeder Spalte ist die zugehérige Verschie- 
bung angegeben. Die Zahl in kursivem Druck 
neben jeder Spalte bezieht sich auf die Nume- 
rierung der Stabe in Abb. 3. 

Wir kénnen nun leicht zeigen, daB das Ein- 
heitsverschiebungsgesetz fiir Tragwerke in fol- 
gender einfacher Form geschrieben werden 
kann: 


R=a S. (6b) 
An Stelle von S kann man einen statisch ver- 


traglichen Zustand S einfiihren und erhalt unter 
Benutzung von (6a) 


ba=E=a'b, (7) 
wo E eine Einheitsmatrix ist. Die interessante Beziehung (7) laBt sich auch direkt herleiten. 
Selbstverstandlich setzt (7) voraus, daB die Anzahl und Richtungen der Krafte R in (6b) und 


Verschiebungen Y in (6a) ttbereinstimmen. Wir bestatigen leicht, da die Gleichungen (5) die 
Beziehung (7) erfiillen. 


Statisch bestimmte Tragwerke Kinematisch bestimmte Tragwerke 
Wenn S die wahren Spannungen infolge R 


sind, dann gilt fiir statisch bestimmte Trag- 
werke immer 


Wenn © die wahren Verzerrungen infolge 1 
sind, dann gilt fiir kinematisch bestimmte Trag- 
werke immer 


S=obR, 


fe li 


(8a) — (8b) 


7 
7 


Abb. 4a. Statisch bestimmtes System. Abb. 4b. Kinematisch bestimmtes System. 


wo die Matrix b sich lediglich aus den Gleich- 


wo die Matrix @ sich lediglich aus den Vertrag- 
gewichtsbedingungen ermitteln lat. Demnach 


lichkeitsbedingungen ermitteln la8t. Demnach 
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ist (egip | ist a=a. 


3 Pa : ; | 
Kin Beispiel eines statisch bestimmten Trag- 


Kin Beispiel eines kinematisch bestimmten Trag- 
werkes zeigt Abb. 4a. 


werkes zeigt Abb. 4b. 
3. Anwendung auf linearelastische Systeme. Wir beschranken uns in diesem und den folgenden 
Abschnitten auf linearelastische Systeme, fiir die also das Hookesche Gesetz gilt. 


a) Die Einheitsgesetze. 
Das Einheitslastgesetz | Das Einheitsverschiebungsgesetz 
Die Anwendung des Prinzips der virtuellen Die Anwendung des Prinzips der virtuellen 
Krafte fiihrt bei linearelastischen Systemen auf | Verschiebungen fiihrt bei linearelastischen Sy- 
eine zweite Formulierung des LEinheitslast- | stemen auf eine zweite Formulierung des Kin- 
gesetzes, die der von (2a) aquivalent ist. Wir | heitsverschiebungsgesetzes, die der (2b) aqui- 
finden valent ist. Wir finden 
iT | C,€ay ==) 6, 6d , (9a) Io Rie fecuy elec dV, (9b) 
wo der erste Ausdruck fiir die wahre Verschie- | wo der erste Ausdruck fiir die wahre Kraft R 
bung r in Abschnitt 2’ a) erlautert ist und € ein | in Abschnitt 2’ a) erlautert ist und o ein Span- 
Verzerrungszustand, der einem statisch vertrag-  nungszustand, der einem kinematisch vertrag- 
lichen Spannungszustand o infolge der gege- | lichen Verschiebungszustand € infolge der gege- 
benen Beanspruchung entspricht, ist [man kann | benen Beanspruchung entspricht, ist [man kann 
demnach € in jedem geeigneten statisch aqui- demnach o in jedem geeigneten kinematisch 


valenten — z. B. einem statisch bestimmten — | Aquivalenten — z. B. einem kinematisch be- 
Teilsystem finden (Abb. 5a)] und 6, der wahre | stimmten — Teilsystem finden (Abb. 5b)| und €, 
Spannungszustand infolge R = | ist. der wahre Verzerrungszustand infolge r = 1 ist. 

Wie im Falle der nichtlinearen Systeme setzt Wie im Falle der nichtlinearen Systeme setzt 


(9a) voraus, daB keine mit r assoziierten Ver- | (9b) voraus, daB keine mit R im Gleichgewicht 
schiebungen in Richtung der mit R= 1 im | befindlichen Krafte in Richtung der mit r = 1 
Gleichgewicht befindlichen Krafte stattfinden. | assoziierten Verschiebungen wirken. 


ia 


Zustand é 


| 
| 


Zustand O,; 


Abb. 5a. Das Einheitslastgesetz fiir lineare Systeme. | Abb. 5b, Das Einheitsverschiebungsgesetz fiir lineare Systeme. 


Zustand &, 


Die Gleichungen (9a) und (9b) zeigen, daf es fiir lineare Systeme méglich ist, die wahren und 
vertraglichen Zustande zu vertauschen. Selbstverstandlich ist es zulassig, in (9a) und (9b) nur 
wahre Zustande zu benutzen und 

lor fojedV (Oo)n| =n (se Od V1 (9d) 
zu schreiben. Die Einfiithrung von statisch und kinematischen vertraglichen (virtuellen) Zustanden 
bedeutet aber, wie bekannt, eine grofe Vereinfachung der Theorie und der praktischen Berech- 
nungen. 
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b) Anwendung auf Tragwerke. 


Das Einheitslastgesetz 
Wenn S wieder die Spaltenmatrix der Span- 
nungen infolge von Lasten R bedeutet, dann 
gilt fiir alle lineare Tragwerke 


SS U0lts (8a) 


Das Einheitslastgesetz (2a) und (9a) nimmt nun 
fiir lineare Tragwerke die gleichwertigen For- 
men 


r=bv=b'v=b'b (10a) 


ein, wo © der Definition von € entspricht. Die 
Verschiebungen 7 in Richtung der R kénnen 
selbstverstandlich wie in (9a) infolge einer belie- 
bigen Beanspruchung entstehen. Fiir die mit 
R im Gleichgewicht befindlichen Krafte gilt die 
in 3a) erwahnte Beschrankung. 


Das Einheitsverschiebungsgesetz 

Wenn U wieder die Spaltenmatrix der Ver- 
zerrungen infolge von Verschiebungen Yr he- 
deutet, dann gilt fiir alle lineare Tragwerke 


O =i (8b) 


Das Einheitsverschiebungsgesetz (2b) und (9b) 
nimmt nun fir lineare Tragwerke die gleich- 
wertigen Formen 

R=a S=asS—eas (10b) 
ein, wo S der Definition von o entspricht. Die 
Krafte R in Richtung von? kénnen selbstver- 
standlich wie in (9b) infolge einer beliebigen 
Beanspruchung entstehen. Fiir die mit 1 asso- 
ziierten Verschiebungen gilt die in 3 a) erwahnte 
Beschrankung. 


Wenn jetzt die Krafte R in (10b) und die Verschiebungen Yr in (10a) an den gleichen Knoten- 
punkten und Richtungen wirken, dann gilt [siehe (7)]. 


ba=bva=E=a b=ab. 


(11) 


Die Nachgiebigkeit des Tragwerks 


Die Steifigkeit des Tragwerks 


Wir bezeichnen die Spannung und Verzerrung eines Elementes p in einem Tragwerk mit S,, 


und U,. Es bestehen folgende Beziehungen: 
©, =fp>S>> (12a) 
wo f, als die Nachgiebigkeit des Elementes p 


bezeichnet wird. Wenn wir an einem Element 
mehr als eine Spannung und Verzerrung vor- 
schreiben, sind S,,, v,, {,, Matrizen. Die Matrix 
[, ist immer symmetrisch (Maxwellsche Rezi- 
prozitatssatz}). 


S,= k,v,. (12b) 
wo k, als die Steifigkeit des Elementes p be- 
zeichnet wird. Wenn wir an einem Element 
mehr als eine Verzerrung vorschreiben, sind 
©; So k,, Matrizen. Die Matrix k, ist immer 
symmetrisch (der zum Maxwellschen duale Rezi- 
prozitatssatz!). 


Das einfachste Beispiel fiir ik und k, ist durch einen Stab in einem idealen Fachwerk gegeben. 
Wenn hier Ss die Stabkraft und OF die Langenanderung des Stabes sind, ist 


l 


ib “Tyee (13a) 


_EA 
an ges 


k (13b) 


wo E, A, | augenfallige Bezeichnungen sind. Die Nachgiebigkeit und Steifigkeit von komplizier- 
teren Elementen ist in der friiheren Arbeit*, die auch die Schubverformung beriicksichtigt, ein- 
gehend untersucht. Im allgemeinen sind die Definitionen der Spannung S p und der Verzerrung U,, 
der Elemente im Kraft- und Deformationsverfahren verschieden. 


Unter Benutzung von (3) kénnen wir nun Unter Benutzung von (3) kénnen wir nun 
= fS (12c) S=kv (12d) 
setzen, wo setzen, wo 
? [MU Sarat 2 O- k, 0---------- 07 
Os. 0. 
f =| 0---07,0----0 (14a) = 0-20, 0 een (14b) 
& S 0 aS 0 
| 0---------- Of, | 0--------- Ok, _ 


' Vel. C. B. Biezeno und R. Grammel, a. a. O., Kap. II, Ziff. 9. 
* Siehe FuGBnote 2 von Seite 174. 
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eine diagonale Matrix ist; f wird als die Nach- 
giebigkeit der ] unverbundenen Elemente des 
Tragwerks bezeichnet. 


eine diagonale Matrix ist; Kk wird als die Steifig- 
keit der ] unverbundenen Elemente des Trag- 
werks bezeichnet. 


Wenn in (12c) und (12d) dieselben Spannungen (Krafte) und Verzerrungen vorgeschrieben 


sind, besteht die Beziehung 


fE=E=k?. (15) 


Anwendung der Gleichungen (10a) und (12c) 
ergibt fiir die m Verschiebungen Yr in Richtung 
und infolge der R 


r=b fbR=b' fbR=b fb=FR. (16a) 
Hierin ist 

F=b'fb=b'fb=b'fb (17a) 
die Nachgiebigkeit des Tragwerks in den vor- 
geschriebenen m Richtungen. Gleichung (17a) 
zeigt, daf} die quadratische Matrix symmetrisch 


ist, wie es auch der Maxwellsche Reziprozitats- 
satz vorschreibt!. 


Anwendung der Gleichungen (10b) und (12d) 
ergibt fiir die m Krafte R in Richtung und 
infolge der r 


R=ad kar=ad'kar=da kar=Kyr. (16b) 
Hierin ist 

K=adka=adka=a’'ka_ (17b) 
die Steifigkeit des Tragwerks in den vorge- 
schriebenen m Richtungen. Gleichung (17b) 


zeigt, da} die quadratische Matrix symmetrisch 
ist, wie es auch der zu Maxwell duale Rezipro- 


zitatssatz vorschreibt!. 


Die Beziehung (15) verallgemeinert sich in einem Tragwerk zu 


FR =E=KPF. (15a) 


Der Aufbau des Tragwerks aus den ver- 
schiedenen Elementen erfolgt in (17a) durch 
die Matrix b, die eine statische Beziehung aus- 
driickt. Damit kénnen wir den Aufbau (17a) 
als eine verallgemeinerte Serienschaltung von 
-Federn oder Elementen betrachten. Abb. 6a 
zeigt ein einfaches Beispiel einer Serienschal- 
tung. 


Ye a 


Abb. 6b. Beispiel der 
Parallelschaltung von 
Elementen. 


Kg ak 


a b 
Fy 35 Fy 


Der Aufbau des Tragwerks aus den ver- 
schiedenen Elementen erfolgt in (17b) durch 
die Matrix @, die eine kinematische Beziehung 
ausdriickt. Damit kénnen wir den Aufbau (17b) 
als eine verallgemeinerte Parallelschaltung von 
Federn oder Elementen betrachten. Abb. 6b 
zeigt ein einfaches Beispiel einer Parallelschal- 
tung. 


Abb. 6a. Beispiel der 


Z Serienschaltung von Ele- 

Yy menten, Die mit Schwarz 
bezeichneten Teile des 

Kragbalkens sind jeweils 


als starr angenommen. 
= F 


b UY 
——_ Y 
Ky = K 


Beispiel. Als Anwendung der Nachgiebigkeits- und Steifigkeitsmatrizen Ff und K betrachten 
wir die Eigenschwingungen eines Tragwerkes, wenn m Massen M, an den m Knotenpunkten wirken. 


Wir bezeichnen mit 


1 Vel. FuBnote 2 von S. 174. 
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die Diagonalmatrix der Massen. Die m Higenkreisfrequenzen @ und die Form r = 1, e'®! der Eigen 
schwingungen des Tragwerkes berechnen sich 


1 
entweder aus! | oder aus 


/FM——E|=0 | MK — a? E | =0 
| yy” | 


und (18a) (18b) 


und 
| [M21K—oP*E |r =9. 
Allgemeinere Massenverteilungen lassen sich ohne weiteres beriicksichtigen. Da die Nachgiebig- 


keit (Steifigkeit) sich nicht nur auf Krafte (Verschiebungen), sondern auch auf Momente (Ver- 
drehungen) bezieht, ist es leicht méglich, die Rotationstragheit in die Rechnungen einzusetzen. 


atc bey 
FM—-,E|r, =0. 


Die Untersuchungen von Schwingungen in Balken und Rahmen, deren F’ oder K (mit Beriick- 
sichtigung der Schubverformung) einfach zu ermitteln ist, ergeben interessante Ubungsbeispiele. 


4, Die Berechnung von Tragwerken. 


Wir untersuchen in diesem Abschnitt eine 
Reihe von Problemen in Tragwerken, deren 
Spannungen nicht durch statische Betrach- 
tungen allein ermittelt werden kénnen. Wir 
bezeichnen derartige Tragwerke als statisch un- 
bestimmt. Die tiberbestimmten GréBen heiBen 
statisch unbestimmte Krafte oder Spannungen. 
Abb. 7a zeigt ein zweifach statisch unbestimm- 
tes Fachwerk. 


Wir untersuchen in diesem Abschnitt eine 
Reihe von Problemen in Tragwerken, deren 
Verzerrungen nicht durch kinematische Be- 
trachtungen allein ermittelt werden kénnen. 
Wir bezeichnen derartige Tragwerke als kine- 
matisch unbestimmt. Die iiberbestimmten Gré- 
Ben heiBen kinematisch unbestimmte Verschie- 
bungen oder Verzerrungen. Abb. 7b zeigt ein 
sechsfach kinematisch unbestimmtes Fachwerk 
(r, und r, vorgeschrieben, U, bis U, unbekannt). 


Abb. 7a. Zweifach statisch unbestimmtes System. 


a) Problem I. 


In einem n-fach statisch unbestimmten? Trag- 
werk seien die m Krafte R gegeben (Abb. 7a). 
Ks sollen die statisch unbestimmte Kraft- oder 
Spannungsmatrix 


Wx eX. (19a) | 


und die Spannungen S bestimmt werden. 
Die vollstandige Kraftmatrix schreibt sich 


{RX}. (20a) | 


e760 Jay ie: 0! 0) =) sie. = elle 


Wir kénnen nun die Spannungen S in der Form 
sehen 


S=b,R+b,X=bR,  (2la) 


wo die Matrizen Bb, und b, aus statischen Be- | 
trachtungen berechnet werden und Db noch | 
unbekannt ist. Das Teilsystem, in dem b, 
ermittelt wird, benennen wir das_ statische 
Grundsystem. 


1 Selbstverstandlich ist die Determinantengleichung in (18 a) bzw. 


enthalten. 


Abb. 7a. Sechsfach kinematisch unbestimmtes System. 


In einem n-fach kinematisch unbestimmten? 
Tragwerk seien die m Verschieb-r vorgeschrie- 
ben (Abb. 7b). Es sollen die kinematisch un- 


_ bestimmte Verschiebungsmatrix 


©9116: 16) 61 s0. ae ome 


U ={U,U, U,}  (19b) 


und die Verzerrungen © bestimmt werden. 
Die vollstandige Verschiebungsmatrix 
schreibt sich 
AT URS (20b) 


Wir kénnen nun die Verzerrungen U in der 
Form sehen 


vo=ar+a,U =ar, (21b) 


wo die Matrizen @) und @, aus kinematischen 
Betrachtungen berechnet werden und @ noch 
unbekannt ist. Das Teilsystem, in dem dy 
ermittelt wird, benennen wir das kinematische 
Grundsystem. 


(18b) in der folgenden Matrizengleichung 


* Die Anzahl der statisch und kinematisch iiberbestimmten Gré®en eines Tragwerks ist selbstverstandlich 


im allgemeinen verschieden; vgl. auch Abb. 7a u. 7b. 
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Um die X zu ermitteln, wenden wir die 
Methode der Einheitslast (10a) auf die Vertrag- 
lichkeitsbedingung in Richtung von X an. Wir 
finden 


bv=b fS=b fb,R +b fb,xX —0 
oder 


X=—DOD,R, (22a) 


D=b,fb, wd D,=bifb,. (23a) 
Somit folgt aus (21a) und (22a) 
b=b,—b, D“D,. 


womit die Spannungen und Verzerrungen im 
Tragwerk bekannt sind. 


Um die Verschiebungen r in Richtung der | 


Krafte R zu bestimmen, benutzen wir wieder 


(10a) und erhalten 
r—bv=b fbR=FR. 
Wir kénnen fir b (das statisch vertragliche 


System) 6, einsetzen. Somit ist die Nachgiebig- 
keit des Systems 


F=f, — Di pp (26a) 


Hierin ist F,) = b, f by die Nachgiebigkeit des 


Grundsystems. 


rechnung des Systems nach dem Kraftverfahren 
nur die Matrizen b), b,, f und die Spalten- 
matrix R bendtigt. 


Beispiel. Fiir das zweifach statisch un- 
bestimmte Fachwerk in Abb. 7a soll die Ma- 
trix b und die Nachgiebigkeit F' unter den Be- 
dingungen bestimmt werden, daf alle Stabe die 
gleiche Nachgiebigkeit 1/E A besitzen. 

Die gewahlten statisch unbestimmten Gré- 
Ben sind in Abb. 7a angegeben. Die entspre- 
chenden Matrizen b, und 6b, folgen aus ele- 
mentaren statischen Betrachtungen. Es ist 


X = {Xy X;} ? 


ee XG 5%, 
ere 0 Ly Ue | 
es ey hie | 2 
0 0 |—1/y2 0, 3 
i coat | De Sea 
fo <0 ion eS 
og > =) _ ays aya 6 
{og 40. | 1 0 Z 
poe V2: ape 0 1 8 
Ore Oe 1 On=ia9 
OFk0-] 0 1 10 


— 
i) 
= 
pe) 

— 


(24a) | 


(25a) | 


Um die U zu ermitteln, wenden wir die 
Methode der Einheitsverschiebung (10b) auf 
die Gleichgewichtsbedingung in Richtung von 


U an. Wir finden 
aS=akv=akar+aka,U =9 


oder 
U=—C1C,r; (22b) 
wo 
C=G,koa, umd C,=a,ka,. (23b) 
Somit folgt aus (21b) und (22b) 
a=a,—a,C'C,, (24b) 


womit die Verzerrungen und Spannungen im 
Tragwerk bekannt sind. 

Um die Krafte R in Richtung der Verschie- 
bungen Y zu bestimmen, benutzen wir wieder 


(10b) und erhalten 


R=aS=adkar=Kr.  (25b) 


| Wir kénnen fiir @ (das kinematisch vertrag- 
_ liche System) Gy einsetzen. Somit ist die Stei- 


figkeit des Systems 
K = K,—C,C1C,. (26b) 
Hierin ist K, =a,ka, die Steifigkeit des 


| Grundsystems. 
Es ist besonders zu beachten, daB die Be- | 


Es ist besonders zu beachten, da die Be- 
rechnung des Systems nach dem Deformations- 
verfahren nur die Matrizen @), @,, k und die 
Spaltenmatrix Y bendtigt. 


Beispiel. Fiir das sechsfach kinematisch 
unbestimmte Fachwerk in Abb. 7b soll die 
Matrix @ und die Steifigkeit KA unter den Be- 
dingungen bestimmt werden, dab alle Stabe 
die gleiche Steifigkeit E.A/I besitzen. 

Die gewahlten kinematisch unbestimmten 
GréBen sind in Abb. 7b angegeben. Die ent- 
sprechenden Matrizen @, und @, folgen aus ele- 
mentaren kinematischen Betrachtungen. Es ist 


U = {U, U, Us; U, Us U;} » 
a, =4@ aus (5b) da Grundsysteme identisch, 


U, U, Us U, Us Us 
ere eee tame (or ( 1 
Ose lar 0 es O20) Choe re 
(Yale eat) 3 
(Pt ae 4 
be 00 f0m) 0 -=te 20 5 
Cam en e025 105 Oral Be 6 
iy2 0° 0 —1//2) 0 —=1jy2 | 7 
(ped nme 0ee 70 50 8 
C= 122 0- Saiy2 0° 5x9 
=o 0 0 1/2 0 —1/72 | 10 
(27b) 
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Hier bezieht sich die Zahl in kursivem Druck neben jeder Zeile auf die Numerierung der Stabe 
in Abb. 3 und Abb. 8. Uber jeder Spalte ist die zugehorige Kraft oder Verschiebung angegeben. 


Abb. 8. Numerierung der Stabe. 


Die Nachgiebigkeit der unverbundenen Ele- 
mente ist 
l 


f= Tea Ey, (28a) 


Die Steifigkeit der unverbundenen Elemente 
ist 


wo der Index 10 die Anzahl der Diagonalelemente in der EKinheitsmatrix angibt. 


Anwendung der Gleichung (24a) ergibt? 


R, R, 
rete See i) 
=e 31 2 
ele eee 5 
83 24 4 
; At 3 5 
Deets ey 21 16, an(29a) 
S12 3 2" 7 
22 esi Sales 
—2472 372] 9 
—28/2 —24/2 | 10 
Fir die Nachgiebigkeit des Tragwerks finden 
wir aus (26a) 
. I 397 136 | i 
eines a 
55 EA Abs op (30a) 


Wir bestatigen leicht, daB die Beziehungen b’ 


b) Problem II. 


In einem n-fach statisch unbestimmten Trag- 
werk seien die Verschiebungen Y vorgeschrieben. 
Es sollen X und S berechnet werden. Die 
Krafte R sind hier unbekannt, kénnen aber aus 
(25a) berechnet werden. Es ist 

R=F"'r. (31a) 
Damit ist dieses Problem auf Problem I zuriick- 
gefiihrt. 


c) Problem III. 


‘Wir bezeichnen als Anfangsverzerrungen H 


k==4h,,, (28b) 
Anwendung der Gleichung (24b) ergibt? 
| Te 
es Soil 
62> oF 21 2 
48 = 81 g 
81) 7 = 332 4 
, 48  — 81 5 | 
a= —-| —57 159 6 (29b) 
335 a 
45 V2 = 5525 7 
—3Iy2 isa} 2a) 8 
—48f2 81y2| 9 
12/2 —104/2 | 10 


Fur die Steifigkeit des Tragwerks finden wir 
aus (26b) 


93 —136 


EA 
K=—— (30b) 


—136 397 


a= E, und F K = E, erfillt sind. 


die auf die unverbundenen Elemente einge-— 


pragten Verzerrungen infolge von Temperatur- 


anderungen, Fertigungsfehlern, Lagerverschie- | 


bungen, Lasten auf Elementen, usw. Ein ein- 
faches Beispiel einer Anfangsverzerrung liegt 


In einem n-fach kinematisch unbestimmten 
Tragwerk seien die Krafte R vorgeschrieben. 
Es sollen U und v berechnet werden. Die Ver- — 
schiebungen Y sind hier unbekannt, kénnen © 
aber aus (26b) berechnet werden. Es ist | 


r=KOR. (31b) 
Damit ist dieses Problem auf Problem I zuriick- 
gefihrt. 


Wir bezeichnen als Zwangsspannungen J die 
auf die Elemente mit festgehaltenen Knoten- 


| punkten eingepragte Spannungen infolge von 


Temperaturanderungen, Fertigungsfehlern, La- 
gerverschiebungen, Lasten auf Elementen usw. 
Kin einfaches Beispiel einer Zwangsspannung 


* Aus Platzersparnisgriinden sind hier die Matrizen C, C1 und D, D~ nicht wiedergegeben. 
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in einem Fachwerk vor, dessen Stab p von der 
Lange |, einer Temperaturanderung O- 
gesetzt ist. Dann ist die Langenanderung des 
Stabes 


ee nO pe 


wo a, der Temperaturausdehnungskoeffizient 
ist. Wenn mehr als eine Verzerrung im Element 
vorgeschrieben ist, wird H,, eine entsprechende 
Spaltenmatrix! H,. 


Der Aufbau des Tragwerks aus den ver- 
zerrten unverbundenen Elementen wird im 
allgemeinen einen zusatzlichen Spannungszu- 
stand benétigen, wenn das System statisch 
unbestimmt ist. 


Wir untersuchen nun folgendes Problem. 
Es seien vorgeschrieben: 1) die Matrix der An- 
fangsverzerrungen der unverbundenen Elemente 


1S fee |: aoe H,} (32a) 
und 2) die m Krafte R =0 an den Knoten- 
punkten. 


Es sollen wieder X und S berechnet werden. 
Die gesamte Verzerrung der Elemente ist die 
Summe der elastischen und der eingepragten 
Anfangsverzerrungen. Somit ist 


vo=fb,X+H. (33a) 

Anwendung von (10a) und (33a) ergibt 
bv=DX+b,H=0 

und 

X=—D"b,H. (34a) 
Damit wird 

S=b,X =—b,D"b | H 

und (35a) 


o=—fb,D'b6,H+H. 

Die Verschiebungen? r der m Knotenpunkte 
infolge von H berechnen sich wieder aus (10a). 
Wir finden 

P= 0,0. b fi’, (36a) 
wo 6b, und 6 in Problem I definiert sind. Die 
zweite Formel in (36a) zeigt, da die Verschie- 
bungen Y eines Systems infolge von Anfangs- 
verzerrungen H ohne die entsprechende Span- 
nungsberechnung ermittelt werden kénnen; es 
geniigt, wenn die wahren Spannungen (also die 


Matrix b) infolge von Lasten R in Richtung— 


von Y bekannt sind. 


d) Problem IV. 


Es sind vorgeschrieben: 1) die Anfangsver- 
zerrungen 


H=(H,} 


1 Vel. FuBnote 2 von S. 374. 


| liegt in einem Fachwerk vor, dessen Stab p vom 
aus- | 


Querschnitt A, einer Temperaturanderung C, 
ausgesetzt ist. Dann ist die Zwangskraft im 
Stabe, wenn die Endknotenpunkte festgehalten 
sind 
Jp=—EA,«,9,, 

wo «, der Temperaturausdehnungskoeffizient 
ist. Wenn mehr als eine Zwangsspannung am 
Element vorgeschrieben ist, wird J, eine ent- 
sprechende Spaltenmatrix! J. 

Die Freilassung der festgehaltenen Knoten- 
punkte des Tragwerks unter Zwangsspannungen 
wird im allgemeinen einen zusatzlichen Ver- 
zerrungszustand hervorrufen, wenn das System 
kinematisch unbestimmt ist. 

Wir untersuchen nun folgendes Problem. Es 
seien vorgeschrieben: 1) die Matrix der Zwangs- 
spannungen in den Elementen mit festgehal- 
tenen Knotenpunkten. 


Pers HI J J;} (32h) 


und 2) die m Verschiebungen Y = 0 der Knoten- 
punkte. 

Es sollen wieder U und v berechnet werden. 
Die gesamte Spannung in den Elementen ist die 
Summe der elastischen und der eingepragten 
Zwangsspannungen. Somit ist 


S=ka,U+J. (33b) 
Anwendung von (10b) und (33b) ergibt 
a,S=CU+a,J=0 
und : 
U=—C'aJ. (34b) 
Damit wird 
vo=a,U =—a,C a J 
und (35b) 


S=—ka,C1a,J+J. 


Die Krafte? R an den m Knotenpunkten infolge J 
berechnen sich wieder aus (10b). Wir finden 


R=aS=a'J, (36b) 


wo @, und @ in Problem I definiert sind. Die 
zweite Formel in (36b) zeigt, daB die Krafte R 
eines Systems infolge von Zwangsspannungen J 
ohne die entsprechende Verzerrungsberechnung 
ermittelt werden kénnen; es geniigt, wenn die 
wahren Verzerrungen (also die Matrix @) infolge 
von Verschiebungen ¥ in Richtung von R be- 
kannt sind. 


Es sind vorgeschrieben: 1) die Zwangs- 
spannungen 


J= {Jp} 


2 Es ist hier angenommen, daf die Verschiebungen r in Richtung der Krafte R von Problem I stattfinden. 
3 Es ist hier angenommen, da die Krafte R in Richtung der Verschiebungen r von Problem I wirken. 
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wie in Problem III und 2) die Verschiebungen 
r=0. 
Die Vertraglichkeitsbedingung 7 = 0 wird 
unter Benutzung der Probleme I und IIT 
r=FR+b'H=0, 
somit 
R=—F0'ff. (37a) 
Die Spannungen S und Verzerrungen  ergeben 
sich nun aus (21a) und (34a) zu 
S=— [bF"bd' +b, Db) H,) 
(38a) 
r=[S+H. 


wie in Problem III und 2) Krafte 
R= 0. 
Die Gleichgewichtsbedingung R = 0 wird unter 


| Benutzung der Probleme I und ITT 


R=Kr-+aJ=0, 


| somit 


r=—Ka J. (37b) 

Die Verzerrungen U und Spannungen S ergeben 
sich nun aus (21b) und (34b) zu 

»o=—— [aka +a,C'aJ, ) 


(38b) 
S=kv4+J. 


Im allgemeinen sind in statischen Problemen die Krafte und nicht die Verschiebungen vorgeschrie- 
ben. Somit sind bei Anwendung der Kraft- bzw. Deformationsmethode die Probleme T und IIT 


bzw. II und IV von besonderem Interesse. 


e) Problem V. 


Es ist in einem hochgradig statisch unbe- | 
stimmten Tragwerk oft vorteilhaft, ein Grund- | 


system, das selbst statisch unbestimmt ist, zu 
wahlen. Wir bezeichnen mit Z die Spalten- 
matrix der statisch unbestimmten GréBen im 
Grundsystem. Die vollstandige Kraftmatrix 
ist dann 


(RZX. 


Wir untersuchen nun folgende Verallgemeine- | 


rung von Problem I. 


Es seien gegeben: 1) die Krafte R und 2) die 
unbestimmten GréBen Z fiir jedes X — 1 und 
R=1. Es sollen X und § bestimmt werden. 


Wir setzen die wahren Spannungen S wieder 
in der Form an 


S=b,R+b,X=bR,  (2la) 


wo die Matrizen 6, und 6, aus statischen Be- 
trachtungen allein berechnet werden kénnen, 
da Z fiir jedes R und X bekannt ist. Wir fiihren 
auch die Spannungsmatrix S, die mit R und X 
nur statisch vertraglich sein braucht, ein. Es ist 


S=b,R+b,X, (39a) 


wo b, (bzw. b,) eine Spannungsmatrix ist, deren 
jede Spalte mit dem entsprechenden R = 1 


tisch vertraglich ist. 


Es ist z. B. méglich, 
by (6,) in einem statisch bestimmten Teilsystem 


zu finden, in welchem Falle die Definition von b, 


mit der von b, im Problem I iibereinstimmt. 


von X ergibt unter Benutzung der zwei alter- 
nativen Formen von (10a) 


bv=bj fS=b, fb,R+6,fb,X —0 


Es ist in einem hochgradig kinematisch un- 
bestimmten Tragwerk oft vorteilhaft, ein Grund- 
system, das selbst kinematisch unbestimmt ist, 
zu wahlen. Wir bezeichnen mit W die Spalten- 
matrix der kinematisch unbestimmten GréBen 
im Grundsystem. Die vollstandige Verschie- 


_ bungsmatrix ist dann 


fr WU}. 


Wir untersuchen nun folgende Verallgemeine- 
rung von Problem I. 

Es seien gegeben: 1) die Verschiebungen Yr 
und 2) die unbestimmten GréBen W fir jedes 
U =1lundr =1. Es sollen U und v bestimmt 
werden. 

Wir setzen die wahren Verzerrungen U wieder 
in der Form an 


v=a,r+a,U =ar, (21b) 


wo die Matrizen @ und @, aus kinematischen 
Betrachtungen allein berechnet werden kénnen, 


da W fiir jedes r und U bekannt ist. Wir fihren 


| auch die Verzerrungsmatrix U, die mit r und U 
| nur kinematisch vertraglich sein braucht, ein. 
| Es ist 


v=ar+a,U, (39b) 


wo @) (bzw. @,) eine Verzerrungsmatrix ist, 
deren jede Spalte mit dem entsprechenden r = 1 


(baw) X= 1) und X= 08 (haws E20) seas | OAS ee en Oe 


matisch vertraglich ist. Es ist z. B. méglich, 
a, (@,) in einem kinematisch bestimmten Teil- 
system zu finden, in welchem Falle die Defi- 
nition von @, mit der von @ im Problem I 


| ubereinstimmt. 


Die Vertraglichkeitsbedingung in Richtung 


Die Gleichgewichtsbedingung in Richtung 
von U ergibt unter Benutzung der zwei alter- 
nativen Formen von (10b) 


a,S=—akv=akayr+a ka,U —0 
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oder 
bo=b fS=b,fb,R+b,fb,X—0. 
Somit kann X wieder durch (22a) 
X=—D"R,D 
ausgedriickt werden, wo jetzt 
D=b, fb, =b fb; = b; fb, 
und 
1-0, | bb, fb, = bi fb, 


ist. Die Spannungsmatrix S und die Nach- 
giebigkeit F ermitteln sich nun aus (21a) und 


(22a) 


(40a) 


oder 
a,S=akv=akar+aka,U =0. 
Somit kann U wieder durch (22b) 


U=—C1Gr (22b) 
ausgedriickt werden, wo jetzt 
C=a,ka,=a,ka, =a, ka, 
und (40b) 


C, =a, ka,=a,ka,=a\ka, 


ist. Die Verzerrungsmatrix © und die Steifig- 
keit KC ermitteln sich nun aus (21b) und (26b). 


(26a). Die Einfiithrung der statisch vertraglichen | Die Kinfiihrung der kinematisch vertraglichen 
Matrizen @), @, kann eine besondere Verein- 


fachung der Rechnung bedeuten. 


Matrizen by, 6b, kann eine besondere Verein- 
fachung der Rechnung bedeuten. 


Wir bemerken abschlieBend, daB Probleme II bis IV ahnlich verallgemeinert werden kénnen. 


5. Anwendung der Kraftmethode auf Systeme mit Ausschnitten. Bei der praktischen An- 
_ wendung der Matrizenstatik in Verbindung mit dem elektronischen Rechenautomaten ist eine der 
wichtigsten Aufgaben, die sorgfaltige und wiederholte Nachpriifung der Elemente (des K oeffizienten- 
schemas) der Grundmatrizen by, 6, und f. Um die Kontrolle zu vereinfachen und leicht iibersehbare 
Matrizen zu erhalten, ist es wichtig, alle Spezialfalle, die eine besondere Betrachtung bendtigen, 
zu vermeiden. Kin derartiger Fall tritt z. B. in einem membranformigen System wie einem Fliigel 
auf, wenn einzelne Elemente zwischen den Neétzlinien fehlen (Tragwerk mit Ausschnitten). Es 
erweist sich, da dann im allgemeinen ein komplizierteres Grundsystem (fiir die Berechnung der 
b,-Matrix) als fiir das entsprechende Tragwerk ohne Ausschnitte gewahlt werden mufB. Auch ist 
es notwendig, in der Nahe der Ausschnitte besondere statisch unbestimmte Spannungssysteme X, 
die das sonst regulare Schema der 6,-Matrix stéren und die Kontrolle erschweren, einzufihren. 
SchlieBlich kann es vorkommen, da einige der Diagonalkoeffizienten in den Gleichungen fir die 
unbekannten X von der gleichen GréB®enordnung wie die iibrigen Koeffizienten in den entsprechen- 
den Gleichungen sind, und da damit die genaue Lésung der Gleichungen auf dem Digitalautomaten 
erschwert wird. 

Um diese Schwierigkeiten in Tragwerken mit Ausschnitten zu vermeiden, ist es vorteilhaft, 
einen Kunstgriff, der vom Verfasser in der bereits erwahnten Arbeit! entwickelt worden ist, an- 
zuwenden. Diese Methode ist auch das ideale Verfahren, um die Neuverteilung der Spannungen 
infolge einer nachtraglichen Einfiihrung von Ausschnitten zu finden, ohne die ganze statische 
Berechnung wiederholen zu miissen. 

Das Prinzip des Verfahrens ist einfach. Um eine regulare (kontinuierliche) Konstruktion und 
dementsprechend ein regulares Schema der Elemente der Grundmatrizen zu erhalten, fithren wir 
so viele zusatzliche Elemente wie notwendig in das System ein. Obwohl die Dimensionen (Quer- 
schnitte oder Dicken) dieser neuen Elemente willkiirlich sind, empfiehlt es sich, diese an die der 
umgebenden Teile des Tragwerkes anzupassen. Selbstverstandlich werden das urspriingliche und 
das neue (kontinuierliche) System unter der gleichen Beanspruchung verschiedene Spannung auf- 
weisen. Es ist aber méglich, in den beiden Tragwerken identische Spannungen zu erzielen, wenu wir 
den neuen (fiktiven) Elementen des kontinuierlichen Systems derartige Anfangsverzerrungen ein- 
pragen, daB die Gesamtspannungen (infolge der gegebenen Beanspruchung und der Anfangs- 
verzerrungen) sich in diesen Elementen auf Null reduzieren. Damit sind die neuen Elemente tat- 
sachlich eliminiert, obwohl das regulare Schema der Matrizen und Gleichungen erhalten bleibt. 
Es zeigt sich auch, daB die Matrizenformulierung uns erlaubt, die notwendigen Anfangsverzerrungen 
und, was noch wichtiger ist, die assoziierten Spannungen sehr einfach aus der Spannungsberech- 
nung des kontinuierlichen Systems unter der ausschlieBlichen Wirkung der gegebenen Bean- 
spruchung abzuleiten. Die wahren Spannungen im urspriinglichen System kénnen jetzt sehr ein- 
fach aus der Spannungsuntersuchung des neuen Systems durch Uberlagerung der Spannungen 
infolge der gegebenen Beanspruchung und der Anfangsverzerrungen bestimmt werden. Man 
beachte, da® nur eine statisch unbestimmte Berechnung notwendig ist: die des kontinuierlichen 
Systems unter der gegebenen Beanspruchung. Abb. 9 erlautert das angewandte Prinzip. 


1 Vgl. FuBnote 2 von S. 174. 
14 
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Eine einfache Uberlegung zeigt, daB die zu dieser Methode des Kraftverfahrens duale Methode 
des Deformationsverfahrens die Berechnung eines Tragwerks ergibt, von welchem einige Elemente 
(unendlich) starr sind. Wir werden in diesem Abschnitt nur das praktisch wichtigere Kraftverfahren 
mathematischentwickeln, verweisen aber auf den nachsten Abschnitt fiir die Ergebnisse des dualen 


Verfahrens. 


eft San 


Stab h soll eliminiert werden Anfangsverzerrung H auf Element h allein Tragwerk ohne Element h 


Abb. 9. Berechnung eines Tragwerks mit Ausschnitten. 


Das urspriingliche Tragwerk mége h Ausschnitte haben, die durch Einfiihrung entsprechender 


Elemente ausgefiillt werden. Es sei nun S die Spannungsmatrix im kontinuierlichen System infolge | 
der gegebenen Beanspruchung. Die Gesamtspannungen S, einschlieBlich der Spannungen infolge | 


der Anfangsverzerrungen H[ auf die h Elemente ermitteln sich gema® (21a) und (39a) aus 


S,=S—b,D"b,H, (41) | 


wo b,, die den h Elementen entsprechende Untermatrix von 6, ist. Wir erhalten nun aus der | 


Bedingung, daB die Gesamtspannung in den h Elementen verschwindet, 


S., =S, —6,, D* b,, H aa 


oder 
H = [b,, D7 6,)7S, . (42) 
Somit sind die Spannungen im urspriinglichen System mit den Ausschnitten 
S, =S — b, D5, [b,, D7 6,,)7 S,,. (43) 


Es ist zu beachten, da wir bei dieser Methode die Inversion der Matrix 
[b,, D+ b;, 


benétigen. Die Anzahl r der Zeilen oder Spalten dieser Matrix ist gleich der Anzahl der Spannungen?> 
S (und damit auch der Anfangsverzerrungen H) in den Elementen h. In einem hochgradig statisch 


unbestimmten System wie in einem Flugzeugfliigel ist deswegen zu erwarten, daB r im allgemeinen — 


viel kleiner als die Anzahl der Unbekannten im urspriinglichen System ist. Damit zeigt sich auch, 


daB im Falle einer nachtraglichen Einfithrung von Ausschnitten in eine bereits berechnete Kon- 
struktion das neue Verfahren viel schneller als die direkte Untersuchung des Systems mit Aus- 
schnitten sein kann. 


Wenn das Tragwerk ausschlieBlich von Kraften R beansprucht ist, kénnen wir (43) in der Form 


ansetzen, wo 
b, = b — b, Db, [b,, D> bi, Jb, (44a) 
1st. 

Das Verfahren erlaubt auch eine einfache Ableitung der Nachgiebigkeit F., des urspriinglichen 
Systems, wenn die Nachgiebigkeit W des kontinuierlichen Systems bereits bekannt ist. Tatsiichlich 
ergibt die Anwendung von (25a) und (36a) auf beide Tragwerke (1, sind die Verschiebungen des 
ersteren Systems) 


a=, R=r+b,H=FR + 6, [b,, Db.) 6, 


oder 
Wp) UEDA fh 

Hierin ist “ 

AF = bj (by, DB,,}b, (45a) 

das Inkrement der Nachgiebigkeit infolge der Ausschnitte. 


Streng genommen miissen wir nur diejenigen Spannungen S in den Elementen h gleich Null setzen, die 


linear unabhangig sind. Sonst wird die obige Matrix singular, Die A ilen i i 1 
Anzahl der linear unabhangigen S in den Bienceen h. eee Caen aoa 


Soper (44) 
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Selbstverstandlich bedeutet die Einfiithrung von Elementen in Ausschnitten eine Erhéhung des 
Grades der statischen Unbestimmtheit des Systems; aber dies ist unbedeutend, wenn ein elek- 
tronischer Digitalautomat benutzt wird. Die Methode hat sich in der Praxis als sehr einfach und 
zeitersparend erwiesen. 

Beispiel. Fiir das in Abb. 3a gezeigte Fachwerk soll die Spannungsmatrix b, und die Nach- 
giebigkeit F fiir die Lasten R, und R, bestimmt werden, wenn der Stab 4 fehlt (siehe auch Abb. 8). 
Alle tibrigen Stabe sollen wieder die gleiche Nachgiebigkeit I/F A besitzen. 

Das kontinuierliche Fachwerk mit dem eingesetzten Element h (4) ist in Abschnitt 4a) berech- 
net. Wir entnehmen den Gleichungen (27a) und (27b) 


by, =[0—1//2] und b, = 2, [83 24). (46) 
Unter Benutzung der in 4a) bestimmten Matrix D™ finden wir 
peak: aD 
[On On = ae ; (47) 
AbschlieBend ermitteln wir b, und AF aus (44a) und (46a): 
Ry Ry 
~ — 20.625 0 7 J) 1 861 249 | 
—165.0 35 2 IN je ae (48) 
34.375 0 3 SETS 1p 
0 0 4 
34.375 0 i Die Nachgiebigkeit F’, ergibt sich aus (45) 
1 3 na 
b= Se eee 0 6 unter Benutzung von (30a) fir F. 
20615 -V 2° <0 i 
110.)'2 55) 2 8 
— 34.375-y2 - 0 9 
55) 2 0- =| 10 


Selbstverstandlich ist die direkte Untersuchung des einfach unbestimmten Fachwerks von Abb. 8 
trivial einfach. Die obige Berechnung ist auch nur als eine Erlauterung des Verfahrens zu verstehen. 


Spanning 


H Verzerrung On, v Verzerrung 
Abb, 10a. Berechnung eines modifizierten Systems nach dem Abb. 10b. Berechnung eines modifizierten Systems nach dem 
Kraftverfahren. Deformationsverfahren. 


6. Die Berechnung von modifizierten Systemen. Wir verallgemeinern nun die Methode von 
Abschnitt 5 auf die Berechnung von modifizierten Tragwerken. Unter einem modifizierten Trag- 
werk verstehen wir ein System, das aus einem urspriinglichen System durch Abanderung der Nach- 
giebigkeit oder Steifigkeit einzelner Elemente erhalten wird. Wir behandeln hier folgendes Problem. 
Unter der Annahme, daB die Spannungsverteilung in einem Tragwerk infolge einer gegebenen Be- 
anspruchung schon ermittelt worden ist, sollen die Spannungen bestimmt werden, wenn einzelne 
Elemente nachtraglich modifiziert werden. Die Lisung dieses Problems nach dem Kraftverfahren 
erfolgt ahnlich wie in Abschnitt 5. Demnach werden den zu modifizierenden Elementen derartige 


14* 
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Anfangsverzerrungen H eingepragt, daB die gesamten Verzerrungen infolge der gegebenen Bean- 
spruchung und der HT die gleichen wie im modifizierten System unter der gegebenen Beanspruchung 
allein sind. Diese Bedingung fiihrt wie im Falle der Ausschnitte auf eine sehr einfache Matrizen- 
formulierung der Spannungen im modifizierten System abhangig von der bekannten Spannungs- 
verteilung im urspriinglichen System. Im folgenden wird auch das Deformationsverfahren auf das 
gleiche Problem angewandt, und die Abb. 10a und 10b erlautern das physikalische Prinzip der 
Lésung nach den beiden dualen Verfahren. Die mathematische Entwicklung bestatigt, daB die 


Methode von Abschnitt 5 nur ein Spezialfall der neuen Methode ist. 
Wir nehmen an, da das urspriingliche System g + h Elemente hat, von denen h nachtraglich 
modifiziert werden. Der Index m bezeichnet GréBen im modifizierten System. 


Das Kraftverfahren 


Die Nachgiebigkeit der unverbundenen Ele- | 
mente ist 


0 hh 
im urspriinglichen System und wird | 


| 


| (50a) 


(49a) 


, 0 
-| Om. fey. 


im modifizierten Tragwerk. 


Das Deformationsverfahren 


Die Steifigkeit der unverbundenen Elemente 


ist 
k, 0 
k = (49b) 
0 k, 
im urspriinglichen System und wird 
0 
ke | (50b) 
0 k, + Ak, 


im modifizierten Tragwerk. 


Die statische Untersuchung des urspriinglichen Systems 


Die bekannten Spannungen im urspriing- 
lichen System infolge der gegebenen Bean- 
spruchung seien mit S bezeichnet. Zusatzlich 
wirken noch die Anfangsverzerrungen H auf die 
unverbundenen Elemente h. Die Gesamtspan- 
nungen berechnen sich aus Problem I und IIT 
zu [siehe auch (41) | 


SS -5,D 1b iT (51a) 


Die Gesamtverzerrungen in den Elementen h 
sind demnach 


©, =f, Snr +H=f, [S,—b,,D 16, H] +H. 
(52a) 


Die statische Untersuchung 

Da wir vorschreiben, daB die Spannungen 

und Verzerrungen in den modifizierten Ele- 

menten fh unter der gegebenen Beanspruchung 

~ allein identisch mit (51a) bzw. (52a) sein sollen, 
ergibt sich 


Cnn = Lf, + Afi] (S, — Ox, D+ 0), H] = %, - 
(53a) 
Unter Benutzung von (52a) folgt dann 
wo 
P = [b,, D* by, + Af;'] (55a) 
ist. Man beachte, daB die Anzahl der Zeilen 


der Matrix P gleich der Anzahl der Spannungen 
S in den h Elementen ist (vgl. S. 188). Durch 
Kinsetzen von (54a) in (51a) erhalten wir 


S,,=S—b, D“bi,P“S, , 
ES — is Se : 


| (56a) 


Die bekannten Verzerrungen im urspring- 
lichen System infolge der gegebenen Bean- 
spruchung seien mit © bezeichnet. Zusatzlich 
wirken noch die Zwangsspannungen J auf die 
Elemente h mit festgehaltenen Knotenpunkten. 


Die Gesamtverzerrungen berechnen sich aus 


Problem I und III zu 
v,, = v—a, Ca, J. (51b) 


Die Gesamtspannungen in den Elementen h 
sind demnach 


S, =k, v,,, +3 =k, [v,—a,,C1 a), J] +d. 


(52b) - 


des modifizierten Systems 


Da wir vorschreiben, daB die Verzerrungen 
und Spannungen in den modifizierten Elementen 
h unter der gegebenen Beanspruchung allein 
identisch mit (51b) bzw. (52b) sein sollen, ergibt 
sich 


San = [k, + AR] [v, —a,, Ca, J] =S,. 


(53b) 
Unter Benutzung von (52b) folgt dann 
J= OF VU,» (54b) 
wo 
Q = [a,,C 7 a, + Ake] (55b) 


ist. Man beachte, daB die Anzahl der Zeilen 
der Matrix Q gleich der Anzahl der Verzerrungen 
v in den h Elementen ist. Durch Einsetzen von 


(54b) in (51b) erhalten wir 
v, = 0—4,C74',0%,, 


56b 
Se [ee Uae | Ge) 
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Damit sind die Spannungen und Verzerrungen 
des modifizierten Systems lediglich durch die 
Spannungsverteilung im urspriinglichen System 
ausgedriickt. 


Wenn die Spannungen S nur von Lasten R 
herrithren (S =DBR), wird die modifizierte 
Matrix b,, (S,, = b,, R) 


b,, =b—b, D“*bi,P1b,. (57a) 


Damit sind die Verzerrungen und Spannungen 
des modifizierten Systems lediglich durch die 
Verzerrungsverteilung im urspriinglichen Sy- 
stem ausgedriickt. 


Wenn die Verzerrungen © nur von Ver- 
schiebungen Y herrithren (Uv = @Y), wird die 
modifizierte Matrix a,, (U,, = @,, 1) 


m 


a, = a, Ce ay, O28 ay, Vv (57b) 


Nachgiebigkeit des modifizierten Systems 


Die Anwendung von (25a) und (36a) auf 
das urspriingliche und modifizierte System 
unter der Belastung R ergibt 


FF =F + AF, (58a) 
AF =b,P-b, . (59a) 


AF ist das Inkrement (positiv oder negativ) 
der Nachgiebigkeit infolge der Modifizierungen. 


Die Anwendung von (25b) und (36b) auf das 
urspriingliche und modifizierte System unter 
der Verschiebung 1 ergibt 


K_=K 2 AK, 
AK =a), 0", . 


AK ist das Inkrement (positiv oder negativ) 
der Steifigkeit infolge der Modifizierungen. 


(58b) 
(59b) 


Spezialfalle 


1) Elimination der h Elemente (Ausschnitte) : 
Ks ist 


Vp P =b,,DDi,. (60a) 


Dieser Fall ist in Abschnitt 5 behandelt und 
entspricht Punkt A in Abb. 9a. 


und 


2) Verstarrung der h Elemente: Es ist 
Af, =—f, P=b,, D* b,,—f; . 
(61a) 

Dieser Fall entspricht Punkt R in Abb. 9a. 


und 


Beispiel. Fir das in Abb. 3a gezeigte Fach- 
werk soll die Spannungsmatrix b,, und die Nach- 
giebigkeit Ff, fiir die Lasten R, und R, bestimmt 
werden. Die Nachgiebigkeit des Stabes 10 ist 
1/6 EA, aber die aller tibrigen Stabe ist weiterhin 
l/E A. Wir betrachten das in 4a) berechnete 
Fachwerk mit gleichen Nachgiebigkeiten fiir die 
Stabe als das urspriingliche System, dessen 
Stab 10 nachtraglich modifiziert ist. Wir er- 
halten aus (27a) und (29a) 


by, = [0 1], 
1 = —, 
b, ==. [— 28/2 — 24/2]. (62a) 
I l Swell 
eee Ghd Ob a6 iA 


Unter Benutzung der in 4 a) bestimmten Matrix 
D*™ finden wir 


=e ae 


50 RA (63a) 


1) Verstarrung der h Elemente: 
Es ist 


Dieser Fall ist dual zu der Methode des Ab- 
schnittes 5 und entspricht Punkt R in Abb. 9b. 


2) Elimination der h Elemente: Es ist 
Ak, = —k, Q =a,,C 1 a,,—k;. 
(61b) 


Dieser Fall entspricht Punkt A in Abb. 9b. 


und 


Beispiel. Fiir das in Abb. 3b gezeigte Fach- 
werk soll die Verzerrungsmatrix @,, und die 
Steifigkeit K,, fiir die Verschiebungen r, und r, 
bestimmt werden. Die Steifigkeit des Stabes 10 
ist 6 EA/l, aber die aller iibrigen Stabe ist 
weiterhin F'A/l. Wir betrachten das in 4 a) be- 
rechnete Fachwerk mit gleichen Steifigkeiten 
fiir die Stabe als das urspriingliche System, 
dessen Stab 10 nachtraglich modifiziert ist. Wir 
erhalten aus (27b) und (29b) 


a,,=[0 000 1//2 0 —1//2], 


1 a — 

@, = 555 [12/2 — 1042]. (62b) 
6EA SHA EA 
ee erate 


Unter Benutzung der in 4b) bestimmten Matrix 
C— finden wir 


339 EA 


Qi= 556 fo (63b) 
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und damit, und damit, 
Ry Ry Wil To 
32:20 COME ame [ —44.47 50.40) aed 
73.04 23.32 2 — 62.85 — 13.64 2 
22.88 — 3.96 3 49.81 — 96.65 3 
91.96 31.68 4 72.93 37.95 4 
; 22 6on ee 1300 5 ; 49.81 — 96.65 5 
Oe 4.84 27.72 6 @,, = 335| —05.81 235.39 6 
32.12)/2 3.9672 | 7 44.4772 — 50.40/2 | 7 
18.04/72 23.32/2 | 8 —31.42//2 160.68/2 | 8 
—22.88//2 3.96/2 | 9 —49.81/2 96.6572 | 9 
_ — 36.9672 —31.68/2 | 10 3.56/72  — 30.83/2 | 10 
AF I 31.4 26.9 | mu 1.3 pas a 
Eek CRE ket (64a) AK == (64 b) 
oe 26.0 sh | SEU 2A 95.8 
Wir ermitteln Ff’, aus (58a). Wir ermitteln K,, aus (58b). 


Wir bestatigen wieder, daB b;, a,, = E, erfiillt ist. 


Der Einflu8 der erhéhten Steifigkeit des Stabes 10 auf die Spannungsverteilung ist, wie zu erwarten, 
verhaltnismaBig gering. 


(Eingegangen am 4, August 1956) 


Anschrift des Verfassers: Professor J. H. Argyris D. E., London S. W.7, Prince Consort Road, Imperial College 
of Science and Technology (University of London). 
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Diagonalwinkel-Zuordnungen im Gelenkviereck 
Von K. Hain 


1. Einleitung. Winkelzuordnungen im Gelenkviereck und in mehrgliedrigen Getrieben sind 
ein wichtiges Aufgabengebiet der Getriebesynthese. Es handelt sich dabei im allgemeinen darum, 
die Veranderungen eines Winkels im Gelenkviereck den vorgeschriebenen Anderungen eines anderen 
Winkels zuzuordnen. Auf diese Weise ist es méglich, gegebene Funktionen mechanisch mit Hilfe 
eines Getriebes darzustellen. AuSerdem bilden Winkelzuordnungen eine wichtige Grundlage fiir 
viele andere getriebetechnische Aufgaben. 

Bisher hat man sich jedoch darauf beschrankt, die Getriebesynthese beziiglich der Winkel- 
zuordnungen insofern auszubauen, als man Winkel im Gelenkviereck bzw. Winkel in vielgliedrigen 
Getrieben einander zugeordnet hat, die zwischen drei benachbarten Getriebegliedern lagen. Auf 
dieser Grundlage sind viele Untersuchungen entstanden. 

Neben den bekannten Verfahren fiir Winkelzuordnungen stellt die Zuordnung von Diagonal- 
winkeln im Gelenkviereck eine neuartige Aufgabe dar, da es sich hier um zugeordnete Winkel 
handelt, die von paarweise verschiedenen Getriebegliedern eingeschlossen werden. Einen um- 
fassenden Uberblick iiber die Synthese des Gelenkvierecks hat R. Kraus gegeben!. Auf die Zuord- 
nung von Diagonalwinkeln ist er jedoch dabei nicht eingegangen. Daf eine solche Synthese auch 
einen praktischen Wert hat, kann an Hand von verschiedenen Nutzanwendungen gezeigt werden. 
So hat z. B. G. Kiper an einer Schneidvorrichtung beschrieben, wie durch eine weitgehende Aus- 
nutzung des Gelenkvierecks dieses in der Lage ist, ein vielgliedriges Getriebe vollwertig zu ersetzen?. 
Weitere Beispiele kénnen u. a. aus dem Hebezeugbau angefiihrt werden, wenn zwischen bewegten 
Getriebegliedern ein 6lhydraulischer Zylinder mit Kolben zur Bewegungseinleitung in das gesamte 
Getriebe dient und ein Drehmoment um einen festen Gestellpunkt zu iiberwinden ist. 


2. Ubersetzungsverhiltnis bei Diagonalwinkel-Drehungen. Das Gelenkviereck von Abb. 1 
besteht aus Gestell J und den drei bewegten Hebeln 2, 3 und 4. Laut Aufgabenstellung sollen die 
Veranderungen des Winkels y den Veranderungen des ihm a 
diagonal gegeniiberliegenden Winkels 6 zugeordnet werden. Daz R 
Samtliche Abmessungen des Gelenkvierecks sind dieser Forde- 
rung entsprechend auszubilden. Bei Anspriichen, die die Még- 
lichkeiten des Gelenkvierecks nicht voll ausschépfen, kann 
man sich bestimmte MaBe vorgeben, wie im folgenden gezeigt 
werden soll. 

Die Gelenkpunkte 12, 23, 34, 14 des Gelenkvierecks sind 
durch die Doppelziffern der in einem solchen Gelenk vereinig- 


ten benachbarten Getriebeglieder gekennzeichnet. Zwischen 2 ike 7h 
den Relativdrehungen der Getriebeglieder 1 und 4. um das RG eo ee ae ee ee 
Gelenk 14 und den Relativdrehungen der Getriebeglieder 2 verlaufenden Diagonaldrehungen D,, und D,,. 


und 3 um das Gelenk 23 1a4Bt sich das Ubersetzungsverhaltnis 

durch Strecken ausdriicken®. Bezeichnet man mit D,, die Drehungen des Gliedes 4 relativ zum 
Glied J und mit D,, die Drehungen des Gliedes 3 relativ zum Glied 2, so sind damit auch die Dreh- 
richtungen festgelegt. Schreibt man 


Dy 
Dyp 
so ergibt sich durch die vertikale Zuordnung der Fufzeichenziffern eine sogenannte K-Achse 12—34. 
Diese Achse wird in Abb. 1 von der Achse 14—23 im Punkt Q geschnitten, wobei die Ziffern der 
letzteren Achse die Gelenke bezeichnen, um die die Drehungen D stattfinden. Liegt der Punkt Q 


zwischen den durch die Drehungen D gekennzeichneten Gelenken 14 und 23, so miissen die Dre- 
hungen D,, und D,, entgegengesetzten Drehsinn haben. Die Drehung D,, erfolgt z. B. im Uhr- 


12—34 


1 R. Kraus, Mitt. Techn. Inst. Tung-Chi-Universitat Woosung, Bd. 2 (1935). 

2 G. Kiper, Die Bedeutung der Relativbewegungen fiir das Konstruieren periodischer Getriebe, VDI- 
Bericht Bd. 5, Diisseldorf 1955 (Vortrage der VDI-Tagung Stuttgart 1953, 5. 115/117). 

3 K. Hain, Landtechn. Forschung 3 (1953), S. 97. 
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zeigersinne, die Drehung D3., also die Bewegung des Gliedes 3 gegeniiber dem Glied 2, mu dann 
im Gegenuhrzeigersinne stattfinden'. Man kann auch eine andere Anordnung der FuBzeichen- 
ziffern fiir die Drehungen D wahlen und wiirde dann eine entsprechend andere vertikale Zuordnung 
dieser Ziffern erhalten. Im folgenden sollen jedoch immer mit dem Winkel y die Drehung Dy 
und mit dem Winkel § die Drehung D,, bezeichnet werden. 

Das Gegenstiick zu Abb. 1 zeigt Abb. 2, insofern als der Punkt Q auBerhalb der Gelenke 14 
und 23 liegt. Mit dieser Lage ist aber auch eine gleichsinnige Drehung des Gliedes 4 um 1] und des 
Gliedes 3 um 2 verbunden. In Abb. 2 verlaufen beide Drehungen D,, und D3, um die Gelenk- 
punkte 14 und 23 im Uhrzeigersinne. 


— aa 
Abb. 2. Gelenkviereck mit gleichsinnig verlaufenden Diagonaldrehungen Abb. 3. Getriebetotlagen bei Zuordnungen von Diagonal- 
D,, und D3;,. drehungen. Grenzlagenstellungen zwischen den Gliedern I und 4. 


3. Totlagenstellungen bei Zuordnungen von Diagonaldrehungen im Gelenkviereck. Nach den 
vorangegangenen Ausfihrungen soll die Bewegungseinleitung entweder im Gelenk 23 oder im Ge- 
lenk 14 und die Ableitung der Bewegung im jeweils anderen Gelenk stattfinden. Es ist bekannt, 
daB es Getriebestellungen gibt, bei denen eine Sper- 
rung stattfindet, also eine Bewegungseinleitung Bae —— 


an einer bestimmten Stelle nicht mehr méglich e Sa ‘he 
ist. Solche Sperrstellungen gibt es auch bei der \ ce k sia 
Zuordnung von Diagonalwinkel-Drehungen. In ’ =< 


Abb. 4. Getriebetotlagen bei Zuordnungen von Diagonaldrehungen. Abb. 5. Zuordnungen der Diagonalwinkel y,, und f,, 
Grenzlagenstellungen zwischen den Gliedern 2 und 3. fiir zwei Getriebelagen, z 


Abb. 3 bilden die Glieder 2 und 3 eine Strecklage, so dafB das Glied 4 gegeniiber dem Glied 1 
seine Bewegungen umkehrt. Damit ist es aber auch nicht méglich, am Glied 4 eine Bewegung 
einzuleiten. Eine zweite solche Stellung ist durch die Lage I (2) (3) (4) gekennzeichnet. Auch 
hier liegen die Glieder (3) und (2) auf einer Geraden. Durch die Lagen (4) und 4 ist fiir das gegebene 
Gelenkviereck auch der gesamte Schwingbereich des Gliedes 4 gekennzeichnet. 


_ 7 Es sei darauf hingewiesen, da die Ziffernfolge bei der Bezeichnung der Gelenke ohne Bedeutung ist 
Sie wurde hier so gewahlt, da die kleinere Ziffer immer an erster Stelle steht. Bei den Drehungen D jedoch 
wurde der Ziffernfolge insofern eine Bedeutung gegeben als D,, die Bewegungen des Gliedes 4 gegeniiber dem 
Ghed J und D,, diejenigen des Gliedes 3 gegentiber dem Glied 2 kennzeichnen. 

Die Drehung D,, wiirde also entgegengesetzt zu D,, und die Drehung D,, entgegengesetzt zu Dy sein. 


| 
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In Abb. 4 sind die beiden Getriebestellungen des Gelenkvierecks eingezeichnet, in denen das 
Glied 4 in den Getriebelagen 4 und (4) mit dem Gestell 1 zur Deckung kommt. Dadurch befinden 
sich die Glieder 2 und 3 baw. (2) und (3) in je einer Grenzlage, und in beiden Lagen findet eine Be- 
wegungsumkehr der Relativbewegung beider Glieder statt. Eine Bewegungseinleitung durch eine 
Drehung um das Gelenk 23 ist in beiden Getriebestellungen unméglich. 


4. Diagonalwinkel-Zuordnungen bei gegebenen Gestellgelenkpunkten. Wenn nach Abb. 5 
_ ein Gelenkviereck sich von seiner Lage 4,4,B,B, in eine zweite Lage A,A,B,B, bewegt hat, hat 
das Glied B)B (in Abb. 1 bis 4 Glied 4) gegeniiber dem Gestell A,B, (frither Glied 1) den Winkel y,, 
durchlaufen. Das Glied AB (frither Glied 3) schlie8t mit dem Glied A,A (frither Glied 2) in der 


5 es ne 
Abb. 6. Konstruktion des Kreises k,, als geometrischer Ort fiir alle Punkte Abb. 7. Brauchbarer Bereich des Kreises k,, 
A, bei Zuordnung der Diagonalwinkel y,, und f,,. Gegensinnige Dreh- von Abb. 6. 


richtung von 7,,. und f,). 


Lage J den Winkel (, und in der Lage 2 den Winkel f, ein. Die Veranderung dieses Winkels ist 


damit 
Pre = Bo oa By ° 
Den Winkel /,, findet man sehr leicht, wenn man beispielsweise in der Getriebestellung 1] um 


den Punkt A, den Kreis k, durch B, zeichnet und diesen zum Schnitt B,, mit dem Kreis k, bringt, 
den man um A, durch B, zeichnet. Dann ist 

By, = S Bi A,B, . 
Diese Konstruktion des Winkels /,, liegt allen folgenden Verfahren zugrunde. 

Ist in einem Gelenkviereck der Winkel 6,, dem Winkel y,, zuzuordnen, wobei auch deren Rich- 
tungen vorgegeben sein sollen, so kann man, da durch diese Forderung die Méglichkeiten des 
Gelenkvierecks noch nicht voll ausgeschépft sind, noch verschiedene Annahmen treffen. So gibt 
man sich nach Abb. 6 beispielsweise die Lage der beiden Gestellpunkte By) und A, sowie die Lange 
des Gliedes B,B und dessen Anfangslage B,B, relativ zum Gestell A) By vor. Mit dem gegebenen 
Winkel y,, kann man nun die Lagen B,B, und B,B, des Gliedes B,B zeichnen. Uber der Strecke 
B, B, zeichnet man das gleichschenklige Dreieck mit dem Spitzenwinkel f,, und den Spitzenpunkt D, 
wobei die gegebene Richtung von /3,, insofern erfiillt sein muB, als der Schenkel DB, in der durch f,, 
gekennzeichneten Richtung nach DB, drehen muB. Weiterhin zeichnet man tiber der Strecke B, A) 
ebenfalls ein gleichschenkliges Dreieck mit dem Spitzenwinkel /,, und dem Spitzenpunkt M. Auch 

‘hier muB die gegebene Richtung von /),, insofern erfiillt sein, als der Strahl MB, nach MA, wandert. 
Zeichnet man nunmehr um den Punkt M den Kreis k,, durch den Punkt D, so ist dieser Kreis der 
geometrische Ort aller Punkte A, fiir die gegebenen Winkel ,, und f,, sowie fiir die angenommenen 
weiteren Getriebedaten. Es erfiillt also z. B. sowohl das Gelenkviereck A)A,B, By als, auch jedes 
andere beliebige Gelenkviereck 4,4, B' B, mit dem Punkt A, baw. A; auf dem Kreis k,, die gestellten 

Bedingungen. 
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Entsprechend Abb. 6 ist der Kreis k,, der geometrische Ort fiir alle Gelenkpunktlagen Aj. 
Stellt man die Forderungen, daB ein stetiger Verlauf der Winkel y und 6 wahrend des vorgegebenen 
Bewegungshereichs vorhanden ist, daB also zwischen den vorgegebenen Lagen keine Bewegungs- 
umkehr stattfindet, so ist hierfiir nur ein bestimmter Bereich des Kreises k,, als brauchbar zu 
bezeichnen. Bei gegebenem Richtungssinn der Winkel y,, und /,, muB ein Punkt Q entsprechend 
Abb. 1 konstruierbar sein, so daB, wie in Abb. 7 gezeigt, nur der Kreisbogen des Kreises ky. zwischen 


Abb. 8, Konstruktion des Kreises k,, als geometrischer Ort fiir alle Punkte A, bei Zuordnung der Diagonalwinkel y,, und fj». 
Gleichsinnige Drehrichtung von y,, und f,,. 


Abb. 9. Brauchbarer Bereich des Kreises k,, von Abb. 8. 


den Punkten A, und (A,) brauchbar ist, wenn die hierdurch gekennzeichneten Getriebelagen den | 
Totlagen der Abb. 3 entsprechen. Selbstverstandlich wiirden auch Punkte A, auBerhalb dieses | 


Bereiches die gestellte Bedingung erfillen, d. h. es wiirde sich das Glied B,B um den Winkel yp 
drehen; dazwischen kénnte aber z. B. eine Lage dieses Gliedes eintreten, die auRerhalb des dave 


den Kreisbogen B,B, gekennzeichneten Bereiches liegt. 


Die Abb. 8 und 9 entsprechen den Darstellungen der Abb. 6 und 7, wobei hier eine gleich- 
sinnige Richtung, namlich eine Drehung in Uhrzeigerrichtung der Winkel Yi2 und f,, angenommen 
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wurde. Nach Abb. 8 liegt dann der Punkt D als Spitzenpunkt des gleichschenkligen Dreiecks 


B,DB, auf der konkaven Seite des Kreisbogens B,B,, wahrend der Punkt M auf der anderen 
Seite der Strecke B,A, wie in Abb. 6 liegen muB. 


In Abb. 6 und 8 waren zu Beginn der Konstruktion noch verschiedene willkiirliche Annahmen 
gegeben worden, die man insofern als Konstruktionsfreiheiten ausnutzen kann, als mehrere Winkel 
y der gleichen Anzahl Winkel 8 zugeordnet werden kénnen. In Abb. 10 sind die Winkel Wy. und 
13 vorgeschrieben, denen die gegebenen Winkel ,, und §,, zuzuordnen sind. Fiir diese Forderungen 
kann man nach wie vor die Lange des Gliedes B,B und seine Anfangslage B)B, zum Gestell A,B, 
beliebig annehmen. Damit ergeben sich die Lagen B,B,, B,B, und B,B, dieses Gliedes mit den 
eingeschlossenen gegebenenen Winkeln y,, und y,;. Uber der Strecke B,B, zeichnet man nach 
den Regeln der Abb. 6 das gleichschenklige Dreieck mit dem Spitzenwinkel fy. und dem Spitzen- 
punkt D,,, in der gleichen Weise iiber B,B, das gleichschenklige Dreieck mit dem gegebenen 
Spitzenwinkel 6,, und dem Spitzenpunkt D,,. Verbindet man B, mit Ay, so miissen auch iiber dieser 


WA 10 
is SR 
Lie if Sie ca ON 


Abb. 10. Konstruktion eines Gelenkvierecks fiir eine zweifache Winkelzuordnung y;,, y,3 2U fy2, 6,3. Gestellpunkt A wandert 
auf Strahl B,X,. 


Strecke die beiden gleichschenkligen Dreiecke 4)M,,.B, und A)M,,B, gezeichnet werden. Fiir die 
Winkel y,, und f,, ist nunmehr der geometrische Ort aller Punkte 4, der Kreis k,, um M,, durch 
D,,. Fiir die Winkel y,; und f,3 ist aber auch der Kreis k,, der geometrische Ort aller Punkte A,, 
wenn dieser Kreis um M,, durch D,, gezeichnet wird. Die Kreise k,, und kj, schneiden sich in 
zwei Punkten A, und Aj(.), so dafs also fiir den angenommenen Gestellpunkt A, zwei Gelenkvier- 
ecke mit den Punkten A, bzw. Aj) die gestellten Bedingungen erfiillen. 

Es gibt nun fir eine zweifache Winkelzuordnung y,,. und y,3; zu /,. und /,; unendlich viele 
Gelenkvierecke, wenn man die Lage des Gestellpunktes A) verandert. Lat man diesen Punkt 
z. B. auf dem beliebig angenommenen Strahl] B,X, wandern, indem man andere Punkte A) an- 
nimmt, so kann man die weiteren Strahlen B,X,, und B,X,, durch M,, und M,, und die Parallelen 
A\M), 21 A,M,;, sowie A\M), zu A,M,, zeichnen. Dann schneiden sich die Kreise k,, und 
kj, um M), durch D,, und um Mj, durch D,; in neuen Punkten A; und Aj). Alle diese A,- 
Punkte bilden eine Kurve 6, die den geometrischen Ort der Punkte A, fiir eine geforderte zweifache 
Winkelzuordnung 71, #13 2U Py9, Pig darstellt. 

Wenn es fiir alle Lagen des Punktes A, auf dem Strahl B,X,) unendlich viele Gelenkvierecke 
mit vorgeschriebener zweifacher Winkelzuordnung gibt, so muf es fiir die Lage des Punktes A, 
auf diesem Strahl B,X, eine bzw. mehrere Lésungen geben, wenn man noch ein drittes Winkelpaar 
Wy4> 644 hinzunimmt. Man kann dann z. B. eine zweite Kurve 6 fiir die Winkelzuordnungen 79, 
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P14 2U Byo5 Bq zeichnen. Die Schnittpunkte dieser zweiten d-Kurve mit der ersten stellen Lésungen 
fiir die dreifache Winkelzuordnung dar. Schwenkt man auferdem noch den Strahl B,X, um den 
Punkt B,, so hat man damit eine weitere Freiheit gewonnen, und alle Schnittpunkte der ent- 
sprechenden 5-Kurven ergeben dann eine Kurve als geometrischen Ort aller Punkte A, fiir drei- 
fache Winkelzuordnungen. Diese neuen Kurven kénnen sich nun wiederum gegenseitig schneiden, 
wenn man noch ein Winkelpaar y,;, 8,; hinzunimmt. Damit sind aber die Méglichkeiten des Ge- 
lenkvierecks fiir vorgegebene Diagonalwinkel voll erschépft. Die gleichen quantitativen Forderun- 
gen kann man auch bei den bisher behandelten Winkelzuordnungen erfillen, was durch die so- 
genannte Wertigkeitsbilanz’ nachgewiesen werden kann. 


Abb. 11. Konstruktion eines Gelenkvierecks fiir eine zweifache Winkelzuordnung yj», y13 ZU B42, By3- Gestellpunkt A, 
wandert auf Kreisbogen r, um B,. 


Fiir eine zweifache Winkelzuordnung kann man in ahnlicher Weise wie in Abb. 10 auch noch 
unendlich viele Lagen des Gestellpunktes A, finden, wenn man diesen Punkt auf einem Kreisbogen 


ro um B, wandern laBt (Abb. 11). Dann miissen auch die Werte M,, und M,, auf je einem Kreis- _ 


bogen r,, und r,, um B, wandern. Das bedeutet aber, daf man das gesamte Dreieck 4)M,,.M,, 
um den Punkt B, verschwenkt. Dann braucht man nur noch um die neuen Punkte M/, und Mj, 
die Kreisbogen k', durch D,, und k\, durch D,,; miteinander zum Schnitt zu bringen und erhalt 
neue Punkte A; und Aj). Alle diese Punkte ergeben eine Kurve ¢. Diese Kurve ¢ stellt also, 
ahnlich wie in Abb. 10 die Kurve 6, den geometrischen Ort aller Gelenkpunktlagen A, fiir die 


zweifache Winkelzuordnung dar. Bei der Forderung nach dreifacher Winkelzuordnung kann man — 


die Halbmesser der Kreise r verandern und erhalt jedesmal eine neue Kurve e, die sich mit der 
vorhergehenden in einem oder mehreren Punkten schneidet. Alle diese Schnittpunkte ergeben 


dann neve Kurven als geometrische Orter fiir alle Gelenkpunktlagen A, bei dreifacher Winkel-_ 
zuordnung. Diese Kurven miissen mit den entsprechenden, nach Abb. 10 ermittelten iiberein- © 


stimmen. Sie entsprechen den im Schrifttum vielfach erwahnten Mittelpunktkurven. 


5. Diagonalwinkel-Zuordnungen bei gegebenen beweglichen Gelenkpunkten. Bisher wurde | 
auBer den Punkten B, und B die Lage des Gestellpunktes A) angenommen und fir diese Annahmen | 


der bewegliche Gelenkpunkt A bestimmt. Geht man den umgekehrten Weg, indem man die B- 
Punkte und den beweglichen Gelenkpunkt A annimmt, so ergeben sich fiir die Lage des Gestell- 
punktes A, noch einfachere Verfahren als bisher. In Abb. 12 ist wiederum die Winkelzuordnung 


Y12» Big Vorgeschrieben, angenommen werden die Lagen B)B, und B)B,, die den Winkel y,, ein- | 
schlieBen, und die Lage des Gelenkpunktes 4,. Dann braucht man nur den Punkt B, um A, und | 


um den Winkel f,) bis B,, zu verschwenken und die Mittelsenkrechte g auf der Strecke B,B,, 


zu zeichnen. Diese Mittelsenkrechte ist der geometrische Ort aller die gestellte Aufgabe erfiillender | 


Punkte Ay. Jedes beliebige Gelenkviereck, z. B. auch A)’ A,B,By garantiert somit die geforderte 
Winkelzuordnung. 


1 R. Kraus, Feinwerk-Technik 56 (1932) S. 57. 
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_Mit Hilfe der Konstruktion von Abb. 12 kann man miihelos auch eine Zweiwinkel-Zuordnung 
beherrschen. Sind also in Abb. 13 die Winkel y,. und y,; den gegebenen Winkeln §,, und {3 zu- 
zuordnen, so kann man immer noch einen beliebigen Punkt A, annehmen und fiir diesen die Kon- 
struktion nach Abb. 12 sowohl fiir den Winkel f,, als auch f,, durchfiihren. Damit erhalt man die 
Punkte B,, und B,, und die Mittelsenkrechten g,. und g,, der Strecken B,B,, und B,B,,. Der 
Schnittpunkt A, der beiden Mittelsenkrechten g,, und g,, erfiillt somit die gestellten Bedingungen. 
Nimmt man auf der Geraden B,A, weitere beliebige Punkte (A,) an, so braucht man nur noch 
um (A,) den Kreisbogen durch B, zu zeichnen, der von den Parallelen durch (A,) zu B,B,, und 
A,B;, in den Punkten (B,,) und (B,,) geschnitten wird. Damit ergeben sich neue Mittelsenk- 
rechten (g;,) und (g,3), deren Schnittpunkt (4y) wiederum die gestellte Aufgabe erfiillt. Alle so 
ermittelten Punkte A, liegen auf einer Kurve 6, und diese Kurve ist der geometrische Ort aller 
Gestellpunkte A), wenn der Gestellpunkt 4, auf einer vorgegebenen Geraden wandert. 

Man kann auch hier, ahnlich wie in Abb. 11, den 
By Punkt A, auf einem Kreis um B, wandern lassen und 
~~ erhalt dann eine Kurve ¢ wiederum als geometrischen 
| Ort fiir alle Punkte A). 


WN 
} “ ae a 3.1, —— 


(92) | | (I) 


Abb. 12. Konstruktion der Mittelsenkrechten g als geo- Abb. 13. Konstruktion eines Gelenkvierecks fiir eine zweifache Winkelzuordnung 
_metrischer Ort fiir alle Gestellpunkte A, bei Zuordnung W909 1g 24 By, Pyg. Gestellpunkt A, wandert auf beliebig angenommenem 
der Diagonalwinkel y,. und 4». Strahl B,A,. 


Fir jede beliebig angenommene Strecke B,A, nach Abb. 13 gibt es also unendlich viele Punkte 
A, auf der Kurve 0, und bei jeder anderen Lage der Strecke B,A, gibt es eine neue Kurve 6, die 
sich mit der alten in einem oder mehreren Punkten schneidet, wenn man ein drittes Winkelpaar 
Wr4 Py, vorschreibt. Alle diese Schnittpunkte bilden wiederum eine Kurve, die hier der im Schrift- 
tum vielfach erwahnten Kreispunktkurve fiir die wblichen Winkelzuordnungen entspricht. Fir 
ein viertes Winkelpaar y,;, /;, kann man die Kreispunktkurven zum Schnitt bringen und erhalt 
dann nur noch eine begrenzte, der Wertigkeitsbilanz entsprechende Anzahl von Lésungen. 


6. Punktlagen-Reduktionen. In gleicher Weise wie bei den bisher itiblichen Winkelzuord- 
-nungen! kann man auch hier Punktlagenreduktionen verwenden, wenn man dafiir sorgt, daB zwei 
—Gelenkpunktlagen des Punktes A, in einem Punkt zusammenfallen. In Abb. 14 fallt beispielsweise 
der Punkt A, mit dem Punkt A, zusammen, wenn man tiber dem Dreieck B,B)B, das gleich- 

schenklige Dreieck B,A,B, mit dem eingeschlossenen gegebenen Winkel /,, zeichnet, wobei selbst- 
verstandlich der vorgeschriebene Richtungssinn beider Winkel zu beriicksichtigen ist. Fiir diese 
_Annahme kann man einen beliebig in der Zeichenebene liegenden Gestellpunkt A) wahlen. Das sich 
_daraus ergebende Gelenkviereck 4,A,B,B) wird immer, wie leicht einzusehen ist, das zugeordnete 
~Winkelpaar ,3, 6,3 garantieren. Das Getriebeglied A,A befindet sich dabei in einer so gekenn- 
-zeichneten Getriebestellung in einer Doppellage, und zwischen den, dieser Doppellage entsprechen- 
den beiden Getriebestellungen mu es eine Grenzlage A)A, erreicht haben. Diese Grenzlage ist 


1 K, Hain, Zeitschr. Instrumentenkunde 63 (1943) S. 170. 
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dadurch gekennzeichnet, daB die Glieder B)B und AB eine Strecklage einnehmen. Da bei der hier 
zur Frage stehenden Aufgabe die Bewegungen des Getriebegliedes 4A nicht vorgeschrieben sind, 
sind diese auch fiir die praktische Ausfithrung gleichgiiltig. Man kann es u. U. sogar als Vorteil 
bezeichnen, da® das Getriebsglied A,A bei der stetigen Bewegung der Glieder B)B und AB hin- 

und herschwingt, denn dadurch braucht man fir 


ee seine Bewegung entsprechend geringeren Platz vor- 
Bee \ zusehen. 
A, p Ams ee Da bei der Vorgabe einer Punktlagenreduktion 


der Punkt <A, beliebig in der Zeichenebene liegen 
kann, braucht man bei der Konstruktion eines Ge- 
lenkvierecks auch auf ihn bei den entsprechenden 
zwei Lagen keine Riicksicht zu nehmen. Ist also in 
Abb. 14 eine dreifache Winkelzuordnung Yj5, Yj, 
Wis ZU Byo, By3, yg Vorgeschrieben und hat man fir 
das Winkelpaar y,3, /,, die Punktlagenreduktion 
A, = A, vorgesehen, so braucht man nur noch an 
A,B, die Winkel f,, und /,, im gegebenen Richtungs- 
sinne anzutragen und erhalt die Punkte B,, und B,,. 
Die Mittelsenkrechten g,. und g,, auf den Strecken 
Abb. 14, Konstruktion eines Gelenkvierecks fiir dreifache Baba aes ByByy OAS SG ue Cestelpaaea | 
Winkelzuordnung mit Hilfe einer Punktlagenreduktion. Ay. In Abb. 14 ist ersichtlich, da8 man nicht nur | 
diese eine Punktlagenreduktion A, = A,, sondern | 


insgesamt sechs derartige Reduktionsméglichkeiten hat, namlich 
Das bedeutet aber, daB es unter Benutzung von Punktlagenreduktionen fiir eine geforderte Drei- | 


winkelzuordnung insgesamt sechs, die Aufgabe erfiillende Gelenkvierecke gibt, von denen man | 
sich das fiir den gedachten Zweck giinstigste heraussuchen kann. 


Institut fiir Landtechnische Grundlagenforschung, Braunschweig-Volkenrode (Leitung: Prof. Dr. W. Kloth). — 
(Eingegangen am 31. August 1956.) 


Anschrift des Verfassers: Ingenieur K. Hain, Braunschweig, Bundesallee 50. 
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Beitrag zur Berechnung in Querrichtung durchlaufender Plattenstreifen 
mit Hilfe Fourierscher Integrale 


Von F. W. Mader 


1. Einleitung. Plattenstreifen mit Navierschen Randbedingungen wurden fiir die verschieden- 
artigsten Lastangriffe berechnet. Eingespannte Plattenstreifen wurden von Haeussler u. a. unter- 
sucht. Im folgenden soll ein Verfahren vorgefiihrt werden, mit dessen Hilfe Plattenstreifen unter dem 
Angriff von Randmomenten beliebiger Art berechnet werden kénnen. Durch die Beriicksichtigung 
der Stetigkeit der Biegeflache iiber der Auflagerlinie eines durchlaufenden Plattenstreifens kénnen 
dann die Stiitzmomente iiber dieser Auflagerlinie gewonnen werden. 


2. Momentenangriff am Plattenstreifenrand. Das lings eines Streifenrandes angreifende 
Moment wird in einen symmetrischen und einen antimetrischen Anteil zerlegt. Da der Rechnungs- 
gang fiir beide Anteile grundsatzlich derselbe ist, wird im folgenden nur der symmetrische Anteil 
weiterverfolgt. Fiir den antimetrischen ist lediglich cos xx durch sin «x zu ersetzen. Das Moment 
wird nun in ein Fouriersches Integral entwickelt, das folgende 
Form aufweist: y 

co 


M, == [ mo) cos ax dx. 
0 
Die Biegeflache mu8 der homogenen Plattengleichung 


SS 
AAw = 0 geniigen und wegen der je zwei Randbedingungen 
an beiden Streifenrandern (Abb. 1) vier Freiwerte enthalten. 
Man wahlt hier zweckmaBigerweise einen Ansatz aus hyper- _ 
bolischen Funktionen von der Form oy, M, 
lee) 
a if (A Gof ay tay B Gin ay hy C Ginay Abb. 1. Randmoment am Rand y = 0. 
0 
tay D€oj ay) cosax dx. 
Er befriedigt die Plattengleichung. Die Randbedingungen lauten fiir den Rand y = 0 
w=0, (1) 
x M 
reas a. 2) 
wobei Punkte Ableitungen nach y, Striche Ableitungen nach x und 
A Et 
= 12(1—) 
die Steifigkeit der Platte bedeuten. 
Fiir den Rand y = b lJauten die Randbedingungen 
w=0, (3) 
we 0". (4) 
Aus (1) erhalt man A = 0, aus (2) B as = , schlieBlich aus (3) und (4) 
2ab 
ae Cof2ab—1 
und 
~~ Cvj ab 
Dee Gino b* 


1 A, Pucher, Beitrage zur angewandten Mechanik (Federhofer-Girkmann-Festschrift) S. 306, Wien 1950; 
E. Haeussler, Umfangsgelagerte Platten, Diss, Hannover 1951. 
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Fiir die weiteren Untersuchungen benétigt man vor allem den Neigungswinkel der Biegeflache 
langs des Randes, an dem das Moment angreift: 


By =w 


mit 
Wy 0 = ee (C+D) cos «x dx 
und } 
c+ ree ea 
damit wird die Randneigung 
00 


1 f-m(a)Gin2ab—2a0b F 
Uy 0 = Fae | oe SOOT b= leah ek ae 
0 


3. Einseitig eingespannter Plattenstreifen. Die Biegeflache eines beliebig belasteten, frei dreh- 
bar gelagerten Plattenstreifens 1aBt sich ebenfalls mit Hilfe Fourierscher Integrale bestimmen. 
Dies soll an einem einfachen Beispiel gezeigt werden. 

Eine konstante Streckenlast p von der Linge 2c laBt sich in folgende Form entwickeln: 


co 


pee 2 [po cos xx dx , 
0 
wobei fiir den beispielsweise gegebenen Lastfall p(«) = 2p/a sin ac bedeutet. Die Biegeflache 
wird in zwei Bezirke unterteilt und zwar w, von y, = 0 bis y, = b, und w, von y, = 0 bis y, = by 
(Abb. 2). Da die beiden Bezirke lastfrei sind, muB die Gleichung der Biegeflache wieder der homo- 
genen Plattengleichung geniigen. Fir jeden der beiden Bezirke wird derselbe Lésungsansatz wie 
in Ziff. 2 gewahlt. Die acht Freiwerte werden mit A,, B,, C,, D,, A,, B,, C, und D, bezeichnet. 
Zu ihrer Ermittlung stehen vier Rand- und vier Ubergangsbedingungen zur Verfiigung. Am Rand — 
yi = O ist w, = 0 und w, = 0, am Rand y, = 0 ist w. = 0 und w, = 0. Fir y, = 5, und y, = b, 
gelten die Ubergangsbedingungen 


W, = We > 
1 mel 
Ww, =U,, 


GW, + Wy, = Px» 


Abb, 2. Linienlast bei y, = b, yz = 6). I =— N (w” + 1 ) 
bedeutet. 
Aus den Randbedingungen erhalt man A, = B, = 4, = B, = 0. Die Ubergangsbedingungen 
ergeben die vier Bestimmungsgleichungen fiir die restlichen Freiwerte. Es wird 


p(%) 4b, Cojab, Sina b + Sina b, Gina b—ab Gina b, 


Gane Cof2ab—1 : 
(pie p(s) Ginab, Ginab 
1 «Noe Cof2ab—1 * 


Die Werte C, und D, erhalt man aus C, und D, durch Vertauschung von b, und b,. Der Neigungs- 
winkel langs des Randes y, = 0 ist wieder 


co 
By. (4, =9) = My, <9 = ee (C, + D,) cos xx da , 


p(%) ob, Gof ab, Sina b—ab Gina b, } 


Cee asa: Cof2ab—1 : 


damit erhalt man die Randneigung zu 


foe) 
: 1 P(a) & b, Coj a b, Sina b—ab Gina b 
Wino = zy | = a Cop2ab 1 * cos ax da. 
6 
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Am eingespannten Rand muB die Neigung null sein, d. h. es mu& dort ein Moment angreifen, 
das die entgegengesetzt gleichgroBe Neigung hervorruft, so daB kommt 


B, + Bm = 9. 


Ist p(x) gegeben, se ist dies die Bestimmungsgleichung fiir m(x), woraus sich das Einspannmoment 
errechnen 1aBt: 
2 p(a) a b; Coj a b, Gina b —a b Gina b, 


Mu) x Gin 2ab—2ab 
Das Einspannmoment ist un 
co 
1 
Vie = [ m(x) cos xx dx 
0 
[oe} 
Tey: * p() & b, Coj a b, Sina b—a b Gina by ae 
= 1 oy Gin2ab—2Qab — ee ee 
0 
4. Zahlenbeispiel. Es sei 
2, 5 
P(x) = =? sin xc. 


Die Abmessungen der Platte seien 6, = b, = 6/2, c = b = 4,0 m. Gesucht sei das maximale Ein- 
spannmoment an der Stelle x = 0, so daB cosax = 1 ist. Dann wird 


co 


Ve 4p [ sn 4a 2 of 2g on 4a —4 Gin 2a ob, 
I a Gin 8 a — 8a 


6 
Das Integral konvergiert ; denn fiir lim « +0 nahert sich der Integrand dem Wert 1,5, fiir lim «~ —> 00 
wird der Integrand null und die Funktion ist itiberall im Intervall stetig und nimmt periodisch 
positive und negative Werte an, deren Maxima mit wachsendem x kleiner werden. Die numerische 
Auswertung des Integrals mit 4p/2 multipliziert ergibt 


M, = — 0,5789 p . 


5. In Querrichtung durchlaufender Plattenstreifen. Hier gelten grundsatzlich die gleichen 
Uberlegungen wie in Ziff. 3. An einem Beispiel einer Zweifeldplatte soll der Rechnungsgang vor- 
gefiihrt werden. Die Abmessungen und die Belastung des einen Plattenfeldes seien die gleichen 

wie in Ziff. 4, das zweite sei durch eine Flachenlast beansprucht. Durch den Ansatz 


w= f[W,(1 + AGofay tay BGinay + C Ginay 


+ ay D Eo} ay) cos xx dx 


wird die inhomogene Plattengleichung befriedigt, wenn 


Poy = = | po) cosaxdx und W, = Arak 
0 


bedeutet. Aus den Randbedingungen ergeben sich die Frei- 
 -werte 


Piss ies 
1 
B — oy > 
(ie (a b — 2 Gina b) (1 — Goj a b) 
— Gof 2Qab—1 , : im 
y=0 o—____* 
Gin « b (1 — oj « b) a x 
D= labo) aa Abb. 3. Uber Re rece cmendce 


15 
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Es wird wieder die Randneigung am Rand y = 0 bestimmt: 


Bop 20) = Wye ie W,, (C+D) cos xx dx 


1 fH) (2b Citta) (les Sones) 
aN [ of «Of dab 1 DUCE ie 
0 


Bei der durchlaufenden Platte muf die Biegeflache tiber der Auflagerlinie, der Linie y, =y3; = 0 
stetig gekriimmt sein, d. h. es muB sein (Abb. 3) 


Bo + Bos + Bat + Bus = ()) 


Unter der Voraussetzung gleicher Biegesteifigkeiten und gleicher Stiitzweiten beider Platten- 
felder ist 8,4 = Bm3 = Pm, und die Bestimmungsgleichung fiir m(x) lautet 


= Pi, = Bo + Bos » 


so dah 
1 7 ma) Sin2ab—2a6 d 
an | ca EON a bg 
0 
1 ~ il f ; ; 
= | sae jt) (xb, Cof xb, Gin ab — xb Sin xb,) 
: 
+ Ee (ab — Gin xb) (1 — Co} xb) cos «x dx 
P(d) (a b, oj a b, Gina b —a b Gina b,) + Ee (a b — Gin a b) (1 — Go} « b) 
WC) . earn; ~a (Gin 2ab—2ab) 


wird. Das Biegemoment tber der Auflagerlinie 1aBt sich hieraus errechnen zu 


foe) 


: 1 
M, = = [ mo) cos ax dx. 
0 
6. Zahlenbeispiel. Mit denselben Abmessungen wie unter 3. und mit c, = ¢3, 


co 


il 2 F 
Pi Pie | Po) cosaxdax, p(x) = —-sin xe 
0 
wird 
[e.e) 
M, a 2p [ sin 4.0 
I o 
0 
Es soll das Moment an der Stelle x = 0 errechnet werden, dann ist cos xx = 1, das Integral kon- 
vergiert, lim « —>0 nahert sich dem Wert 5,5. Die numerische Auswertung ergibt 


M, = 1,4324 p. 


7. Zusammenfassung. Die Beanspruchungen und Verformungen eines Plattenstreifens unter 
dem Angriff eines Randmomentes wurden untersucht. Das Randmoment wurde in ein Fourier- 
sches Integral entwickelt. Die gesuchten GréBen ergeben sich als bestimmte Integrale, die numerisch 
ausgewertet werden. Insbesondere wurde die Neigung der Biegeflache am betrachteten Rand 
untersucht. Fiir eine gegebene Belastung, die ebenfalls in Fouriersche Integrale entwickelt wurde, 
wurde am betrachteten Rand die Neigung der Biegeflache bestimmt. Aus der Bedingung einer 
horizontalen Tangente der Biegeflache bei starrer Einspannung, bzw. der Stetigkeit der Biege- 
flache iiber der Auflagerlinie bei durchlaufenden Platten konnte durch Vergleich der Randneigun- 
gen das Kinspann- bzw. Stiitzmoment ermittelt werden. Anhand von Beispielen wurden der Rech- 
nungsgang erlautert und Ergebnisse gewonnen. 


b, Coj a b, Ginab—b Gina b, | (ab — Gina b) (1 — Gof a b) 


~~ Sin2ab—2ab o® (Gin 2a b—2a8) cos ax dx. 


(Eingegangen am 4. September 1956.) 


Anschrift des Verfassers: Dipl.-Ing. Fritz W. Mader, Hannover, Kortumstr. 12. 
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Uber die Geschwindigkeitsverteilung in Wirbeln* 
Von A. Timme 


1, Einleitung. Bei vielen Strémungsvorgangen lassen sich eigenartige kreisende Bewegungen 
des strémenden Mediums beobachten, die Sifenbar einer heaondlecen Cescrmabioneit Laie 
In der beschreibenden Physik werden solche Flissigkeitsbewegungen mit Wirbel bezeichnet, 
Volksmund oft auch Strudel genannt!. Schon im Altertura fader man bei Homer eine Beschr bite 
des Strudels Charybdis, und im Mittelalter befaBte sich vornehmlich Leonardo da Vinci mit wirbel- 
artigen Ablésungserscheinungen?. 

Will man spezielle Strémungsprobleme, bei denen Wirbel auftreten, wie z. B. die Karmdnsche 
WirbelstraBe oder viele Ubergangsstadien zwischen laminarer und turbulenter Strémung theore- 
tisch naher erfassen, so ist eine genaue Kenntnis der Geschwindigkeitsverteilung in einem einzelnen 
Wirbel von besonderer Wichtigkeit. Im Rahmen des Forschungsprogramms des Hermann Fét- 
tinger-Institutes, das sich vornehmlich mit diesem Gebiet befaBt, sollten im AnschluB an die 
theoretischen und experimentellen Arbeiten von U. Domm? iiber die Stabilitat von WirbelstraBen 
neue Untersuchungen hieriiber durchgefiihrt werden. 


Eine Karmansche WirbelstraBe entsteht, wenn ein zylindrischer Kérper quer zu seiner Langs- 
achse angestrémt wird. Es lésen sich bei bestimmten Anstrémgeschwindigkeiten, die Reynolds- 
schen Zahlen Re>40 entsprechen, hinter dem Kérper periodisch abwechselnd rechtsdrehende 
und linksdrehende Wirbel ab, die sich stromab zu der charakteristischen doppelreihigen Wirbel- 
straBe anordnen (Abb. 1 und 2). Dieser Vorgang laBt sich im Versuch verhaltnismaBig leicht 
darstellen, so daB es nahe lag, gerade bei der Kérmanschen WirbelstraBbe die Verhaltnisse im Wirbel 
und die Wirkung von Wirbeln aufeinander zu untersuchen. 


2. MeByerfahren. a) MeBprinzip. Ein zylindrischer Kérper wurde von einem Schleppwagen 
durch einen mit Wasser gefiillten Kanal gezogen. Die hinter dem Zylinder entstehende Karmdnsche 
WirbelstraBe wurde durch Bestreuen der Wasseroberflache mit Aluminiumpulver kenntlich 
gemacht und von zwei iiber dem Kanal befindlichen Kameras in zeitlich festgelegter Reihenfolge 
photographiert. In diesen Photographien konnte aus der Lange der Bahnkurven der Aluminium- 
teilchen, dem Abbildungsmafstab und der KameraverschluBzeit fiir jeden Ort der Betrag und die 
Richtung der Teilchengeschwindigkeit bestimmt werden. Mit der Annahme, da®B die Aluminium- 
teilchen wahrend der Beobachtungszeit alle Bewegungen der Wasseroberflache mitmachen, erhalt 
man somit auch die Geschwindigkeitsverteilung im Strémungsfeld an der Wasseroberflache. 


b) MeBanordnung. Der fiir die Versuche beniitzte Kanal war 2m lang, 50 cm breit und 
40 cm tief. An der einen Schmalseite hatte er ein Uberlaufwehr, damit man nach jedem Versuch 
bequem die Wasseroberfliche erneuern konnte. Neben dem Kanal lief der Schleppwagen auf 
einem separat gelagerten Breitflanschtrager. An einem seitlichen Ausleger des Wagens konnten 
Zylinder verschiedenen Durchmessers befestigt werden. Wahrend der Meffahrt betatigte der 
Wagen zwei Schleifkontakte, die elektromagnetisch eine Stoppuhr fiir die Geschwindigkeits- 
bestimmung ein- und ausschalteten. Mit zwei weiteren Kontakten wurden die Verschliisse der 
Photoapparate ausgelést. 

Als wesentliche Vorbedingung fiir die Reproduzierbarkeit der Untersuchungen war der gleich- 
maBige Lauf des Schleppwagens zu fordern. Vorversuche mit groBdimensioniertem Gleichstrom- 
motor; Gewichtsantrieb mit Reibungs- oder Strémungsbremse oder ihre Kombinationen mit 
einem Synchronmotor konnten nicht befriedigen. Lediglich der Betrieb mit einem geniigend 


_kraftigen, nicht selbstanlaufenden Synchronmotor zeigte bei stroboskopischer Uberpriifung eine 
-ausreichende Konstanz der Wagengeschwindigkeit im Bereich von 5 cm/s. Fiir noch kleinere 


Geschwindigkeiten gelangte ein selbstanlaufender Synchrondoppelmotor der Type AEG DSSE 1/4 s 


zum Einbau. 
* Von der Fakultat fiir Maschinenwesen der Technischen Universitat Berlin genehmigte Dissertation. 
Berichter: Prof. Dr.-Ing. R. Wille, Prof. Dr.-Ing. H. Amtsberg; eingereicht am 5. 4.56, mdl. Priifung am 23. 7.56. 
1 A. Sommerfeld, Mechanik der deformierbaren Medien, S. 34, Leipzig 1955. 
2 Hdb. d. Exp.-Physik Bd. IV, 2. S.7 ff. 
3 U, Domm, Ing.-Arch. 22. 6 (1954), S. 400. 
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V r otors 
erschiedenartige Fehlerquellen, wie z. B. Drehschwankungen des OBO OE a , 
; 1 “* ee 5 o 

auf unsaubere Herstellung der Zihne am Rotor zurickzufihren waren, Kigenschwingungen 


di p 
my Systems Motor-Schleppseil (0,3 mm Durchmesser, Stahldraht) — Wagen, Teilungsfehler und 


des 


Abb, 2, Wirbel aus einer Kdrmdnschen WirbelstraBe, MaBstab 1,67: 1. 


Spielschwankungen des Flanschgetriebes (3000:4) am kleinen Motor, konnten durch geeignete 
Mafinahmen cum grano salis beseitigt werden. 
Die beiden Plattenkameras der Typen 


1. Schneider-Kreuznach Doppel Anastigmat Isconar f = 4,5 13,5 cm, 


2. Rhodenstock ee re Eurynar f = 6,3 13,5 cm 


XXV. Band 1957 Timme: Uber die Geschwindigkeitsverteilung in Wirbeln 207 


waren an ein verstellbares Geriist tiber dem Kanal montiert. Ihre VerschluBéffnungszeit lieB sich 

aus der Strichlange eines am Schleppwagen befestigten Wegzeigers bestimmen. Um aber von 

zufalligen Geschwindigkeitsschwankungen des Wagens wahrend der Aufnahme unabhangig zu 
sein, wurde spater eine schnellaufende Synchronuhr mit von innen leuchtendem Zifferblatt und 
Zeiger mitphotographiert. Dazu stand die Uhr seitlich am Kanalrand. Durch ein Linsensystem 
und einen dicht tiber der Wasseroberflache angebrachten Spiegel wurde das verkleinerte Bild 
der Uhr in die Objektebene der Kameras hineinprojiziert. Schwierig war die Justierung des 
Systems deshalb, weil durch die grofe Parallaxe der beiden nebeneinander angebrachten Kameras 
die bendtigte Spiegelflache leicht bei dem inneren Apparat Teile des MeBraumes abdeckte. 


Mit dieser Anordnung war es nicht nur méglich, die Verschluéffnungszeit T genau zu be- 
stimmen, sondern auch noch die Zeitdifferenz At zwischen den beiden Aufnahmen zu iiberprifen. 
Die VerschluBausiésung erfolgte, wie schon erwahnt, elektromagnetisch. Die Zeit At registrierte 
eine zweite Stoppuhr ebenfalls elektromagnetisch. Um die Anzugszeit der Auslisemagnete, die 
bei beiden Apparaten verschieden sein konnte, auszuschalten, schlossen die Magnete erst in dem 
_ Augenblick einen Relaiskontakt, indem der Verschlu8 tatsachlich betatigt wurde. Die Justierung 
war jedoch sehr schwierig, wie sich aus Uberpriifungen mit der Synchronuhr ergab, zumal die 
_Compurverschliisse sowohl bei der Stellung M als auch T benutzt wurden, so daB schlieBlich nur 
noch der SynchronuhrmeBwert verwandt wurde. Dicht iiber der Wasseroberflache befanden sich 
noch zwei MeSBmarken als Langenmaftstab fiir die Auswertung photographischer Vergréerungen. 


c) Die Messung und ihre Auswertung. Etwa 5 bis 10 Stunden vor Beginn der Messung 
wurde in den Kanal frisches Wasser nachgefiillt, bis der Wasserspiegel das Uberlaufwehr iiber- 
spulte. Hierdurch erhielt man eine einwandfrei saubere Oberflache. Danach durfte die Wasser- 
-masse nicht weiter gestért werden. Die MeBfahrt selbst begann mit dem Augenblick der Be- 
streuung. Der Wagenantriebsmotor war so einzuschalten, da die Auslésung der Photoverschliisse 

zwischen 10 bis 40 sec nach dem Aufbringen der Bestreuung erfolgte. Zu einem spateren Zeitpunkt 
_konglomerieren die Aluminiumteilchen, so da man kein feinmaschiges Bahnliniennetz erhalt. 
War der Versuch stérungsfrei abgelaufen und das Plattenpaar auswertbar, so wurden Vergrébe- 
-rungen im Mafstab 1:1 zur Wirklichkeit hergestellt. In vorherigen Testaufnahmen waren im 
_Zusammenhang die Kameraobjektive, das VergréBerungsobjektiv und die Schrumpfung des 
Photopapiers durch den TrocknungsprozeB iiberpriift worden. Sowohl hierbei als auch bei den 
_spateren MeBaufnahmen wurde die Parallelitat der Plattenebene zur Objektebene mittels Wasser- 
_waagen kontrolliert. Trotz einer geringen kissenférmigen Verzeichnung der Objektive von 0,3% 
bzw. 0,07°% und einer linearen Verzerrung von 0,5°% bzw. 0,06°% wurden diese Fehler als vernach- 
_lassigbar klein angesehen. 

Zur Kontrolle der Schrumpfung des Papieres, die von der Temperatur der Trockenpresse ab- 
hangt, wurde spater im Kontaktverfahren ein durchsichtiger MaBstab mitkopiert. Fiir die Einzel- 
auswertung geniigte dies jedoch nicht, sondern dazu mute jeder Wirbel extra im Mafstab 5:1 
zur Wirklichkeit (die Platte etwa l5fach) vergréBert werden. Sodann wurde in die VergréBerung 
in OriginalgréBe als Koordinatensystem die StraBenachse und senkrecht dazu Achsen mit Zenti- 
-metereinteilung durch die scheinbaren Wirbelmittelpunkte eingezeichnet. Nach Ubertragung 

dieser Achsen in die EinzelvergréBerungen konnten nun die Strichlangen der Teilchen entlang 
den Querachsen ausgemessen werden. 

Wegen der Unsymmetrie des Rotationsfeldes wurden nur Teilchenbahnlinien durch die Achse 
oder in unmittelbarer Nahe beriicksichtigt. Im Laufe fortschreitender Versuchstechnik muBte 
festgestellt werden, daB die Papierschrumpfung durch Feuchtigkeits- und Temperaturschwan- 
kungen im Laborraum zu merkbaren Fehlern fihrte. Infolgedessen konnten die PapiervergréBe- 
rungen nur noch zur Orientierung der Koordinatenachsen benutzt werden, wahrend die eigentliche 

Auswertung direkt unter dem VergréBerungsapparat erfolgte. 


3. Mefergebnisse. Den im folgenden Kapitel zu interpretierenden MeBergebnissen legen 
die Erfahrungen aus weit iiber 100 MeBanlaufen zugrunde, welchen ungezahlte fehlgeschlagene 
Vorversuche vorausgingen. Es wurden etwa 70 Wirbelprofile genauer ausgewertet, wobei die Aus- 
-wertung eines Wirbels mit dem Profil zur Zeit t, und dem Profil zur Zeit t, mehrere Stunden Dunkel- 
kammerarbeit erforderte. Der notwendige Zeitaufwand lieB eine gewisse Einschrankung in der 
Variation der Versuchsparameter geboten erscheinen. 

a) Bereich Re = 200. Im ersten Versuchsabschnitt wurden die Messungen mit folgenden 
Versuchsbedingungen durchgefiihrt: Zylinderdurchmesser D =2cm, Zylindergeschwindigkeit 
U = 1,018 cm/s. Es ist iiblich, die eine Karmansche Wirbelstrabe charakterisierende dimensions- 
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lose Kennziffer, die Reynolds-Zahl, mit diesen Werten zu bilden: 


UD 
KG = 5 


v 


Fiir das untersuchte Gebiet ist damit 


ine — 200% 


Die Bildung der Wirbel hinter dem Zylinder und ihre Anordnung zu einer doppelreihigen Wirbel- 
straBe ist selbstverstandlich bei Anlegung des hier erforderlichen scharfen MaBstabes ein stérungs- 
empfindlicher Vorgang. Haufig war die StraBe aus unbekanntem Anlafs unregelmaBig. Dadurch 
waren mehr als 2/, aller Aufnahmen nicht auswertbar. Die Ursache dieser Stérungen kann nicht 
mit Sicherheit angegeben werden, sie scheint jedoch von der Wartezeit zwischen zwei MeBanlaufen 
abhangig zu sein. Bei den besten Auf- 
Umfangsgeschw. |cm/sek| nahmen wurde bis zu 19 Stunden die 
ee Beruhigung der Wassermasse ab- 
gewartet. So geniigt vielleicht schon 
ge eine geringe Reststr6mung innerhalb 
der untersuchten Wassermassen, um 
die WirbelstraBe unter Umstanden 


Sali erheblich zu stéren. 
Radius|em| Als Beispiele fiir die durch die 
a > 7 Auswertung erhaltenen Wirbelprofile 
; seien hier nur je ein Profil eines jun- 
22 gen Wirbels (Abb. 3) — das ist ein 
Abb. 3. Geschwindigkeitsverteilung in einem Wirbel Re, = 200. Wirbel in der Nahe des Zylinders — 
Erster Wirbel hinter dem Zylinder. ‘ und eines alteren Wirbels (Abb. 4) 


wiedergegeben. In ihren wesentlichen 
Merkmalen sind sie sich alle 
untereinander ahnlich. 
Auffallend ist zunachst die 
Unsymmetrie, welche auf die 
Eigenbewegung der Wirbel, 
te eee haufig mit Nachlaufgeschwin- 
digkeit bezeichnet, zuriickzu- 
fiihren ist. Dadurch gibt die 
Stelle mit der Geschwindigkeit 
Null, der scheinbare Mittel- 
punkt, nicht den wahren Wir- 
belmittelpunkt an. Sondern 
dieser ist in der Mitte zwischen 
den beiden Stellen maximaler 
Abb. 4. Geschwindigkeitsverteilung in cinem Wirbel Re, = 200. Umfangsgeschwindigkeit an- 
Fiinfter Wirbel hinter dem Zylinder. zunehmen, wo er sich mit dem 
arithmetrischen Mittel der 
maximalen Umfangsgeschwindigkeiten fortbewegen wird. Von dieser Stelle aus mégen im folgenden 
auch die Geschwindigkeiten und Radien gemessen werden. 


G 


{ Umtangsgeschw.{cm/sek| 
G6 


Radius{cm] 


In der Nahe des Wirbelmittelpunktes wachst die Umdrehungsgeschwindigkeit etwa linear mit 
dem Radius r, um nach einem Ubergangsgebiet ungefahr mit der Potenz 1:r wieder abzufallen. 
Die zeitlichen Anderungen, die durch die beiden Aufnahmen hintereinander festgestellt werden 
konnten, zeigten allgemein einen Abfall des Geschwindigkeitsanstieges im Wirbelmittelpunkt. 
Die maximalen Umfangsgeschwindigkeiten nehmen mit der Zeit ab, wahrend der zugehérige Radius, 
den man auch mit Kernradius zu bezeichnen pflegt, gréBer wird. Es bietet sich somit das Bild 


eines laminaren Wirbels, der sich auf Grund der Zahigkeitswirkung ausbreitet und damit einem 
Ausgleichsgesetz unterliegt. 


b) Bereich Re = 1000. Im zweiten Versuchsabschnitt wurde derselbe Zylinderdurchmesser 
D =2cm verwendet, die Zylindergeschwindigkeit auf U = 5,24 em/s erhéht. Somit war die 


Reynolds-Zahl 
Re = 1048 , 
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Diesmal erwies sich der Vorgang der StraBenbildung als noch wesentlich stérungsempfindlicher 
als im Bereich der kleinen Re - Zahl. Sehr haufig wuchsen kleine Stérungen oder UnregelmaBig- 
keiten derart an, da®B schlieBlich die ganze StraBe zerfiel. Auch die ausgewerteten Wirbelprofile 
zeigten Stérungsmerkmale, die in dem vorhergehenden Abschnitt nicht beobachtet wurden. 

In den beigefiigten Wirbelprofilen (Abb. 5 bis 8) lassen sich folgende charakteristische Er- 
scheinungen beobachten: Auf dem linken Ast von Abb. 5a (Zeit t,) zeichnet sich eine Ausbuchtung 
ab, die sich bei Abb. 5b (Zeit t,) wieder dem Gesamtprofil eingefiigt hat. In Abb. 6a und Abb. 6b 
ergibt sich der zeitlich umgekehrte Verlauf. Also Stérungen, die auftreten und wieder abklingen 
kénnen (dieser Effekt ist aus der Drehung des Wirbels allein noch nicht zu erkliren, wie man 
aus einer zahlenmaBigen Nachrechnung ersehen kann). 

Das Profil in Abb. 7 ist zu beiden gemessenen Zeiten stérungsfrei. Weiterhin konnte verschie- 
dentlich beobachtet werden, wie ein Wirbel von regelmaBigem Aussehen sich in drei Einzelwirbel 
auflést. Den umgekehrten Vorgang zeigt Abb. 8. In Abb. 8a ist der Wirbel noch sehr gestért, 
um dann in Abb. 8b wieder die Charakteristik eines sehr jungen Wirbels — wenigstens in Kern- 

-nahe — zu gewinnen. Sieht man von diesen Stérungserscheinungen ab, so gelten auch fir die 
Wirbel dieses Re-Bereiches die oben angegebenen Merkmale eines Ausbreitungsvorganges. 


4. Vergleich der MeBergebnisse mit der hydrodynamischen Theorie. Eine exakte Lésung 
der Navier-Stokesschen Differentialgleichung konnte von C. W. Oseen! und G. Hamel? fiir den 
'“Spezialfall der ebenen Strémung angegeben werden, sofern rotationssymmetrische Randbedin- 
-gungen vorliegen. Fiir die reine Wirbelstrémung wurde als Anfangsbedingung die Geschwindig- 
_keitsverteilung des Potentialwirbels 


Cee C0 (1) 


290 Ti 


-eingesetzt. Damit ergab sich fiir die reibungsbehaftete Fliissigkeit 


fe eae C= 0. (2) 


Diese Gleichung beschreibt qualitativ tatsaichlich das zeitliche Verhalten eines Wirbels so, wie es 
auch in den Experimenten beobachtet werden konnte. Daher erschien es sinnvoll, diese Gleichung 
auch fiir einen quantitativen Vergleich mit den gemessenen Wirbelprofilen heranzuziehen. 


Nun sprach U. Domm?® im Anschlu8 an seine Stabilitatsrechnung die Vermutung aus, daB die 
fiir das zeitliche Wachsen der Wirbel verantwortliche GréBe y einen gréBeren Zahlenwert als den 
der laminaren Zahigkeit haben miiBte. H. B. Squire’ berechnete die Geschwindigkeitsverteilung 
‘eines Wirbels unter Beriicksichtigung turbulenter Schwankungen und erhielt das laminare Ge- 
-schwindigkeitsgesetz, bei dem das y um eine additive GréBe erhéht war. 
| Um einen solchen Effekt zu untersuchen, hat man in (2) neben den beiden GréBen T' und t 
auch den Wert von y als Unbekannte anzusehen. Mit einem Geschwindigkeitsprofil allein ist die 

Bestimmung dieser drei GréBen noch nicht méglich, erst die Aufnahme eines zweiten Wirbel- 
profils nach einer Zeit At erméglicht die Berechnung auf Grund folgender Beziehungen: 


) Der Geschwindigkeitsanstieg in der Nahe des Wirbelmittelpunktes (r = 0) ergibt sich zu 


; 0 Cy ie a if! 3 
) Eis “oo Bvt’ (3) 
hus den Geschwindigkeitsanstiegen zu den Zeiten t, und ¢, folgt mit t, —t, = At 

= 8a At (4) 


y Vee Ve 


Poe ‘Por 
Den Radius r, der maximalen Umdrehungsgeschwindigkeit ¢,,,, erhalt man nach kurzer numerischer 
m 
-Rechnung mit 


ite 1,2564 (5) 
4yvt 


W. Oseen, Ark. f. Mat., Astron. och Fys. 7 (1911); Hydrodynamik, S. 82., Leipzig 1927. 
Hamel, Jahresb. Deutsch. Math.-Vereinig. 25 (1917), S. 34. 


ehe FuBnote 3 von Seite 205. : 
I B. Squire, National Advisory Committee for Aeronautics 16, 666 F.M. 2053 vom 18. 3. 1954. 
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cebildet. Mit Hilfe dieser drei Bezichungen kann man die wirksame kinematische Zahigkeit sowohl 
force (6) als auch durch (4) und (5) bestimmen, ebenso kann die Wirbelstarke mit (4) und (6) 
oder (6) und (7) oder (4) und (7) berechnet werden. Weiterhin 1aBt sich die Wirbelstarke noch aus 
dem Stromlinienbild durch graphische Integration um den Wirbelkern recht sicher ermitteln. 
Die auf die angegebene Weise berechneten Wirbelstarken wiesen auf prinzipielle Unterschiede 
von Wirbeln aus dem Re-Bereich um 200 und dem Bereich um Re = 1000 hin. Im Re-Bereich 200 
waren alle Wirbel fast gleich stark, und die fiir denselben Wirbel nach den verschiedenen Methoden 
ermittelten Werte unterschieden sich im allgemeinen nur um ein geringes. Im Mittel betrug die | 
Wirbelstarke | 
I =5+ 0,8 cm/s” bei Re = 200. 
Demgegeniiber ergaben sich bei der Bestimmung von /’ im Bereich Re = 1000 von Wirbel zu 
Wirbel wesentliche Unterschiede, ttber die jedoch keine GesetzmaBigkeit gefunden werden konnte. 
So war es beispielsweise méglich, daB der erste Wirbel hinter dem Zylinder die Starke I = 20 ema 
der nachstfolgende J’ = 10 cm?/s und der dritte /’= 15 cm?/s hatte. Auch die Wahl der Be- 
stimmungsgleichungen brachte Differenzen bis zu 30% mit sich. So kann man hier dem Mittelwert 
[=17,5+7,5 cm2/s bei Re = 1000 


nur die Bedeutung einer numerischen Mittelung von MeBergebnissen beimessen. 


Die gleiche Tendenz war auch bei der Berechnung des y- Wertes zu finden. Im Bereich Re = 1000 
waren die Unterschiede in den Einzelmessungen wieder so groB, da man jedem Wirbel nur ein 
individuell giiltiges » zur Zeit der Messung zusprechen kann. Der gemittelte y-Wert betrug 


vy = 0,09 cm?/s bei Re = 1000, 


in den Grenzen 


vy = 0,02 bis »y = 0,2 cm?/s 
gegeniiber dem mit einem Standardviskosimeter gemessenen 
Viaminor = 0,01 cm?/s_ fiir Wasser bei 20° C. 
Dagegen stimmten im Bereich Re = 200 die berechneten y-Werte innerhalb der Fehlergrenze mit 
dem laminaren Wert tiberein: 
y = 0,01 + 0,003 cm?/s = Yjaminar bei Re = 200. 
Dies gilt jedoch nicht fiir verhaltnismaBig junge Wirbel, d.h. fir Wirbel in Zylindernahe (etwa 
erster und zweiter Wirbel hinter dem Zylinder). Hier erhéht sich der berechnete Wert bis zu 
y = 0,03 cm?/s bei Re = 200. 


AnschlieBend sei zur Beurteilung des Gewichts von Einzelmessungen noch das Ergebnis einer 
Fehlerrechnung mitgeteilt, die fiir einen speziellen Wirbel durchgefiihrt wurde. Als Fehler wurden 


5% der den Kurven zu entnehmenden Werte angesetzt, mit denen der maximale Fehler der Rechen- 
gréBen im jeweils ungiinstigsten Fall 


ieee 60% [vgl. (4)], 
= = + 65% [vgl. (6)] , 

Ar) =+ 40% [vgl. (7)] 
T? 


betragt. Lediglich bei der graphischen Integration der Wirbelstarke liegt der MeBfehler bei 


AML 
= = a Bh 


wobei hier allerdings ein systematischer Fehler hinzu kommt durch die raumliche Ausbreitung 
der Wirbelstarke. Ist diese soweit fortgeschritten, daB sie sich mit den Nachbarwirbeln iiber- 
schneidet, wird die so gemessene Wirbelstarke zu klein ausfallen. Fiir die angegebenen Mittelwerte 
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diirfte der Fehler jedoch wesentlich niedriger sein, da sie sich jeweils aus einer gréBeren Anzahl 
ausgewerteter Wirbel ergeben. 


5. Deutung der Ergebnisse. Der Vergleich der experimentell gefundenen Wirbelprofile mit 
dem Geschwindigkeitsverteilungsgesetz von Hamel-Oseen ergab Ubereinstimmung in den drei 
wesentlichen Merkmalen des realen Wirbels, dem zeitlichen Kleinerwerden des Cesehindiet i: 

-anstieges im Wirbelmittelpunkt, dem Sinken der maximalen Umfangsgeschwindigkeit a dem 
_Wachsen des zugehérigen Radius mit der Zeit. Die hiermit berechneten /- und y-Werte stimmten 
fiir den betreffenden Wirbel unabhangig von den fiir die Berechnung herangezogenen Merkmalen 
tiberein, wenn auch nicht exakt, so doch gréBenordnungsmaBig. Daraus laBt sich schlieBen, dab 
(2) das Verhalten eines realen Wirbels im Prinzip richtig beschreibt, und daB zweitens die feler 
chung (2) auch bei Anwesenheit turbulenzartiger Stérungen giiltig ist, sofern man fiir das Symbol 
y einen geeigneten Wert einsetzt. In der Navier-Stokes-Gleichung hat » die Bedeutung der kine- 
“matischen Zahigkeit gemaB dem Newtonschen Schubspannungsansatz 
_ du 
T=V7O0 a 
wobei » die Eigenschaft eines Stoffbeiwertes hat. 


Aus einer Veranderung des »-Wertes folgt, da der Newtonsche Ansatz bei Einsetzen turbulenz- 
artiger Stérungen seine Giiltigkeit verliert, bzw. daB dem y keine echte Stoffeigenschaft zukommt. 
Soll der Newtonsche Ansatz per definitionem weiter gelten, wie es in der Turbulenztheorie itiblich 
ist, so erhalt man zusatzlich zu der laminaren Zahigkeit die ,,scheinbare Zahigkeit“. 


Das Verhalten der Wirbel kann man mit dieser Festlegung folgendermaBen deuten: 

Der Wirbel entstehe durch Ablésen und Aufrollen der Reibungsschicht des Zylinders. Beide 
Vorgange bergen den Keim ungeordneter Bewegung in sich, so daB der gerade entstandene Wirbel 
mit kleinen Stérungen behaftet ist, welche eine scheinbare Zahigkeit hervorrufen. Im Bereich 
Re = 200 klingen diese Stérungen rasch ab, so da bei alteren Wirbeln nur noch die laminare 
Zahigkeit wirksam ist. Dagegen sind im Bereich um Re = 1000 die Stérungen entweder von 
vornherein gréfer, oder sie wachsen mit der Zeit noch an. Diese kleinen Geschwindigkeitsschwan- 
kungen kénnen weiterhin zu endlichen Stérungen entarten, wie sie auf den Wirbelprofilen Abb. (5), 
(6) und (8) zu beobachten sind. Die endlichen Stérungen kénnen dariiber hinaus zu einer Auf- 
spaltung des Wirbels in mehrere Einzelwirbel fiihren. Die eigentliche Ursache dieser endlichen 
Stérungen und eine Gesetzmafigkeit tiber ihre Entstehung konnte experimentell noch nicht 
gefunden werden. 

Zur Erfassung des Wesens der endlichen Stérungen sei ein Gedankenmodell gestattet!: Ein 
Wirbel von einer bestimmten Starke an setze sich aus mehreren kleineren, sozusagen Elementar- 
wirbeln zusammen. Das Gefiige sei nur lose und kann sich gegeneinander verschieben. Dies macht 
sich in den Ausbuchtungen der Profile bemerkbar, welche anwachsen und abklingen oder auch zur 
Auflésung in die Elementarwirbel fiihren kénnen. 

Wenn nun tatsachlich solchen Teilbetragen der Wirbelstairke eine Verschiebung gegeneinander 
méglich ist, so ist damit auch ein zusatzlicher Impulsaustausch notwendigerweise verbunden. Das 
bedeutet, daB zu dem normalen Impulsaustausch auf Grund der Thermomolekularbewegung noch 
eine zusatzliche Austauschgré8e hinzukommt. Damit erhalt man einen der vollturbulenten Stré- 
mung entsprechenden Ausdruck fiir den Austausch? 
du 
ay 


9 


A=oP 


nur daB in diesem Falle die Zahl der Austauschvorgange kleiner, die Struktur gréber ist und die 
Individualitat der Fliissigkeitsballen langer gewahrt bleibt. Demnach miBte der Wert fiir die 
wirksame kinetische Zahigkeit ansteigen, jedoch nicht in dem Mafe wie in der vollturbulenten 
Bewegung. Auferdem ist danach zu erwarten, daB das gemessene y von der Individualitat des 
einzelnen Wirbels, d. h. von dem zufalligen Verhalten der Elementarwirbel wahrend der Mefzeit, 
und vielleicht auch von ihrer Anzahl abhangt. Die Messungen im Bereich Re = 1000 mit der von 
Wirbel zu Wirbel verschiedenen Wirbelstarke und dem rund 10fach gréBeren y-Wert lieBen sich 


in eine solche Vorstellung ohne Zwang einordnen. 


1 A, Roshko, Calif. Inst. of Technology. NACA T. N. 2913 (1953). — L. Schiller, A. Naumann, Ing.- 


Arch. 9 (1940) S. 335—. 
2 L. Prandtl, Fiihrer durch die Strémungslehre, S. 112ff. Braunschweig 1949. 
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Ebenso wiirden sich auch die MeBergebnisse der Wirbelstarke einfiigen, da die Zahl der Elemen- 
tareinheiten im Gesamtwirbel nur im statistischen Mittel festgelegt sein mége. 

Hierbei soll der Ausdruck ,,Elementarwirbel nicht so aufgefaBt werden wie etwa die Quanten 
der Mikrophysik mit ihrer konstanten, unteilbaren GréBe, sondern vielmehr als Gebilde gréBSerer 
Stabilitat als die des Gesamtwirbels. 

Damit erhebt sich die Frage von der Stabilitat eines Wirbels. Falls eine dimensionslose Kenn- 
ziffer das Verhalten eines Wirbels ausreichend festlegt und diese Kennziffer die Reynolds-Zahl ist, 
so ware diese zweckmaBigerweise so zu definieren: 


Re A epm _ 21meym (8) 


wv v v 


Gilt weiterhin die Hamel-Oseen-Gleichung (2), so ist 
Di pe wy WTI 


m-pm 


und damit 


Re, =0,23 =. (9) 


Wenn nun wahrend des Entstehungsprozesses zu viele Elementarwirbel zusammengefihrt werden, 
oder mit anderen Worten, wenn die Wirbelstarke zu groB wird, so daB der Wirbel durch eine geringe 
auBere Stérung zum Zerfall gebracht und turbulenzahnliche Strémung erzeugt werden kann, so 
liegen ahnliche Verhaltnisse vor wie z. B. bei der laminaren Rohrstrémung in der Nahe der kriti- 
schen Reynolds-Zahl. 

In Analogie dazu kénnte man nach einer kritischen Reynolds-Zahl fragen, bei der ein Einzel- 
wirbel gerade noch seinen laminaren Charakter bewahrt. Nach den bisher vorliegenden Messungen 
miiBte diese gréBenordnungsmaBig gleich 

Rewicir = 0,23 paz = 570 (10) 
sein. Dieser Wert mu sicher noch, sofern unsere Uberlegungen iiberhaupt zu Recht bestehen, bei 
Vorliegen weiteren Versuchsmaterials genauer untersucht werden. 


6. Kritische Untersuchungen iiber das Auswertungsverfahren. Das in Ziff. 4 beschriebene 
Auswertungsverfahren enthalt zwei wesentliche Vereinfachungen, deren Einflu8 auf die gewonnenen 
Ergebnisse im folgenden abgeschatzt werden soll. Die eine Vereinfachung besteht in der Verwen- 
dung des Hamel-Oseenschen-Geschwindigkeitsgesetzes, welches den wahren Verhaltnissen vermut- 
lich nicht geniigend nahe kommt. Die zweite liegt darin, daB das Induktionsfeld der Nachbar- 
wirbel nur bei der Bestimmung des Wirbelmittelpunktes beriicksichtigt wurde. 


a) Einflu8 des Wirbelmodells auf die MeBergebnisse. Als Grundlage fiir die Auswer- 
tung der gemessenen Geschwindigkeitsverteilungskurven diente, wie schon erwahnt, die von 
Hamel und Oseen berechnete Lisung der Navier-Stokesschen-Gleichung fiir den Einzelwirbel mit 
der Randbedingung, da®B zur Zeit t = 0 die Geschwindigkeitsverteilung eines Potentialwirbels 
vorliege. Diese Liésung ist insofern etwas unbefriedigend, als sie keine Riickschliisse iiber die 
Entstehung und das Verhalten eines sehr jungen Wirbels zulaBt, da die Existenz eines 
Potentialwirbels physikalisch unméglich ist (im Mittelpunkt unendlich groBe Geschwindigkeit 
und unendlich negativer Druck). Im Zusammenhang mit parallel laufenden Demonstrations- 
versuchen, bei denen die Wirbelablésung hinter dem Zylinder vom mitfahrenden System aus 
gefilmt wurde, entstand die Vorstellung fiir zwei erweiterte Wirbelmodelle!. 

1. Kreiswirbelmodell. Seitlich vom Zylinder lése sich eine Wirbelschicht ab, die sich 
stromabwarts kriimmt und nach Abreifen der Schicht vom Zylinderrand einen geschlossenen Kreis 
bildet. In diesem Augenblick bestehe der Kreisumfang aus einer gleichmaBigen Folge von Elemen- 
tarwirbeln mit der Zirkulationsdichte y. Nach W. Miiller® ist die Wirbelstarke fir die langs 
eines Kreises r = a gleichmaBig ausgebreitete Wirbelschicht 


2a (x—acosg)? + (y—asing)? 


CE ba Pees 4vt 
eral ¢ dp . (11) 


eens Naturwiss. 37 (1950), S. 193; H. Kaden, Ing.-Arch.2 (1931), S. 140; H. Drescher, Z. Flugwiss. 4 


? W. Miiller, Einfithrung in die Theorie der zahen Fliissigkeiten, S. 115ff., Leipzig 1932. 
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Setzt man x=rcos # und y=r sin #, so erhalt man nach geeigneter Umformung 


r2 4 a? 
oe te Be ee eaek Jolis yi) (12) 
und mit Hilfe des Stokesschen Satzes die Geschwindigkeitsverteilung 
2 ; ya . aie TG 
a tee aad ke Delia) ar (13) 
j J 


Dieses Integral ist geschlossen nicht lésbar*, und um das numerische Ergebnis méglichst allgemein 
giiltig darzustellen, sei folgende Normierung eingefiihrt: 


2a art (a nell ‘ a 
a aot ht 


| Mit diesen dimensionslosen Koordinaten lautet dann die Hamel-Oseensche Geschwindigkeits- 
verteilung 


V= 


(14) 


P 


1 : 
ee) (15) 
und fiir das Kreiswirbelmodell 
R 
Pe ket 2 RAydR miv-V — o,f 2”! 16 
R o(t ) mi Sad ee ee (16) 


0 


Fihrt man die Bessel-Funktionen tiber in Reihendarstellung und integriert gliedweise, so erhalt man 

1 — (R24 42 uct 7 as TRV es Ae IROUE op PAA 
Mani. * re ee) ar (© me leer ga |e sie ments a 
(17) 


In dieser Darstellung 148t sich leicht ibersehen, daB mit verschwindendem Parameterwert A das 
Wirbelprofil wieder in die Hamel-Oseensche Geschwindigkeitsverteilung tibergeht, was dem Ansatz 
nach auch zu erwarten war. Denn 


A = 0 stellt die Verhaltnisse fiir r 


sehr groBe Zeiten dar, oder den 
trivialen Fall, da8 der Ringhalb- 
messer von Anfang an a = 0 war. 
In Abb. 9 sind fiir verschiedene 
Parameterwerte die numerisch 
berechneten Loésungen aufge- 
tragen. 

Fir groBe A-Werte weist 
das Geschwindigkeitsprofil drei 
Wendepunkte auf, und im Wir- 
belzentrum geht der Geschwin- 


R 
06 V- i Re A 7 fs RAaR — 
0 


digkeitsanstieg fiir kleine Zeiten R 
gegen Null (Abb. 12). Dies sind 

zwei markante Charakteristika, ag, AEE _, Via rans ae 

die an Hand der vorliegenden vars p ap 3 yout? Vivt 
MeBergebnisse iberprift werden Abb. 9. Theoretische Geschwindigkeitsverteilung im Kreiswirbel. 


* Inzwischen ist aus der deutschen Ubersetzung des Lehrbuches von Kotschin! eine funktionentheo- 
retische Lésung der Gleichung (13) bzw. (16) bekamnt, in der die Reihenglieder von (17) durch Bessel-Funktionen 


zusammengefaBt werden. Man erhalt 


1 1 ay WR 
fir RSA: V=pe(*t Dey 12 AR), 
cw 


[o.e) 
1 Lee A\k 
fir R= A: v= zl Cs) > (2) nea] 
k=0 


1 N, J. Kotschin, I. A. Kibel u. N. W. Rose, Theoretische Hydrodynamik, Bd. 2, S. 320, Berlin 1955. 


“= 


kénnen. Es ist jedoch in keinem Falle auch bei sehr jungen Wirbeln eines der beiden Merkmale, 


wenigstens innerhalb der MeBgenauigkeit, bemerkt worden. Uber den Verlauf des Anstieges im — 


Nullpunkt mége erst im Anschlu8 an das zweite Wirbelmodell diskutiert werden. 


2. Flachenwirbelmodell. Wieder lése sich seitlich vom Zylinder eine Wirbelschicht ab, die 
sich stromabwarts kriimmt und diesmal zu einer Spirale aufrollt. Hat die Spirale den Durch- 
messer 2a erreicht, so reiBe die Wirbelschicht am Zylinder ab, und die von der Spirale bedeckte 
Flache bilde angenahert eine Kreisflache, die mit einer gleichmabigen Dichte von Elementarwirbeln 
belegt sei. Setzt man wieder die Zirkulationsdichte y, die jetzt jedoch die Dimension sec habe, 
und den Kreishalbmesser a, so erhalt man fiir die Wirbelstarke 


a 2x (x — a cos gy)? + (y —a@ sin ¢)? 


ieee | [ae ae da dp (18) 


Sart 


a 72 + a? 


Hig oe | Peete ip (i 7 ae (19) 


4vt 2Qvt 


Unter Benutzung des Stokesschen Satzes ist dann die Geschwindigkeit 


Wife 2 j Ay fie 
p= [wrdr=zl | [race ax Joli gyi) daar (20) 


und nach Einfiihrung der dimensionslosen Koordinaten, wobei analog 


2V4ve 
VO Gry 
zu setzen ist, ergibt sich 
1 48 
V =a | [RAeHH I ((2 RA)dAAR. (21) 
606 
Verwendet man die Integrationsformein 
as 
2e-(44+8) FG2 RA) = | e ' I(Rt)J(At) ede (22) 
0 
und beachtet, daB 
A 
A 
| A J,(At) dA = 7 J(A t) (fiir R entsprechend) , (23) 
0 
so ist 
[o,0) G 
3 2 f a 4 
— | an Eines (24) 


0 
Ww . . . . 40) = 
obei t einen dimensionslosen Integrationsparameter bedeutet. Auch hier ist eine geschlossene 


Integration nicht méglich. Nach Auflésung der Bessel-Funktionen in unendliche Reihen und glied- 


weis ion i ABi 
er Integration ist es zweckmaBig, das Summenschema nach neuen Variablen umzuordnen 
und man erhalt | 


B= Aa Re, 


= : i 
ReS Bs), mit | C= AR, (25) 


n=0 


216 Timme: Uber die Geschwindigkeitsverteilung in Wirbeln Ingenieur-Archiv 


XXV. Band 1957 Timme: Uber die Geschwindigkeitsverteilung in Wirbeln 217 
wobei f(B,C), = 
1 1 B 1 B 1 BS 1 
= | pee ) 
1(y) ar She ay le x) 
CAA ede: 1 B 1 B® 
oles At 5 ot 6 se i 
4- oiieimomsTneMtemis sis; rey a. ¢. Mer (op You ‘elais) f6, shuey so, redone, sluice, “aleve fe 
Gee 1 1 B 1 B 1B 
| aw = 3: és xe : ae = ean eee 
ae oe. i TE Dae) Wel Ine 3! | i (2>b) 
> Bae, oy co (-1)” B™ 
— ni(n+i)! eo me) int © 
foe} lee) 
, = 1)" cn Be 
oe > > Geol Gin om ly | 


bedeuten mége. Auch hier erkennt man sofort, daB fiir 4 = 0 (B = R?, C = 0) die Hamel-Oseen- 
sche Geschwindigkeitsverteilung vorliegt. 

Aus den in Abb. 10 aufgetragenen Kurven ist zu ersehen, daB die so erhaltenen Wirbelprofile 
bei gleichen Parameterwerten A wie im Kreiswirbelmodell immer der Hamel-Oseenschen Geschwin- 
digkeitsverteilung ahnlich bleiben. 

Da nun die gemessenen Wirbel- 
profile ebenfalls der Hamel-Oseen- | ; 
t— V=2/'24 7 (at) (Rt)dt — 


V 


schen Geschwindigkeitsverteilung 
entsprechen, diese aber fiir sehr 
kleine Zeiten aus obenerwahnten 
Griinden sicher falsch ist, laBt sich 
schlieBen, daB man mit der Vor- % 
stellung des Flachenwirbels den 

wahren Verhaltnissen auch bei 

sehr jungen Wirbeln sicher naher 92 


06 


kommt. 

Es galt daher zu untersuchen, p 
welchen EKinflu8 die endliche Aus- j 7 - 5 ee aan 
dehnung des Wirbelkernes gegen- 
ttber dem  Hamel-Oseenschen Vac, Vt _ 2Vbvt ele fisetee 
Punktwirbel auf den Impuls- aoe a7 Ee Ve 


austausch- Vorgang im Innern des Abb. 10. Theoretische Geschwindigkeitsverteilung im Flachenwirbel. 
Wirbels ausibt, d.h. mit anderen 

Worten, wie sich die nach dem Flachenwirbelmodell zu berechnende kinematische Zahigkeit y von 
der Hamel-Oseenschen unterscheidet. Unter Vernachlassigung der Glieder héherer Ordnung erhalt 
man fiir den Geschwindigkeitsanstieg im Nullpunkt 


a? . 
/ i ( ~ an) 9 

= ay era (26) 
und hieraus die kinematische Zahigkeit mit der Anderung der Geschwindigkeitsanstiege zu den 
Zeiten t, und t, mit 


9 
a” 


1 il 
4 At Co, 2.0 a? _ ( Cog, 27 4) ; 
Por ea AT a es 


Man erkennt aus den Klammerausdriicken unter dem Logarithmus, daB der Kernradius nicht 
beliebig groB sein kann, sondern héchstens 


(27) 


_— 


ip 


, 
2 It Cp, 


(28) 


a 
a = 


Dies ist jedoch nur zur Zeit t, = 0 mit konstantem Anstieg bis zum Kernradius méglich (dann ist 
gerade [' = $ ¢d3 =2ma%c,). Zur Darstellung der Funktion » =» (a, Ac,,, At, I’) seien alle 


218 Timme: Uber die Geschwindigkeitsverteilung in Wirbeln Ingenieur-Archiv | 


/ 


itbrigen Werte auf den Wert 1] normiert, und ¢,,,/¢,,, = 0,7 gesetzt,-was etwa den allgemeinen | 
MeBgréBen entspricht. Dann ist 


1 
alc), a\] a\2\|° Q?) 
mney] P(e) 

Aus der graphischen Darstellung Abb. 11 ist zu ersehen, daB bis zu a/a) = 0,7 die Abweichung 
von y kleiner als 2%, bei a/ay = 0,8 Kleiner als 5% bleibt, und erst bei a/ay = 1, ein Wert also, 


der in realer Form nie existieren wird, die kinematische Zahigkeit noch nicht einmal den doppelten 
Wert annimmt. 


normiert auf 


Se ail 


04 


02 
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Abb. 11. Kinematische Zahigkeit im Flachenwirbel als Funktion des Kernhalbmessers a. 


Mit dieser Rechnung konnte gezeigt werden, da die nach Hamel-Oseen berechnete Austausch- 
gréBe fir Impuls in weitem Bereich auch fiir den Flachenwirbel giiltig ist. Und weiterhin, daB 
der mégliche Fehler nur negativ auftritt, d. h. die GréBe y kann nur zu klein, nie zu groB bestimmt 
werden. 
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08 
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Abb. 12. Geschwindigkeitsanstieg im Wirbelmittelpunkt als Funktion der Zeit t, 


AnschlieBend sei noch ein MeBergebnis vermerkt, welches sich nicht ohne weiteres in das 
Schema des Flachenwirbels einordnen ]aBt. Bei mehreren, der Lage nach noch sehr jungen Wirbeln 
ist zweifelsfrei innerhalb der MeBgenauigkeit ein Wachsen des Geschwindigkeitsanstieges im Null- 
punkt gemessen worden. Tragt man wie in Abb. 12 fiir einen als Beispiel herangezogenen Wirbel 
mit / =5 cm?/s, a = lem, »y = 0,01 cem?/s den zeitlichen Verlauf des Geschwindigkeitsanstieges 
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im Nullpunkt fiir die drei besprochenen Wirbelmodelle auf mit 


; Je 20 
(Hamel-Oseen)  ¢,, = Pay ara 
a 25 
(Kreiswirbel) = ; - : re - ane ; (30) 
a 25 
(Flachenwirbel) Cas = = 6 — A im) = 0,8 f — e) 7) 


so sieht man, da8 nur im Falle des Kreiswirbels ein Wachsen des Geschwindigkeitsanstieges erfolgen 
kann. Da jedoch die gemessenen Wirbelprofile die andern Merkmale des Kreiswirbels nicht auf- 
weisen, lat sich vermuten, da die urspriingliche Wirbelstarkeverteilung wohl flichenhaft an- 
geordnet ist, die Wirbelstarke zum Rande des Kreises jedoch zunimmt. Ein derartiges Bild lieBe 
sich unschwer in die Vorstellung von dem zeitlichen Ablauf der Aufrollung der Wirbelschicht ein- 


fiigen (Abb. 13). 
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Abb. 13, Wirbelstarkenverteilung zur Zeit t = 0. ali i 
Vin 
Abb. 14. Normierte maximale Geschwindigkeit V,, und zugehériger Radius R,, 
als Funktion des Kernhalbmessers A. ok 7 4 


0 7 2 


Eine Abschatzung des urspriinglichen Kernhalbmessers, die allerdings nur bei einem allein- 
stehenden Wirbel zulassig ist, sei hier kurz angedeutet: Die normierte maximale Umfangsgeschwin- 
digkeit V,, und der zugehérige Radius R,,, die sich aus den Gleichungen (17) und (25) explizit 
schlecht darstellen lassen, werde durch eine Naherungsfunktion aus den numerischen Beispielen 
Abb. 9 und 10 bestimmt. Beide GréBen sind dann nur noch Funktionen des normierten Kern- 
radius A (Abb. 14). 

Die gemessenen Werte c,,, und r,, normiert man, indem man den unbekannten Wert j4vt=x 
frei wahlt und mit den so gebildeten 

2nV4vt 2% 


rm Vm 
eon, T om pe R,, = penoae 


in die graphische Darstellung Abb. 14 eingeht. Der gewahlte Wert x ist iterativ solange zu ver- 
andern, bis die jeweils gefundenen A-Werte iibereinstimmen. Dann ist der urspriingliche Kern- 
halbmesser 


G=2A. 
Versuchsweise nach diesem Verfahren durchgefiihrte Rechnungen ergaben fiir verschiedene Wirbel 
gréBenordnungsmabig 
a=0,9+ 0,4cm. 
16 


220 Timme: Uber die Geschwindigkeitsverteilung in Wirbeln Ingenieur-Archiv 


Hierzu muB nochmals betont werden, daB es sich nur um eine Abschatzung handelt, da das Ver- | 
fahren nur fiir alleinstehende Wirbel abgeleitet war. So erhielt man z. B. fiir denselben Wirbel zu_ 
den Zeiten t, und t, immer etwas unterschiedliche Werte fiir den urspriinglichen Kernradius. 
Den normierten Kernhalbmesser A kénnte man noch auf andere Weise bestimmen, wenn man 
die kritische GréBe e 
Com ps m 
Koes = Sia (31) 


Coy Tm 
bildet. Es stellt sich heraus, daB K fiir den Hamel-Oseen-Wirbel eine Invariante, fiir den Kreis- 
und Flichenwirbel jedoch eine Funktion von A ist (Abb. 15). Mit dem aus den MeBgroBen c,,,,, 

c,, und r,, berechneten K kénnte man aus der graphischen 
i Darstellung Abb. 15 das A bestimmen, sofern die MeBgenauig- 


i et keit dazu ausreicht. Ist aber jeder Wert nur mit einem 
Fehler von 5°% behaftet, so ist 
| AK Apia Nee etn 
— aoe 20 my I 15° 
und bei einem beispielsweise gemessenen 
i | K = 0,65 + 0,1 = 0,55—0,75 
i 1 erstreckt sich der Fehlerbereich schon tiber alle fiir den 
| Flachenwirbel berechneten K-Werte. 
iz | b) Die Wirkung der Nachbarwirbel auf den Auf- 
| wirbel. Alle bisherigen Betrachtungen gehen von der Vor- 
| 
3 H 4, 
Dm Vn ! . 
Tie aE I 
nae 
a i 
/ marae 77 
/ 2 ee 
f i: Kreiswirbe/ | 1 a 
1 wi | 
7 — Flichenwirbel L 
res A 
955 7 2 Zylinder 1. bis 4 Wirbe/ hinter Zylinder 


Abb. 15. KenngréBe K als Funktion des Abb. 16. Kdrmdnsche WirbelstraBe. Lage der Achsen SS’ 
und PP’. 


Kernhalbmessers A, 


aussetzung aus, daB der zu untersuchende Wirbel, im folgenden mit Aufwirbel bezeichnet, allein 
in einer ruhenden Fliissigkeit existiert. Fiir einen Aufwirbel im Verband einer WirbelstraBe ist 
das nun offensichtlich nicht der Fall, vielmehr wirkt die Umgebung — das sind die Nachbar- 
wirbel — derart auf ihn ein, daB der Wirbel sich mit einer Eigengeschwindigkeit in Richtung 
auf den Zylinder bewegt. Uber den Mechanismus der Induktionswirkung der Nachbarwirbel 
kann man von vornherein keine exakten Angaben machen, da als Liésung der Navier-Stokes- 
schen Differentialgleichung fiir Wirbelbewegungen nur der Spezialfall des Einzelwirbels bekannt ist. 

Ist die Bewegung eines Wirbels in einem gewissen Bereich, dem ,,Wirbelkern“., mit der eines 
starren Kérpers vergleichbar — und nach den bekannten Experimenten mit gefarbten Wirbel- 
kernen erscheint einem diese Annahme plausibel —, so miiBte die Induktion sich nur auf den 
Wirbelmittelpunkt, analog dem Schwerpunktsatz der Mechanik, auswirken und der Wirbelkern 
als ganzes die Geschwindigkeit annehmen, die die Nachbarwirbel an der Stelle des Aufwirbelmittel- 
punktes erzeugen. 

Diese Frage ]aBt sich experimentell entscheiden, wenn man das Geschwindigkeitsprofil des Auf- 
wirbels nicht nur auf der Achse senkrecht zur WirbelstraBe (SS’ in Abb. 16) aufnimmt, sondern 
zusatzlich noch auf der Achse parallel zur WirbelstraSe (PP’ in Abb. 16), durch den wahren Wirbel- 
mittelpunkt. Handelt es sich um die Rotation eines starren Kérpers, so muf in beiden Fallen 
der Geschwindigkeitsanstieg im Mittelpunkt der gleiche sein. Aus Abb. 17 und 18 ist zu entnehmen, 
da diese Annahme offensichtlich nicht zutrifft. Auf der senkrechten Achse ist der gemessene 
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Geschwindigkeitsanstieg um 26% kleiner als auf der Parallelachse, auch ist der Unterschied zwischen 
den beiden mittleren Com Mit 34,5% betrachtlich. Somit kann man auch den innersten Bezirken 
eines Wirbels kaum die Eigenschaften eines starren Kérpers zusprechen. 

In der Potentialtheorie der idealen Fliissigkeiten, die sich von der Eulerschen Gleichung fir 
rotorfreie Strémung ableitet, gilt allgemein das Superpositionsprinzip, nach dem man aus zwei 
partikuléren Integralen der Gleichung durch Addition ein allgemeineres bilden kann. In unserem 
Falle hieBe das, daB sich die resultierende Geschwindigkeit an jedem Ort vektoriell aus den Ge- 
schwindigkeitsanteilen der Nachbarwirbel und des Aufwirbels zusammensetzt. Nun ist aber einmal 


; op [cm/sek] 


iy 


. [em] 
5 4 8 A ey. 2 4 § 


wahrer Wirbelmittelpunkt os 
scheinbarer =» 


O4|- —Geschwindigkertsanstieg 
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ee 0,331sek 
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Abb. 17. Gemessenes Wirbelprofil auf der Achse SS’ senkrecht zur Wirbelstrabe. 
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Abb. 18. Gemessenes Wirbelprofil auf der Achse PP’ parallel zur Wirbelstrafe. 


der untersuchte Bereich des Aufwirbels keineswegs rotorfrei, und zum anderen gilt die Hamel- 
Oseensche Lésung der Navier-Stokes-Gleichung nur bei rotationssymmetrischen Anfangsbedin- 
gungen. Trotzdem soll versuchsweise das Geschwindigkeitsprofil auf den beiden Achsen nach dem 
Superpositionsprinzip berechnet werden mit der Annahme, daf die Konvektionsglieder zahlen- 
maBig keinen allzu groBen Einflu8 ausiiben werden. 

Dafiir ist es zweckmaBig, nach dem Hookerschen Vorgehen! alle Nachbarwirbel durch 
Potentialwirbel zu ersetzen, da in gréBerer Entfernung der Unterschied zwischen Real- und Poten- 
tialwirbel vernachlassigbar klein ist. Dann wird nach v. Karman das Geschwindigkeitspotential 


1S, G. Hooker, Proc. Roy. Soc. A 154 (1936), S. 67. 
16* 
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der WirbelstraBe durch 
Aa, HG 
ie sin —- (2 —%) 
re ae (32) 
sin iE (2 + 2) 


dargestellt. Mit dem Aufwirbel an der Stelle 


I » It 
ie ales (33) 
ist das Geschwindigkeitsfeld 
ah 
d RE a (7) 
(62) . v 
Ae. al a ee ime Th 2x (34) 
i Sin ("] + cos Ga ) 
und auf der senkrechten Achse SS’ mit 
za=ztiy (35) 
die Geschwindigkeit in x-Richtung 
Wo} Ga 
ie l 
pe (36) 


An der Stelle des Aufwirbels wird der Potentialwirbel durch einen Hamel-Oseenschen Punktwirbel 
(1) und (2) ersetztt: 


“Ct : 
=r l I Agena: ( ~ an) 
ree ein 7) — on ) 2ar' Qnr 
7 Co} ca a = r? 
a 4vt 
l an, Hite win OEE lar” i (37) 
Gin Cr] — Gin (F y) 
wobei 
h 
PEE oe (38) 
ist. Wird in (37) y durch r gemaB (38) ersetzt, so folgt 
ip Cig Cr} pee 
es chee (39) 


f 1— 61077} ein (Tr) Gol (TT pane 


Fiir kleine r entwickelt man mit Gin r X r und Coj r ~ 1 + 1/2 r? die Gleichung (39): 


ea(7) 
ell r 
Up 2arbig(T)+ 27 Senda Reo 
=a t6 | 
r Dt 1 ; " 
7 | 40 
1+ 7 30(7)r 8auVt ( ) 
coals 
=f r 


—— 
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In dem Beispiel aus Abb. 17 und 18 waren die dazugehérigen Werte 
h= 1,9 cm J == Siem" |6 (es gia 


4] 
1 =9,1 cm Ayt = 2,5 cm? ee 


berechnet. 


Die mit den Werten berechnete Kurve konnte mit der gemessenen besser in Ubereinstimmung 
gebracht werden, wenn man 


T= 5,713 em" ss (41a) 
setzte. Dann lautet (39) 
a — 1,094 0,913 a3 
% = T—1,738 - Sin (0,69 r) — €oj (0,697) or i Oo) 
und (40) 
= 0.912 
us = sp 0,365 r. (43) 
w[cm/sek] 


gemessen 
——-—berechnet 


-G2 
cosh ach, Bere 
T'=5730m’/sek dvt=2o5cem yal. (7 vce 8 
U sinh (22)- sinty[2%y) ear 
h=19em L-grem 2 =9209 : : 


a a 3 
Abb. 19. Gemessenes und berechnetes Wirbelprofil auf der Achse SS’ senkrecht zur WirbelstraBbe. 


In Abb. 19 sind das berechnete und das gemessene Wirbelprofil fiir die Achse senkrecht zur Wirbel- 
straBe SS’ eingetragen. Nach dem gleichen Verfahren soll nun die Geschwindigkeitsverteilung in 
y-Richtung auf der zur StraBe parallelen Achse PP’ berechnet werden. Gleichung (35) lautet jetzt 


z=x+ia (44) 
und entsprechend (36) 


i a (45) 


Nach zweckmaBiger Umformung ist 


ne: Sof? (7) ge (7 x) = 


; 2 Gin? (T) sin ‘) - cos? (7) 


Wieder wird der Aufwirbel durch einen Realwirbel ersetzt [vgl. (37)]: 


n cor (7) cos x) os 275 
vt 
Or I zt h\ 2 200 oer es (47) 
2 Sint | i )p sn ( i | | cost ( i z) 
mit 
ae (48) 
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Wird x durch (48) ausgedriickt, erhalt man 


22 2 

wo (7 ale " es 

pase T" 4vt 
ee os jh (49) 


2 2 
2 Sin? (7) cos ( i rf | sn*( . " oe 


Geht r—> 0, wird sin (ar) % ar, cos (ar) Y 1 —1/2 a? r* und (49) 


h\ 2 
ry hl eg 
—T aH ( l ies ae jE ) 
Oe (ah\2n2r?  4n?r? | Qar 8avt 
2 cin ( i p 2 
— Gof? & ‘) 
Lr eee 
2 Gint ("Par bane 2a sae (50) 
— Gof? | 
= ie 4 as 
[ein (7) + 1]2ay ee nee, 
(ough oe Ir J 
> 8axt® 
Mit den in (41) und (41a) angegebenen Zahlen lautet (49) 
a 0,942 sin (0,69 r)_ |. 0,913 25 (51) 
0,991 [1 — cos (0,69 r)] + sin? (0,69 r) r 
und (50) 
Vo = — 0,365 r. (52) 


Abb. 20 zeigt das berechnete und das gemessene Wirbelprofil ebenso wie Abb. 19 in bemerkens- 
werter Ubereinstimmung, sofern man die grundsatzlichen Bedenken gegen das Verfahren und die 
Vereinfachungen des Ansatzes beachtet. Setzt doch das Kdérménsche Geschwindigkeitspotential 


v{cm/sek] 


Uf xt. 
5 7-4-3 
/ 
gemessen 
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pw tosh? (22 zh) «sin 4 Bie 
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i 
Abb. 20. Gemessenes und berechnetes Wirbelprofil auf der Achse PP’ parallel zur WirbelstraBe. 


eine beiderseits unendliche lange WirbelstraBe in aequidistanter Anordnung mit konstanter Wirbel- 
starke voraus. In Wirklichkeit ist die StraBe nicht unendlich, wird die StraBenbreite h mit zu- 
nehmendem Abstand vom Zylinder gré®er und sind die Starken der Wirbel nur in einem Bereich 
von +16% (s. S. 212) untereinander gleich. 
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Da das Superpositionsprinzip dem tatsichlichen Mechanismus verhaltnismaBig gut entspricht, 
ist hiermit die wesentlichste MeBgréBe, der Geschwindigkeitsanstieg im Wirbelmittelpunkt, zu 


untersuchen. Nach (40) betragt er auf der senkrechten Achse SS’, auf welcher allgemein die Aus- 
wertung erfolgte, 


duy Sn re if Jf ; I ah 
Rl=*=gere faye ee Z = ts ("}. 23) 
Wird die Gleichung nach yt aufgelist: 
— 1p 
Ue aero aoa a : (54) 
8 x (uy — >—s— 
( Qn(Z+ >) 
so lat sich y aus den MeBwerten zu den Zeiten t, und t, bestimmen [vgl. (4)]: 
Tf 1 il 
y —= — — = - ~ —_——— —— = ——— 
8 a At ae ie 7 [- 
one tr || 8! Sa(Z ery 
, , ‘ag LP SO (] 
a | cae es mit Z’ = ——_— Wea 55 
8 w At up, Ups + (Ug, + Uge) Z’ + 22 rs 2 72 . (59) 


Aus einer zahlenmaBigen Durchrechnung an Hand der vorliegenden Mef&daten ergibt sich der Unter- 
schied zu dem mit (4) berechneten Wert zu 


Tat (5 ! 1 L ni) = + 21%: 


7 a 
es Uo1 Uoe U1 Ue 


Das bedeutet, daB das oben angegebene Auswertungsverfahren mit einem systematischen Fehler 
behaftet ist, der sich aber im Vergleich zur sonstigen MeSgenauigkeit nicht erheblich auswirken 
dirfte. Da das Vorzeichen des Fehlers eindeutig ist, 14Bt sich sagen, die berechnete kinematische 
Zahigkeit wird auf Grund der Induktion der Nachbarwirbel um einen begrenzten Betrag zu gro} 
gemessen. 


7. Zusammenfassung. Es wird iiber eine experimentelle Arbeit aus dem Hermann Fottinger- 
Institut an der Technischen Universitat Berlin-Charlottenburg berichtet, in welcher die Geschwin- 
digkeitsverteilung in Wirbeln einer Kdérmdnschen WirbelstraBe bestimmt wurde. Die Messung 
erfolgte durch Photographieren der markierten Fliissigkeitsoberflache. Unter Zugrundelegung 
des Hamel-Oseenschen Geschwindigkeitsgesetzes fiir den realen Wirbel mit . 


konnte die wirksame kinematische Zahigkeit »y und die Zirkulation [’ berechnet werden. Fiir 
Wirbel im Bereich einer kritischen GréBe wurde die kinematische Zahigkeit y rund 10fach gréBer 
als mit einem Standardviskosimeter bestimmt. Unterhalb der kritischen Gréfe ergab die Zahig- 
keit nur kurz nach der Entstehung der Wirbel einen erhéhten, danach jedoch den normalen lami- 
naren Wert. Als Deutungsméglichkeit der MeBergebnisse wird ein turbulenzartiger Impulsaus- 
tauschvorgang in Ubereinstimmung mit Autoren theoretischer Arbeiten angegeben. 

In einem weiteren Abschnitt wird das Auswertungsverfahren beziiglich des verwandten Wirbel- 
modells und der Induktion der Nachbarwirbel einer kritischen Betrachtung unterzogen. Es konnte 
gezeigt werden, daB das Verfahren vermutlich mit systematischen Fehlern behaftet ist, diese aber 
eroBenordnungsmaBig unterhalb der beobachteten Effekte liegen. 


(Eingegangen am 6. September 1956) 


Anschrift des Verfassers: Dr. A. Timme, Hermann-Fottinger-Institut fur Strémungstechnik, 
Berlin-Charlottenburg, Hardenbergstr. 34. 
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Berichtigung 
zu meiner Arbeit S. 90 dieses Bandes ,,Uber die Biegung 
des orthotropen Plattenstreifens durch Hinzellasten“ 


Von S. Woinowsky-Krieger 
Die erste Gleichung (2) muf richtig 


ae at 
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(Eingegangen am 14. Marz 1957.) 
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Praktische Winke zum Studium der Statik. cranaesen. 
Anwendungen. Rechenkontrollen. Von Dr.-Ing. habil. Ernst Kohl, 0. Professor an der Tech« 
‘nischen Hochschule Braunschweig. Mit 208 Abbildungen. VIII, 224 Seiten Gr.-8°, 1957. 


Ganzleinen DM 18,— 


Inhaltsiibersicht; Die Stabilitatsbedingungen ebener Tragwerke und ihre Anwendung. — Die Aufgabe 
der Statik. — Die Lésung der Gleichgewichtsaufgabe. — EinfluBlinien. — Die Lisung der Formanderungs- 


‘aufgabe.—Querbelastete Balken mit Langskraft. — Statisch unbestimmte Tragwerke. — Stabilitat raum- 


licher Tragwerke. — Sachverzeichnis. 


. Bei dieser Neuerscheinung handelt es sich nicht um ein Lehrbuch der Statik im iiblichen Sinne. Die kurz 


1 


gehaltenen »Praktischen Winke*t wenden sich in erster Linie an den Studierenden, geben aber auch dem 
im Beruf stehenden Statiker mit ihren zahlreichen, an Beispielen durchgefiihrten Rechenkontrollen wert- 


-volle Hinweise. Die Grundlagen sind klar herausgestellt. Ausfiihrlich behandelte Beispiele mit ver- 


schiedenen Lésungsméglichkeiten und deren Vor- und Nachteile machen mit den Anwendungen vertraut. 
Statisch unbestimmte Tragwerke sind nach den KraftgréSen- und Ror Ae wna inns be- 
rechnet. Die zu lésenden Aufgaben werden von verschiedenen Seiten aus beleuchtet. 


Statik der Tragwerke. Von Dr.-Ing. habil. Walther Kaufmann, 0. Professor an 
der Technischen Hochschule Munchen. (Handbibliothek fiir Bauingenieure. Ein Hand- und Nach- 
schlagebuch fir Studium und Praxis. Begr. von Robert Otzen. Neue Reihe herausgegeben von Prof. 
Dr.-Ing. Dr.-Ing. E. h. Ferd. Schleicher.) Vierte erginzte und verbesserte Auflage. Mit 367 Abbildungen. 
VII, 327 Seiten Gr. -8°, 1957. ; Ganzleinen DM 31,50 
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: Inhaltsibersicht: Erster Abschnitt: Allgemeine Grundlagen. Begriff hd Aufgabe der Statik. Die 


auBeren Krafte. Die inneren Krafte. Die Gleichgewichtsbedingungen des starren Kérpers und das Prinzip 
der virtuellen Verriickungen.  Statisch bestimmte und statisch unbestimmte Systeme. Die HinfluBlinie. 


Die Grundgleichungen der Statik des stabférmigen Trigers. Die Grundlagen der Fachwerktheorie, — 


Zweiter Abschnitt; Momente, Quer- und Normalkrafte an statisch hestimmten Stabwerken. Der ein- 
fache Balken. Freithager, Balken mit tberkragenden Enden und Gerbertriger. Dreigelenkbogen und ver- 
wandte Systeme. — Dritter Abschnitt: Ermittlung der Spannkrafte statisch bestimmter Fachwerke.- 
Statische Verfahren fiir das ebene Fachwerk. Die kinematische Methode. Raumliche Fachwerke. — 
Vierter Abschnitt: Die elastischen Formianderungen. Das Prinzip der virtuellen Verriickungen. Die 
Satze von der Gegenseitigkeit der elastischen Formanderungen. Der Castiglianosche Satz vom Differential- 
quotienten der Formanderungsarbeit. Die Biegungslinie. Vollstandige Darstellung der Formanderung 
ebener Systeme. — Fiinfter Abschnitt: Theorie der statisch unbestimmten Systeme. Einfiihrung. Das 
statisch bestimmte Hauptsystem. Die Elastizitatsgleichungen fiir die statisch wnbestimmten GréBen. 
Auflésung der allgemeinen Elastizitatsgleichungen. LHinfiihrung statisch unbestimmter Hauptsysteme. 
Reduktionssatz. Aufstellung von Elastizitatsgleichungen mit nur einer Unbekannten. Einfiihrung von 
FormanderungsgréBen als Unbekannte. Berechnung statisch unbestimmter Stabwerke durch Momenten- 
ausgleich. — Sechster Abschnitt: Statisch unbestimmte Tragwerke. Der durchlaufende Triger. Der 
beiderseits eingespannte Trager. Der Trager auf elastischer Unterlage. Rahmen. Bogentrager. Durch 
einen einfachen Balken versteifte Gelenkbégen und Ketten. Durch einen iiber drei Offnungen laufenden 
Vollwandtrager versteifte Kette. Dreifach statisch unbestimmter Bogen iiber drei Offnungen. Briicken- 
tragerrost tit drillsteifen Haupttragern. Einfiibrung i in die Theorie zweiter Ordnung (Verformungstheorie). 
— Sachverzeichnis. 
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Mathematische Statistik. Von Dr. B. L. van der Waerden, Professor der Mathe- 
matik an der Universitit Ziirich, (Grundlehren der mathematischen Wissenschaften, Band 87.) Mit 
39 Textfiguren und 13 Zahlentafeln. IX, 360 Seiten Gr.-8°. 1957. DM 46,—; Ganzleinen DM 49,60 


Inhaltsiibersicht: Allgemeine Grundlagen. — Wahrscheinlichkeiten und Haufigkeiten. — Mathe- ~ 


matische Hilfsmittel. — Empirische Bestimmung von Verteilungsfunktionen, Mittelwerten und Streu- 
ungen. — Fourier-Integrale und Grenzwertsitze, — GauSsche Fehlertheorie und Students Test. — 


Die Methode der Kleinsten Quadrate. — Schatzung unbekannter Konstanten. — Auswertung von beob- . 
achteten Haufigkeiten. — Bio-Auswertung. — Priifung yon Hypothesen durch Tests, — Anordnungs- — 


tests. — Korrelation. — Tafeln. — Leitfaden. 


Das vorliegende Buch ist aus einer langjahrigen Beschaftigung mit den praktischen Anwendungen hervor- 
gegangen. Seit der Studentenzeit wurden dem Verfasser immer wieder statistische Fragen von Volks- 
-wirtschaftlern, Medizinern, Physiologen, Biologen und Ingenieuren vorgelegt. Durch Nachdenken und 
Literaturstudium lernte er immer bessere Methoden kennen. Diese Methoden werden im vorliegenden 
Buch begriindet und auf méglichst lehrreiche Beispiele aus’ den Natur- und Sozialwissenschaften an- 
gewandt. Der Verfasser hofft, dem Leser so manche Irrwege zu ersparen. ie 


Die Beispiele sind nicht aus der Theorie heraus konstruiert, sondern der Praxis entnommen; daher waren 


bei manchen Beispielen ausfiihrliche Erlauterungen notwendig. Die mathematischen Grundbegriffe sind 
so kurz wie méglich, aber doch verstandlich dargestellt. Manchmal waren langere theoretische Aus- 
fiihrungen notwendig, aber wo immer moglich, I fiir schwierige Beweise ayf gute exes Lehr- 
biicher verwiesen. } 


Wahrscheinlichkeitstheorie. Von Dr. Hans Richter, 0. Professor fiir Mathes . 


matische Statistik und Wirtschaftsmathematik an der Universitat Miinchen. (Die Grundlehren der 


mathematischen Wissenschaften, Band 86.) Mit 14 Textabbildungen, XII, 435 Seiten Gr.-8°. 1956. 
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